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Hochenergiestreuung im Regge-Limes 

Diese Arbeit umfafit zwci Untcrsuchungcn, die bcidc damit zu tun habcn, Rcggc- 
Theorie von QCD ausgehend zu verstehen. Gegenstand der ersten ist, wie das 
Proton an das Odderon koppelt, ein Reggeon, das keine Ladung tragt und unge- 
radc untcr Ladungskonjugation ist. Die zwcitc bctrifFt Unitaritatskorrckturcn zur 
BFKL-Gleichung, die das Pomeron beschreibt, ein Reggeon mit den Quanten- 
zahlen des Vakuums. 

Im ersten Teil dieser Arbeit wird der Oddcron-Beitrag zur elastischcn Strcuung von 
Protonen an Protonen und Antiprotonen bcrechnet. Um den Einflul3 der Proton- 
struktur auf die Odderon-Proton-Kopplung zu untersuchen, wird ein geometrisches 
Protonenmodell konstruiert. Durch Vergleich mit experimentellen Daten wird die 
durchschnittliche Gr6i3e eines Diquark-Clusters im Proton ermittelt. Zwei weit- 
ere Odderon-Proton-Impaktfaktoren aus der Literatur werden durch Vergleich mit 
dem Experiment getestet. 

Der zweite Teil enthalt die Berechnung von vier-Pomcron-Vertizes, die im Rahmen 
der Generalised Leading Logarithmic Approximation auftreten. Diese Naherung 
dient zur Unitarisierung von Streuamplituden, die den Austausch des BFKL- 
Pomerons beschreiben. Eine Anzahl grundlegcndcr Funktionen, aus dencn solche 
Vcrtizes bestehen, werden systematisch im Impuls- und Ortsraum bchandclt, imd 
ihrc Transformationseigcnschaften unter konformen Transformationcn im transvcr- 
salen Ortsraum werden hcrgelcitet. Die Vertizes werden als eine Zahl von Intc- 
gralen ausgedriickt und auf die Form einer konformen Vierpunktfunktion gebracht. 



High energy scattering in the Regge limit 

This thesis comprises two investigations, both connected with the attempt to un- 
derstand Regge theory in the framework of QCD. The first is about how the 
odderon, a Reggeon carrying no charge which is odd under charge conjugation, 
couples to the proton. The second concerns unitarity corrections to the BFKL 
equation which describes the pomeron, a Reggeon with the quantum numbers of 
the vacuum. 

In the first part of this thesis, the odderon exchange contribution to clastic proton- 
proton and proton-antiproton scattering is computed. A geometrical transverse- 
space model of the proton is constructed to investigate the influence of a possible 
diquark cluster in the proton on the odderon-proton coupling. The average size 
of this cluster is determined by comparison with experimental data. Furthermore, 
the validity of two odderon-proton impact factors from the literature is tested. 
The second part consists in the derivation of four-ponicron vertices. These ver- 
tices occur in the Generalised Leading Logarithmic Approximation used to uni- 
tarise scattering amplitudes describing the exchange of a perturbative pomeron. A 
number of basic functions of which such pomeron vertices are composed is treated 
systematically in momentum and position space. Their conformal transformation 
properties in impact parameter space are derived. The vertices are expressed as 
a number of integrals and cast into the form of conformally invariant four-point 
functions. 
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Introduction 



This introductory chapter gives an overview of the theories and physics issues 
surrounding the two investigations contained in this thesis. The basic notions 
of high energy scattering and of descriptions of the strong interaction are briefly 
portrayed. This introduction is intended to give some background and put the 
topics of this thesis into their context, wliicli is why technical details, notably of 
the BFKL equation, are left for later chapters. 



1.1 Background and Motivation 



Quantum Chromo dynamics is now well established as the fundamental theory of the strong inter- 
action. The quark model which it is based on provides the basis of our understanding of matter. The 
description of the quarks' interaction by the gauge theory that is quantum chromodynamics has been 
successfully applied to a very wide range of problems. This interaction varies smoothly between the 
two extremes of confinement, which makes it impossible to isolate quarks, and asymptotic freedom, 
which makes the quarks inside scattering hadrons behave like free particles at small distances. 

Considering how successful and how satisfactory in view of the aim to understand nature from 
first principles QCD is, it is not surprising that it has eclipsed the attempts to describe the strong 
interaction which preceded it. One of those is Regge theory [Re59, Re60, Col77]. It is based on Lorentz 
invariance, unitarity and analyticity of the scattering matrix. Developed before quarks and gluons were 
recognised as the fundamental degrees of freedom of the strong interaction, it is a theory of hadrons. 
It has been quite successful in describing hadronic scattering processes in the so-called Regge limit, 
for asymptotically large energies and momentum transfers of the order of a hadronic mass scale. 

In Regge theory, the exchange of particles in a scattering process is described by the singularities of 

the scattering amplitude in the complex angular momentum plane, the so-called Regge poles and cuts. 
By way of crossing symmetry, they can be related to the masses and spins of existing hadrons. By that 
token, every hadron is a Regge particle, or Reggeon. However, some Reggeons do not correspond to 
any known hadrons. One such is the pomeron which is the pole with the largest real part and therefore 
dominant at asymptotically high energies. It carries the quantum numbers of the vacuum, and the 
simplest model for it consists of two gluons in a colour singlet state [Low75, Nus75]. Another is its 
charge-conjugation-odd partner, the odderon [LN73]. There are some indications that both pomeron 
and odderon might be related to glueballs. 

To this day, concepts of Regge theory are used widely both in phenomenology, eg [JarSO, Jar82, 
DL84, GLN90, DL92] and theory, eg [Gri68, BT76, BS78, FL93, Li94a, FK95, Pes97]. The existence of 
the pomeron is universally accepted, and it is the topic of a wide range of phenomenological work, for 
instance [DL89, DL92, Bra..92, BP95, FL98, KLOO, MNOO, KKLOl, St02]. The odderon, by contrast, 
has never been measured experimentally beyond doubt and consequently is still a contentious topic. 
A brief overview of odderon physics will be given in the following; for an in-depth review, I refer 
the reader to [Ew02]. There are few processes in which the odderon is the dominant or indeed only 
contribution. One such process is the diffractive photo- or electroproduction of rjc or other heavy 
pseudoscalar mesons at HERA [CKMS97, BBCVOl, EIKS98]. But the cross-sections are estimated 
to be too small to measurable in the foreseeable future. Another odderon- dominated process is the 
production of light pseudoscalar or tensor mesons, the cross sections for which are expected to be 
higher [KN98, Ry98, BDD..99, BDDNOO]. But it has not been observed experimentally [OlOl, GolOl, 
Be02, Be02a]. This discrepancy might be due to difficulties with the non-perturbative models which 
are required to describe effects at such low scales. 

However, the processes on which the search for the odderon has concentrated for a long time are 
proton-proton and proton-antiproton scattering, a fact which the large number of publications on the 
topic attests, for instance [DL84, GNL85, GLN90, Za89, LR90, DGJ92, BN99, LTOO, DGMPOO]. These 
processes provide the hitherto only experimental evidence for the odderon. The elastic differential cross 
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sections of these two processes differ qualitatively: In proton-proton scattering, it shows a dip structure 
at a squared momentum transfer t of about t = l GeV^, while in proton- antiproton scattering, it only 
flattens off at this point. Odderon exchange can account for this difference since it carries odd charge 
conjugation parity and hence contributes with different signs to both processes. What is more, the 
odderon is the dominant C = — 1 contribution arising from Regge theory. It is expected to have a 
Regge intercept of slightly below 1 , while C-odd meson trajectories typically have intercepts of around 
0.5. Since the scattering amplitude of a Reggeon exchange with intercept ao is proportional to s"o+" *, 
the latter contribute much less at high energies. This difference in the differential cross sections has 
been measured at the Intersecting Storage Rings at CERN [Bo.. 74, Nag.. 79, AS80, Erh..85, Bre..85]. 
The experimental evidence is marred by the fact that the statistics for pp scattering is low and that 
the difference only shows in a few data points (see Section 2.1). A good description of the data in 
the framework of Regge theory was given by Donnachie and Landshoff [DL84] . It is characterised by 
fair accuracy and considerable predictivity, having only a small number of free parameters. A few 
years later, it was extended [DL86] to higher-energy proton-antiproton data measured at the CERN 
SPS [Ber..86, Ber..87]. Proton-proton data at even higher energy from the Tevatron [Am. .90] provide 
no further check on the fit since they do not extend far enough in —t to show the dip. In the Donnachie- 
Landshoff fit, the odderon-exchange contribution produces the dip in proton-proton and its absence 
in proton-antiproton differential cross sections. Another description of pp and pp scattering is due to 
Gauron, Nicolescu and Leader [GLN90]. It uses the so-called maximal odderon which corresponds to 
a different type of Regge singularity from other odderon contributions. As in the Donnachie-Landshoff 
fit, the odderon contribution is instrumental in the description of the dip region. However, the maximal 
odderon itself already shows a rich structure with a dip. One finding of the work presented in this 
thesis is that each odderon contribution is compatible only with its own fit, that exchanging odderon 
contributions between fits is not possible (see Section 4.3.3). But no successful description of the data 
without an odderon has been found so far. A measurable quantity which is more sensitive to odderon 
exchange than the cross sections is the ratio of the real to imaginary part of the forward scattering 
amplitude. But also here no evidence for the odderon has been found [Aug. .93]. 

If the thin experimental evidence for the odderon were all there is to it, one might have discarded 
it as a misguided concept. But other concepts of Regge theory remain very valid today. Furthermore, 
the odderon can in fact be derived from perturbative QCD. It is described by the Bartels-Kwiecinski- 
Praszalowicz (BKP) equation [Ba80, KP80, Jar80]. Odderon exchange amounts to a simultaneous 
exchange of at least three gluons in a C = — 1 state. Since such an exchange is clearly possible in QCD, a 
failure to find the odderon at all would be a heavy blow to QCD. Odderon research in QCD has received 
a boost from the discovery that the perturbative odderon is equivalent to an integrable model, the XXX 
Heisenberg model with spin zero [FK95, Ko95]. A prime interest in odderon physics has always been 
the determination of its Regge intercept, see for instance [GLN93, GLN94, AB94, AB97, Br98, Br98a]. 
Since the discovery of an additional conserved charge of the odderon [Li94] , work has concentrated on 
finding the spectrum of the associated operator [Ko97, Ko98, MW95, VLOl, DKMOl]. Recently, this 
spectrum has been determined both by Korchemsky et al. [KKM02, DKKM02] and by de Vega and 
Lipatov [VL02] . To date, there are still unresolved differences between the results of the two groups 
which centre on the question which eigenvalues of the operator are physical. Separately, explicit 
solutions of the BKP equation have been found by two groups. One, by Janik and Wosiek [WJ97, 
JW99], has an intercept slightly below one (w 0.96). It is of the general form derived by Lipatov 
and collaborators [Li90, GLN91, Li93], in which an analytic auxiliary function is left open. The other 
solution, by Bartels, Lipatov and Vacca [BLVOO], has an intercept of exactly one. This means that its 
contribution to forward scattering would not decrease with energy. The Bartels-Lipatov- Vacca (BLV) 
solution was contentious at first since it is not of the general form found by Lipatov, at least not with 
an analytic auxiliary function. But now it looks like being accepted as legitimate by the community. 
This is not just of academic theoretical interest, nor just a question about a slightly different intercept. 
Odderon solutions conforming to Lipatov's general form cannot couple to some scattering particles, 
for instance a photon fluctuating into an r}c. The BLV odderon, however, can. Therefore, with the 
advent of this new solution a number of new processes in which the odderon may play a role becomes 
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available. 

A solution of the BKP equation, which describes the odderon's propagation, is only one ingredient 
for calculating the odderon contribution to a given process. Besides, impact factors describing the 
odderon's coupling to scattering particles arc required. Pom,eron impact factors have always aroused 
some interest [As.. 90, KSG95, Go96] and have just been computed in next-to-leading order [CROO, 
BCGK02, Gie02, FIK02]. Less work has been done on odderon impact factors [FK79, LR90, DL91, 
CKMS97]. While some impact factors can be calculated pcrturbatively (for instance the odderon- 77c 
impact factor [CKMS97]), impact factors for proton scattering require non-perturbative ansatze most 
of which have never been confronted with experimental data. As part of this work, the odderon- 
proton impact factors of Levin and Ryskin [LR90] and Kwiccinski ct al. [FK79, CKMS97] were used to 
calculate the elastic differential cross section for proton-proton scattering. The odderon was described 
simply as three gluons in a C = — 1 colour singlet state. The non-odderon contributions were taken 
from the Donnachie-Landshoff fit [DL84] . This allows to fix the coupling constant to be used with these 
impact factors. My results imply that predictions of the rjc production cross section by Kwiecihski 
and two other group using his impact factors [EIKS98, BBCVOl] have to be revised downwards, see 
Section 3.4. 

The main interest of the first part of this thesis (Chapters 2, 3 and 4) is however a different 
one. It has been known for some time that the structure of the proton has considerable effect on 
the oddcron-proton impact factor. In the extreme case in which two of the proton's valence quarks 
are clustered in a point-like diquark, the odderon-proton coupling even vanishes [Za89]. The odderon 
contribution for a diquark cluster of finite size in the proton was computed non-perturbatively by 
Riiter and Dosch [RD96, Rue97]. They used it to determine the ratio of the real to imaginary part 
of the forward scattering amplitude for proton-proton and proton- antiproton scattering, respectively. 
The perturbative odderon contribution to these processes for a finite-sized diquark is presented in this 
thesis, in Chapter 2. By computing the differential cross section in the dip region for proton-proton 
scattering, the average size of the diquark cluster is obtained from experimental data (Chapters 3 
and 4). The geometrical nature of the concept of a diquark cluster necessitates the use of a position- 
space approach to high-energy scattering. I choose the one developed by Nachtmann [Na91], which is 
otherwise used mainly with non-perturbative models, notably the Model of the Stochastic Vacuum, 
see for instance [Do87, DS88, Si88, DFK94, RD96, BN99, BDD..99, SSDP02, SSDP02a]. My result 
of a small diquark size will benefit non-perturbative calculations: Many of them treat the proton as 
a colour dipole [DFK94, BN99, SS02], which is legitimate since soft gluons cannot resolve a small 
diquark. 

Important though many results of Regge theory remain in present-day high-energy physics, our 
understanding of the Regge limit from the first principles of QCD is still far from complete. This is 
because high-energy hadronic processes tend to involve large parton densities and are dominated by 
soft scales. Our knowledge of non-perturbative QCD is not yet sufficient for a rigorous treatment of the 
vast majority of high-energy processes. However, some processes are dominated by one hard scale and 
therefore allow the application of perturbation theory. One such is the scattering of heavy quarkonia. 
The process itself is not experimentally feasible. But since photons tend to fluctuate into qq states 
in high energy scattering, the scattering of highly virtual photons is expected to provide a suitable 
substitute. This is predicted to be measurable [BDL96, BDEL97, BHS97, BHS97a] at a future e+e" 
coUider like the planned Next Linear Collider [NLC, Ku..96]. A situation similar to the process just 
mentioned occurs in the presence of hard forward jets in deep inelastic or in hadronic scattering, called 
forward jets and Mueller-Navelet jets, respectively [Mue91, MN87]. They have been investigated in 
depth by a variety of groups, for example [BDL92, KMS92, Ta92, BBRK92, NPR96, BDD..96, ElOO, 
MuOO, ADF..01, ASD..01, EIMOl, BCVOl]. Another process involving a hard scale is deep inelastic 
electron-proton scattering at small values of the Bjorken scaling variable x. In this case the hard scale 
is provided by the virtual photon which mediates the interaction on the electron's side. 

Even in the presence of a hard scale, perturbation theory in high-energy QCD is not trivial. Due 
to the large energies available, the phase space of intermediate gluons becomes very large, oc log s. In 
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QCD, this logarithmic factor can compensate the smallness of the couphng constant, with the effect 
that more complicated Feynman graphs are not always suppressed. Therefore an infinite number of 
graphs have to be resummed to obtain a consistent leading-order result. This is called the Leading 
Logarithmic Approximation. The resummation was performed for the first time in QCD by Balitsky 
and Lipatov [BL78], using earlier work by Lipatov and collaborators [KLF77]. The result of the 
resummation is described by the well-known BFKL equation. The corresponding Reggeon exchange is 
called the BFKL, perturbative or hard pomeron. It represents a bound state of two reggeiscd gluons 
and was found to correspond to a Regge cut in the complex angular momentum plane. A rcggcised 
gluon is a colour octet state composed of several QCD gluons [Li76]. It was discovered in 1985 that 
the BFKL equation is conformally invariant in transverse position space [Li86]. This was a finding 
of great importance since it allowed the solving of the equation in the non-forward direction. It also 
led to a novel interpretation of BFKL evolution as a two-dimensional conformally invariant quantum 
mechanics in which the complexified angular momentum is the energy-like and rapidity the time-like 
variable. 

Since the invention of the BFKL pomeron it has become clear that it violates a basic principle of 
field theory: unitarity. The Froissart-Martin theorem, derived from this principle, states that hadronic 
cross sections can rise at most with the squared logarithm of energy. The growth of the BFKL- 
pomeron-exchange cross section is power-like, with an exponent of In 2Afc/7i' ~ 0.5. But the BFKL 
pomeron represents only a leading-order approximation. Higher-order corrections tame this power-like 
growth and restore unitarity. 

A complete resummation of higher-order graphs relevant for pomeron exchange is even today 
thought to be prohibitively difiicult. One has to be content with adding specific classes of graphs 
which are chosen to restore compliance with the Froissart-Martin Theorem and hence unitarity. For 
restoring the unitarity of the overall scattering matrix one has to add ladder graphs with more than 
two reggeised gluons in the t channel. These additional terms are called unitarity corrections. The 
amplitudes (or propagators) of these larger bound states of reggeised gluons are described by the 
Bartels-Kwiecinski-Praszalowicz (BKP) equations [Ba80, KP80, Jar80]. They can be greatly simpli- 
fied by treating only the large limit. Then a nearest-neighbour interaction takes the place of the 
n-particle interaction between the reggeised gluons. This amounts to a one-dimensional solid, and 
was found to be equivalent to the XXX Heisenberg model with spin zero, a completely integrable 
model [Li94, Li94a, FK95]. This limited approach of demanding unitarity only for the overall process 
is followed by Lipatov [Li94a]. 

A more ambitious and more complete approach is due to Bartels [Ba91, Ba93, Ba93a]. It demands 
unitarity not just in the main process but also in all possible subprocesses. This requires the inclusion 
of graphs in which the number of reggeised gluons in the t channel is not constant. Starting from two 
reggeised gluons coupling to the scattering particles, amplitudes of more than two reggeised gluons 
emerge. It was found that the three- and five-gluon amplitudes are in fact just superpositions of 
amplitudes with fewer gluons [BW95, Ew98, BE99], a phenomenon known as reggeisation. The four- 
gluon amplitude also contains a reggeising part [BW95] . But besides that a vertex transforming two to 
four gluons occurs, which gives rise to a new form of four-gluon amplitude. The conformal invariance 
of this new vertex was shown in [BLW96]. The six-gluon amplitude, investigated in [Ew98, BE99], 
consists of a completely reggeising part, a partially reggeising part, and a new 2^6 vertex. The 
completely reggeising part is a superposition of two-gluon amplitudes, the partially reggeising part 
a superposition of four-gluon amplitudes. The 2^6 vertex again gives rise to a new amplitude. 
If all transitions to a larger number of reggeised gluons were conformally invariant, these unitarity 
corrections would amount to a conformal field theory in 2 -|- 1 dimensions. In analogy to the conformal 
quantum mechanics suggested by BFKL, rapidity would serve as the time-like variable and the complex 
angular momentum as the energy-like variable. 

Recent work [Pes97, Ko99, JK99, GKK02] suggests an even more intimate connection between 
string theory /conformal field theory and unitarity corrections. Some of this work involves pomerons 
as the fundamental degrees of freedom and pomeron vertices. Vertices of BFKL pomerons can be 
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obtained by projection of the reggeised-gluon transition vertices. The simplest of these, transforming 
one to two pomerons, has been investigated in [Lot96, Pes97, Ko99, Br99, BRV02] and found to have 
the form of a conformal three-point function. There is some confusion about the existence of higher 
vertices transforming one to n > 2 pomerons. Of course higher vertices could be constructed by 
iterating the 1 — > 2 vertex. But such composed vertices would not be local in rapidity since the 
intermediate states of pomerons would have to be allowed to evolve. Pcschanski [Pes97] has given 
an interpolated expression for local 1 — > n pomeron vertices which, like the 1^2 vertex, allows an 
interpretation as string- theoretic Shapiro- Virasoro amplitudes. Since then, Braun and Vacca [BV99] 
have proven that such local vertices do not exist in the dipole approach of Mueller [Mue94, Mue95]. 
This approach provides a derivation of the BFKL equation which is much easier than the original one 
involving gluon ladders. However, it is not thought to be equivalent to the gluon-ladder approach to 
all orders, not least because it does not allow the exchange of an odderon. While Peschanski also starts 
out from dipoles, the expression he generalises is much closer to the gluon-ladder form of the result 
obtained by Lotter [Lot96] . 

Some of the questions about higher pomeron vertices are addressed in the second part of this work. 
Chapter 5 consists of preparing the groundwork for treating reggeised-gluon and pomeron vertices in 
a systematic way. It deals with a function G in terms of which the 2^4 and 2-^-6 reggeised gluon 
vertices can be expressed, and its components. The properties of G under conformal transformations 
prove the conformal invariance of the 2^6 vertex and also of higher vertices if they have an analogous 
form as suggested in [EwOl]. This represents an important step towards proving the existence of a 
conformal field theory of unitarity corrections. In Chapter 6 the 1^3 pomeron vertex is derived in 
the gluon-ladder approach and cast into the form of a conformal four-point function. Its existence does 
not contradict the result of Braun and Vacca but rather reveals a discrepancy between the gluon-ladder 
and the dipole approach. Their equivalence extends only to the 1^2 pomeron vertex but not to 
higher orders. 
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1.2 High-energy scattering 

When a number of particles get close to each other, they interact, possibly break up and fly off in 
directions different from that of their initial trajectories. This is called a scattering process. To simplify 
its description, one assumes that the scattering particles approach each other from a macroscopically 
large distance, interact in a microscopically small region, and that the products of the interaction 
move away again to a macroscopically large distance. Before and after the interaction, the particles 
are thought of as so far apart that they are effectively unaware of each other. The scattering operator 
S is defined as the operator which transforms the initial multi-particle state in the infinite past to the 
final multi-particle state in the infinitely remote future. 

hm \t) = \f)=S\i)=S lim \t) 

t—*+oo t— > — oo 

Since the time-evolution operator is linear, so is the scattering operator. In a suitable base, it can 
therefore be written as a matrix, the scattering matrix, with elements denoted by Sji. More often 
than iS, one uses the T matrix. It is obtained from iS by subtracting the unit matrix and removing 
some kinematic factors. It therefore can be thought of as describing the effect of the interaction, since 
a scattering matrix equal to unity amounts to no interaction. 

Here Pi and Pf are the total four-momenta of all particles in the initial and final states, respectively. 
The T matrix element is also called the scattering amplitude and denoted by T. It is not usually 

expressed as a function of the initial and final states. Rather, kinematic variables which characterise 
the states are used as arguments of T. Easily the most important kinematic variable is s. It is the 
square of the centre-of-mass energy, or equivalently the square of the invariant mass of the scattering 
system. 

s = Pi=Pj 

In "high-energy" scattering, one treats the limit s — > c», or in experimental reality the case where s is 
significantly larger than the squares of all the particle masses, s 3> Vn. 

Let us now turn to the case of two particles scattering and two (possibly different) particles emerging 
after the scattering process. This situation is displayed graphically in Figure 1.1. The diagram can be 
read in two ways: from left to right or from top to bottom. The left-to-right view corresponds to the 
reaction 

a + b — > c + d. {s channel) 
The top-to-bottom view corresponds to the reaction 

a + c — > b + d. {t channel) 
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b d 
Figure 1.1: A two-particle to two-particle scattering process. 



b and c become antiparticles since they now propagate backwards in time. This process has the centre- 

of-mass energy (pa ~ Pc)"^ =■ t. Looking at the the f-channcl process instead of the original s-channel 
reaction is called "crossing" . There is a third kinematic variable u and a corresponding third process 
in which it is the squared centre-of-mass energy: 

a + d — > b + c {u channel). 

These three kinematic variables are very important in any of these processes. They are called Man- 
dclstam variables. They are not independent but add up to the sum of the mass squares of the four 
particles. Here is their definition: 

\2 



•« = {Pa +Pb) 
t = {Pa -Pcf 
U={Pa- Pdf 



(1.1) 



In the s-channel process, s is the squared centre-of-mass energy, t is the squared four-momentum 
transfer: If one assumes the particle a to emerge as particle c, the four-momentum p^ — Pc has been 
transferred to b (which becomes d). This is also the negative square of the invariant mass of the 
(virtual) particle mediating the interaction. One can show that in s-channel processes, t is always 
negative. After crossing, the variables exchange their roles. In the f-channel process, t is the squared 
centre-of-mass energy and s is the squared four-momentum transfer. Now s is negative and t is positive. 
(It has to be taken into account that in the crossed process, the momentum of c is —pc-) 

Interestingly, the same scattering amplitude applies in all crossed processes. It just has to be 
continued analytically to the region in which the Mandelstam variables have the right sign for the 
desired channel. This is called crossing symmetry. 

Let us now turn from kinematics to calculating observables for s-channel processes. The observables 
of most universal interest are cross sections. The most widely used differential cross section is the 
derivative of the cross section with respect to t. It is easily computed from the scattering amplitude: 

da _ T{s,tf 
H ~ 16 7rs2 ■ 



(1.2) 



The total cross section could of course be obtained by integrating over the differential cross section. 
But there is a simpler way which makes it possible to obtain it directly from the scattering amplitude. 
This is by using the optical theorem. It is a special case of the Cutkosky rules and relates the forward 
elastic scattering amplitude to the total cross section. In the high-energy limit, the optical theorem 
has the following form: 



2lmT{s,0) = 2s(Jtot{s). 
This allows to compute the total cross section. 



(1.3) 



1.3. The standard model, the strong interaction and QCD 
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1.3 The standard model, the strong interaction and QCD 

The standard model is a deseription of the three interactions relevant in elementary particle physics — 
the electromagnetic, weak and strong interaction. This description is based on quantum field theory. 

The electromagnetic and weak interactions are jointly described by the Weinberg-Salam model. It 
is basically a gauge theory with four massless gauge bosons and a scalar particle, the Higgs. By fixing 
the vacuum expectation value of the Higgs field, one induces spontaneous symmetry breaking. As a 
consequence, a gauge theory with three massive {W^, Z^) and one massless (7) gauge bosons emerges. 
Even though the Higgs boson has not yet been detected, this model is very successful. 

The strong interaction is described by quantum chromodynamics (QCD). Its foundation is the 
quark model, the idea (supported by experiment) that hadrons are composed of smaller point-like 
particles which cannot be isolated, the quarks. Quarks have two quantum numbers peculiar to them: 
flavour (which accounts for different types of hadrons) and colour (which is the gauge charge of QCD). 
There are three different colours, and to date six flavours of quarks have been found. Observable 
particles are superpositions of systems of quarks of all three colours, which renders them colour- 
neutral, "white". The gauge bosons of QCD are called gluons. They are massless and come in the 
eight existing colour-anticolour combinations, excluding white. 

Making predictions on the basis of QCD is much harder than in QED, for two reasons: One, gluons 
couple to each other. There is a three-gluon and a four-ghion vertex in QCD. Two, QCD exhibits con- 
finement. The most obvious effect of this property is that colour-carrying particles cannot be isolated. 
Its deeper reason is that the strong coupling constant is large at small momenta. (Correspondingly, it 
is small at high momenta, which is called asymptotic freedom.) Taken together, these two points mean 
that gluons have a tendency to split up into more gluons (or quark pairs) and that more complicated 
Feynman graphs are not always significantly suppressed. 

Problems concerning the strong interaction at low energies cannot be solved with perturbative 
QCD. Additional models or lattice calculations are required. Perturbation theory does not work at 
low energies. Perturbative calculations of the evolution of the strong coupling constant lead to a 
divergence at the scale Aqcd ~ 250 MeV. The expression for this evolution probably starts to lose its 
meaning below two to four times that value. Non-perturbative QCD is not the topic of this thesis, but 
the limitations of perturbation theory have to be kept in mind in any QCD calculation. 

1.4 Before QCD: Regge theory 

1.4.1 Regge theory and the scattering matrix 

Before the development of QCD as a quantum field theory of the strong interaction, attempts were made 
to derive the properties and behaviour of hadrons from reasonable assumptions about the scattering 
matrix. The resulting theory is called Regge theory after its inventor, T. Regge [Re59, Re60]. (See 
[C0I77] for a thorough treatment.) 

The assumptions about S are the following: 

1. <S is Lorentz-invariant. We have already made this assumption in the first section when we wrote 
the scattering amplitude as a function of the Lorentz invariants s and t. 

2. <S is unitary, <S<S^ = <S^<S = 1. This is the conservation of probability, ie the probability for a 
specific initial state to end up in any final state at all must be one. 

3. <S depends analytically on Lorentz invariants. Unitarity takes precedence over this, ie singularities 
which are necessary for unitarity are allowed. 

The last point requires some explanation. It may seem strange that unitarity requires singularities 
in the scattering amplitude. The reason lies in the optical theorem (1.3) which relates the forward 
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elastic scattering amplitude with the total cross section. The Cutkosky rules from which it is derived 
are a consequence of unitarity of the <S matrix. By way of the optical theorem, contributions to the 
total cross section (from any reaction) relate to contributions to the imaginary part of the forward 
elastic amplitude. Therefore one can learn something about the scattering amplitude by looking at all 
the final states which might arise from two scattering particles. 

Two particles can merge into one only at discrete energies. They can create two different particles 
only at energies higher than the combined mass of the lightest pair of particles allowed by quantum 
number conservation. Therefore there are intervals on the real s axis where the total cross section, 
and therefore the imaginary part of the forward scattering amplitude, is zero. This allows to apply 
the Schwarz reflection principle from the theory of complex functions: When a function is analytic in 
a given domain and real on a straight line contained in that domain, then it takes values which are 
complex conjugates of each other at points which are mirror reflections with respect to that line. In 
this context this means that 

T{s\f) =T{s,t)* 

in its domain of analyticity. (The straight line is one of the intervals on the real s axis where T is real, 
so the mirror image of s is s*.) 



single-particle resonance poles 




two-particle production cut 



three-particle production cut 



Res 



two-particle production cut 



three-particle production cut 



Figure 1.2: Singularity structure of the scattering amplitude in the s plane. It is symmetric to a line 
parallel to the imaginary axis, but not necessarily to the axis itself. 

Since an analytic function has a unique power expansion in its domain of analyticity, T would have 
to be real on the whole real s axis if it were to be analytic there. But that would mean atot = for 
all energies, which is clearly not usually true. Therefore T has to be singular at large parts of the real 
axis. The non-zero multiple-particle production amplitudes are represented by branch cuts (since they 
extend all the way to infinity from the n-particle-production threshold), while single particle production 
is represented by isolated singularities, usually poles. It can be shown from crossing symmetry and 
the relation between the Mandelstam variables that T has a structure symmetric with respect to a 
parallel to the imaginary axis, so that the singularity structure shown in Figure 1.2 emerges. 



1.4.2 Complex angular momenta 

Let us consider two-particle to two-particle scattering. It is well known from quantum mechanics 
that one can express the scattering amplitude in terms of the spherical harmonics with zero angular 

momentum, Yj o, which are Legendrc polynomials in the cosine of the scattering angle, cos^. This is 
called partial wave expansion. In the s channel, cos 9 = 1 + 2s/t and one can write: 



oc 

T(s, t) = ^(2/ + 1) ai{s) Pi{l + 2s /t) , 

1=0 



(1.4) 



1.5. Pomeron and odderon 
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where a; are the coefficients of the expansion, the partial wave amphtudes. This expression can be 
rewritten according to Cauchy's integral theorem as a contour integral in the complex I plane over a 
function with singularities at I E N>o. The value of the integral is the sum of the residues at these 
poles, which are constructed to equal the terms in the sum in (1.4). This is called Sommerfeld- Watson 
transform. One obtains: 

T(s, t)= idl (2/ + 1) ail, s) P(l, 1 + 2s/t) . (1.5) 

/ sin TTt 

The integration contour extends to infinity and circles the positive real axis clockwise. The argument 
I of a and P is no longer written as an index to indicate that it is a now complex variable and that 
these functions have been analytically continued over the whole complex I plane. 

The analytic continuation of the partial wave amplitude a is not unique. Adding a function with 
zeros at all non-zero integers to it does not change the integral. But such a function has to be a 
sine with a period of two divided by an integer. Such a sine function would rise exponentially on the 
imaginary axis. So there is a unique a{l,s) which rises more slowly than exp(7r|Z|) as I approaches 
infinity. However, some contributions to the partial wave amplitude oscillate as (—1)'. This would be 
analytically continued to a cosine, which also rises faster than our limit. Since adding a sine function 
cannot correct that on both the positive and negative imaginary axis, a different route is chosen. The 
analytic continuation is done for even and odd angular momenta separately, and a{l, s) is replaced by 



77=±1 ^ 

T] is called the signature of the partial wave. 

To investigate the form of the scattering amplitude in the high energy (Regge) limit, one transforms 
the integration contour so that the contour integral vanishes in this limit. What remains are the 
residues of the integrand's poles the contour swept over during the deformation. These are called 
Rcggc poles. Using an asymptotic formula for the Legendre polynomials, one obtains for each of the 
poles a contribution 

T(s,i)ocs"W, (1.6) 

where a{t) is the pole's position in the complex I plane depending on t. In view of this form, one can 

say that at high energies the scattering amplitude is dominated by the Regge pole with the largest 
real part. The expression (1.6) is usually interpreted as the exchange of a "Reggeon", an object with 
angular momentum equal to a{t). a{t) is called its Regge trajectory. 

To gain some understanding of Regge trajectories, one continues them to the region t > 0. It was 
found that mesons plotted in a diagram of angular momentum against square mass are arranged in 
straight lines. Based on this fact it is usually assumed in Regge theory that Regge trajectories are 
linear functions also for t < 0, that is 

a{t) — + at . 

ao is called a Reggeon's Regge intercept, a' its Regge slope. From (1.6) one derives that the differential 
cross section for a process in which a single Reggeon is exchanged is 

— OC s2ao+2a't-2 

dt 



1.5 Pomeron and odderon 



Both pomeron and odderon are not particles but "objects" which can be exchanged in scattering 
processes. They are well described by Regge theory, and hence are Reggeons. This section will give a 
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brief introduction to them and to how to understand them in terms of QCD. For a detailed treatment 
see [FR97] or [DDLN02]. The state of the art of odderon research has been reviewed very recently 
in [Ew02]. 

The existence of the pomcron is universally accepted nowadays. It is based on the following 
path of reasoning; The Pomcranchuk theorem [Pom58] states that total cross sections of scattering 
processes in which charged particles are exchanged have to vanish asymptotically for high energy. Here, 
"charge" can mean any kind of charge, including electric charge, isospin, colour charge etc. . However, 
one observes experimentally that at the highest energies for which measurements are available, total 
cross sections continue to rise slowly. This made Pomeranchuk postulate an object, the pomeron, 
which carries no charge whatsoever and therefore has the quantum numbers of the vacuum. Pomeron 
exchange can therefore account for the continuing rise of total cross sections. 

The existence of the pomeron is well-founded now. Pomeron exchange was found to play a signifi- 
cant role also at lower, more accessible energies. (The original argument for the pomcron's existence 
is hampered by the fact that no experimentally measured energy can reasonably be claimed to be 
"asymptotically high", however high it may be compared to other experiments.) 

The odderon is the partner of the pomeron with odd parity and charge conjugation quantum 
number, P = C = —1. All its other quantum numbers are those of the vacuum. There is much less 
experimental evidence for its existence than for the pomeron. The best evidence comes from elastic 
proton-proton and proton- antiproton scattering and will be presented in Section 2.1. Nobody has been 
able to describe the difference between the differential cross sections in these two processes without 
an odderon-exchange contribution, but it is not reassuring that the odderon contributions used by 
different authors are not compatible (see Section 4.3.3). 
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Figure 1.3: Pomeron and odderon in terms of QCD gluons. Even though the same Feynman graphs 
can contribute to both pomeron and odderon, the colour structure is of course different. 

Even though pomeron and odderon were initially described by Regge theory, they can of course be 
understood in the framework of QCD as well. Both are represented as the exchange of a bound state 
of several gluons. In lowest order, the pomeron consists of two t-channel gluons, and the odderon of 
three. Figure 1.3 presents this diagrammatically. Further legitimacy is lent to this representation by 
the fact that some known glueballs lie on the t > Q part of the Regge trajectories of pomeron and 
odderon. 

The description of the odderon in the QCD framework provides the most important rationale for 
investigating it. There may be some doubt about whether the odderon is observed experimentally. 
But there is no doubt that the C = — 1 multi-gluon exchanges which make up odderon exchange do 
exist in QCD. If their existence was to be disproved experimentally, it would have serious implications 
for QCD, which is our best candidate to date for a theory of the strong interaction. 



1.6. BFKL and the violation of unitarity 
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Figure 1.4: The reggeised gluon is a sum of ladders of which the vertical parts are themselves reggeised 
gluons. The rungs are ordinary QCD gluons. The vertices are sums of QCD vertices which take into 
account rungs crossing each other. 

1.6 BFKL and the violation of unitarity 

As described in the previous section, pomeron exchange amounts to the exchange of a bound state 
of gluons. The pomeron is more easily described in terms of reggeised gluons [Li76]. These are gluon 
ladders not unlike the pomeron itself, but in a colour octet state, ie with the quantum numbers of a 
gluon. They are graphically represented as gluon lines with slanted lines through them, as shown in 
Figure 1.4. 

Sums of such ladder diagrams have an important property: At high energies, the phase space for 
the additional rungs (oc logs) compensates the additional couplings, so that all the ladder terms have 
the same order of magnitude. This means that in a consistent calculation, ladders with any number 
of rungs have to be summed up. In addition, there are genuinely higher-order graphs in which there 
are more t-channel reggeised gluons or the rung gluons split up into quark pairs. Summing up all 
leading-order ladder graphs is called the Leading Logarithmic Approximation and results in a sum of 
the following form: 

y^ctn {as logs)" . 

n 

In the case of the pomeron, the result of the resummation satisfies the BFKL (Balitsky-Fadin- 

Kuraev-Lipatov) equation, an integral equation [KLF77, BL78]. It is presented in some detail in 
Section 5.2.6. The BFKL equation has an important property; It is invariant under global conformal 
transformations in transverse position space. Only this made it possible eventually to solve it in the 
non-forward direction, ie for t 0. This also makes it possible to sec the BFKL pomeron as a state of 
two reggeised gluons evolving according to conformally invariant two-dimensional quantum mechanics. 
Rapidity serves as the time- like variable because the rung gluons are strongly ordered in rapidity.^ 

The BFKL equation has one drawback, however: The scattering matrix for BFKL pomeron ex- 
change is not unitary. The growth of the associated cross section with the centre-of-mass energy is 
power-like rather than logarithmic, which violates the Froissart-Martin unitarity bound. This can be 
remedied by adding a minimal set of (genuinely) higher-order terms. The result is then called the 
Generalised Leading Logarithmic Approximation (GLLA). It involves including some terms and re- 
jecting others. The selected terms are those which have the strongest energy dependence and therefore 
dominate at high energies. 

There are two approaches to restoring unitarity. One restricts itself to restoring unitarity to the 
overall scattering matrix [Ba80, KP80, Li94a, FK95, Ko95]. This requires additional Feynman graphs 
in which an arbitrary number of reggeised gluons (but at least two) are exchanged in the t channel. 

^This is because everything else is suppressed by the propagators. The pomeron has to bridge tlic large difference 
in rapidity between the scattering particles. Because the gluon propagators punish large momenta (reggeised gluon 
propagators axe similaj: to QCD gluon propagators, with an additional factor), the situation in which the momenta of 
the gluon rungs evolve monotonously between the opposing scattering particles is preferred. 
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The more ambitious approach initiated by Bartels seeks to restore unitarity also to the scattering 
matrices of subprocesses obtained by cutting up Feynman graphs [Ba80, Ba91, Ba93, Ba93a]. This 
requires not only more reggeised gluons in the t channel, but also vertices which change the number 
of t-channel gluons. This approach is sometimes called the Extended Generalised Leading Logarith- 
mic Approximation (EGLLA). Analogously to BFKL and non-extended GLLA as two-dimensional 
quantum mechanics, this can be interpreted as a two-dimensional quantum field theory, again with 
rapidity as the time-like variable. Its fundamental degrees of freedom are reggeised gluons which in- 
teract through the number-changing vertices. If one could prove that all elements of this field theory 
are conformally invariant, it would be a two-dimensional conformal field theory, about which much 
is known from the side of mathematical physics. By solving the theory, one could gain much better 
understanding of higli-oncrgy scattering. 

BFKL and its generalisations have traditionally concentrated on the scattering of heavy quarkonia 
(ideally) rcsp. of highly virtual photons (in reality). In onium-onium scattering, the pomeron couples 
to a quark loop. It has been shown [BW95, Ew98, BE99] that the coupling of several reggeised gluons 
in a colour singlet state to the quark loop can be described by a superposition of two-reggeised-gluon 
amplitudes. Therefore in this process the initial condition for the evolution in rapidity can be assumed 
to be a two-gluon amplitude. In the EGLLA approach, this amplitude can give rise to amplitudes 
containing a larger number of reggeised gluons. The basis of this approach is a tower of coupled 
integral equations describing these amplitudes. 

The simplest way in which the number of t-channel gluons can change is by way of a phenomenon 
known as reggeisation. This refers to the fact that the reggeised gluon is a superposition of graphs 
which themselves contain reggeised gluons, sec Figure 1.4. This allows a gluon to "split up" into 
several. There is no vertex involved, it only constitutes a regrouping of the QCD gluons contained in 
the reggeised gluon. The resulting n-reggeised-gluon amplitude can then be written as a superposition 
of two-gluon amplitudes. The EGLLA amplitudes with an odd number of gluons investigated so far 
(three and five gluons) were found to reggeise completely [BW95, Ew98, BE99]. 

The amplitudes with an even number of reggeised gluons are more complex. So far only the four- 
and six-gluon amplitude have been investigated. The four-gluon amplitude contains both a reggcising 
part and a new amplitude which cannot be expressed as a superposition of amplitudes with fewer 
reggeised gluons [BW95]. The new part arises from the initial two-gluon amplitude by way of a 2 ^ 4 
reggeised gluon vertex. The structure of the four-gluon amplitude can be represented graphically 
(diagrams from [BE99]): 

\ ^ 



AA 



+ 




The "splitting up" of the reggeised gluons in the left diagram represents the reggeisation. The ellipse 
in the right diagram is the 2^4 vertex. The rectangle at the top of both diagrams represents the 
quark loop — with attached photon lines — to which the pomeron couples. 

The 2 — > 4 reggeised gluon vertex has since been thoroughly investigated. It was found to be 
conformally invariant [BLW96] . It was projected on pomeron wave functions to obtain a 1 ^ 2 pomeron 
vertex [Lot96, Ko99]. This pomeron vertex was found to allow an interpretation as a Shapiro- Virasoro 
amplitude from string theory [Pes97]. The reggeising part of the four-gluon amplitude can also be 
projected onto two pomerons. However, it is sub-leading in the large Nc approximation [BRV02]. 

The six-gluon amplitude also contains a reggeising part and a completely new part. In addition, 
there is a partly reggeising part which is a superposition of irreducible four-gluon amplitudes (ie those 
resulting from the 2 — > 4 gluon vertex). The six-gluon amplitude can like the four-gluon amplitude 
be projected onto pomerons, and a 1 ^ 3 pomeron vertex can be obtained. For the irreducible and 
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the partly reggeising part, this will be done in Chapter 6. This vertex is remarkable for being local 
in rapidity. By contrast, if two 1 — > 2 pomeron vertices occur after each other, the intermediate four- 
gluon state evolves in rapidity, and as a result the transition from one to three pomerons takes some 
interval on the rapidity axis. 

Whether a local 1^3 pomeron vertex exists has boon contentious. Braun and Vacca [BV99] have 
shown that 1 — > n pomeron vertices for n > 2 do not exist in Mueller's dipole approach. This approach 
is a method of deriving the BFKL equation which is much simpler than the original one involving 
gluon ladders [Mue94, Mue95]. Mueller considers a scenario in which scattering colour dipoles (onia) 
recm-sively split up into more dipoles which finally interact by exchanging a single QCD gluon. The 
dipole splitting leads to the same expressions as the superposition of gluon ladders in the traditional 
derivation of the BFKL equation. However, the two approaches are not equivalent to all orders. For 
instance, there is no odderon in the dipole approach: The squared amplitude for single-gluon exchange 
cannot describe the exchange of a C = — 1 object. The existence of the 1^3 gluon vertices derived 
in Chapter 6 in the gluon ladder approach constitutes a further discrepancy. 
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The perturbative odderon in pp and 

pp scattering 



This part is primarily concerned with the influence of the proton structure 
on the odderon-proton coupling. To that end, elastic proton-proton and proton- 
antiproton scattering is investigated. A geometric model of the proton in transverse 
position space is constructed. The odderon-proton coupling depending on the 
size of a possible diquark cluster in the proton is computed. By adding different 
Reggeon-exchange contributions from a fit by Donnachie and Landshoff to the 
perturbative odderon contribution based on the proton model, it is compared to 
experimental data and the average size of a diquark cluster in the proton is fixed 
within limits. 

The same method — substituting a perturbative odderon contribution within 

the framework of the Donnachie-Landshoff fit - is used to test the validity of two 
other models for the odderon-proton coupling, one by Levin and Ryskin and one 
by Pukugita and Kwiecinski. The coupling constant to be used with these models, 
which is the main parameter in these calculations, is obtained from experiment. 

Lastly, a fit for proton-proton and proton-antiproton elastic scattering by Gau- 
ron et al. is compared to the Donnachie-Landshoff fit. In particular, the respective 
odderon contributions of the two fits are exchanged, which allows to make a state- 
ment on the universality, or otherwise, of the odderon-exchange contributions. 



Chapter 2 

The perturbative odderon 
contribution to elastic pp scattering 



2.1 The odderon in elastic pp and pp scattering 

The odderon plays an important part in elastic proton-proton and proton-antiproton scattering. As 
one can see from Figure 2.1, the differential cross section of the two processes differs qualitatively. In 
elastic proton-antiproton scattering, the differential cross section falls continuously with rising squared 
momentum transfer —t. At one point, the curve flattens, but it keeps on falling. The differential cross 
section for proton-proton scattering, on the other hand, shows a slight dip at the same point in —t. 
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Figure 2.1: The qualitative difference between elastic pp and pp scattering. The plot shows the 
differential cross section for ^/s = 53 GeV. Data from [Bre..85]. 



This is obviously an interference effect. It comes from different contributions to the total scattering 
amplitude interfering constructively in one case (proton-antiproton) and destructively in the other. 
Since the dominant contribution to elastic scattering at high energies is pomeron exchange (C = H-l), 
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the contribution responsible for the difference must be odd under charge conjugation (and hence couple 
to proton and antiproton with different signs). 

Odderon exchange is a C = — 1 contribution. It is the common consensus that it is the odderon 
which is responsible for the difFcrcncc between the proton-proton and proton-antiproton clastic dif- 
ferential cross sections. This belief is based not so much on rigorous proof but on the fact that all 
authors who have successfully described this state of affairs include a significant odderon contribution. 
Disturbingly, though, the different odderon contributions seem to work only in the framework of one 
specific fit, see Section 4.3.3. 

Another point of caution concerns the experimental data, or more precisely the lack of data. 53 GeV 
is the only energy at which both proton-proton and proton-antiproton elastic differential cross sections 
have been measured. The difference between the two sets of data is between 3 a and 4 a at one point, 
and 2 cr or less at a few others. This is all the evidence there is for this effect, and the best experimental 
evidence for the existence of the odderon. 

The relevant scale for investigating the dip region is the squared momentum transfer at the dip, 
—tKi 1.4 GeV^. We regard this as sufficiently large for perturbation theory to be applicable. All the 
same it would be desirable to supplement out investigation with a non-perturbative calculation. For the 
time being, the perturbative calculation to be presented is the first step to completely understanding 
the phenomenon of the dip. 

2.2 Computing the odderon-proton coupling using a 
geometric model for the proton 

2.2.1 The geometric model 

In our geometric model, the proton is assumed to be composed of three quarks arranged in a symmetric 
star as depicted in Figure 2.2. The angle a between the two lower quarks is a free parameter of the 
model. For small a the two quarks at the bottom form a diquark cluster of the size rfdiquark- We will 
therefore also call this model the "quark-diquark model" . a = corresponds to an exactly point-like 
diquark. 




Quark 



diquark 



Figure 2.2: The quark-diquark model for the proton 
In addition to the geometry just described, a wave function is ascribed to the proton: 

m = \ll^^M-RVs'), (2.1) 

where R is the distance from the centre of the proton. 5 is a parameter determining the size of the 
proton. In calculations of the electromagnetic form factor it was determined to be 0.96 fm [Pau99]. 
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The size relevant for the strong interaction ought not to be worlds away from that. Non-perturbative 
QCD calculations suggest the range 0.7. . .0.9 fm [DNPWOl, Rue97]. We use S = 0.8 fm. 

The mean diquark size in transverse space, ie the typical distance between the two quarks of the 
diquark, can be easily computed from the geometry and the wave function. It is^ 



An important reason for using this model is a common practice in non-perturbative calculations: 
Often the proton is treated as a colour dipole. The diquark plays the part corresponding to that of 
the antiquark in a meson, which has the same colour structure. This practice is justified if two of 
the proton's quarks are clustered in a diquark so small that a soft ghion cannot resolve them. Using 
this model in a perturbative calculation allows us to determine the diquark size and thus to test that 
hypothesis. 

2.2.2 High energy scattering in position space 

Performing calculations with the geometric model presented in the previous section requires treat- 
ing high energy scattering in position space. This is done using a framework developed by Nacht- 
mann [Na91, Na96]. It was created for evaluating non-perturbative models but is as suitable for a per- 
turbative calculation. Due to its complexity, this section can give only a brief sketch of Nachtmann's 
formalism. For a more elaborate and didactic introduction and applications in non-perturbative QCD, 
see [Na96]. 

Quark-qucirk scattering 

One starts out by treating the elastic scattering of two quarks (even though isolated quarks do not 
exist in QCD). Later we will advance to colourless quark clusters, ie hadrons. The first step is to relate 
the S matrix to a four-point function in position space. The ^-matrix element can be transformed 
into a Green function in position space with the help of the LSZ reduction formalism: 



Here the pi are the momenta of the incoming and outgoing quarks, f{x) = {i^'^d^ — m)ip{x) and 
is the wave function renormalisation. 

The four point function {T'tp{x3)tp{xi)%l){xi)%[){x2)) contained in the expression on the right hand 
side of Equation 2.3 can be expressed as a functional integral over the quark and the gluon fields: 



{Ti;{x3)^{xi)i}{xi)i){x2)) = / Vi,V^VBi;{x3)^{xi)^{xi)i){x2)eM-iStn\\ciCTy\, (2.4) 



where S'fuu qcd is the full QCD action. Since the action is quadratic in the fermion fields, the fermion 
integration is Gaussian. Therefore the quark fields can be integrated out, which gives a determinant 

-'^This is the diquark size in transverse space. If one assumes the geometric star model to be the projection of a 
three-star rotated in three-dimensional space, one can calculate the mean three-dimensional diquaxk size by multiplying 
(2.2) with v^, which is the average of the function | sin | which occurs in the projection. 




(2.2) 




X ( '^u{ps)f{Xi)u{pi)f{XA)f{xr)u{px)f{X2)u{p2) ) ■ 



(2.3) 




While more accessible to the imagination, this three-dimensional diquark size is of limited use since many calculations 
which depend on the smallness of the diquark are done in transverse space. 



2.2. Calculation with a geometric model 
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of the Dirac operator and several quark propagators: 

( Til){x3Mx4)tp{xi)il){x2) ) = jvB det[-i{i^D - m)] ■ (2.5) 

• [Sf{x3,xi;B) Sf{x4, X2;B) + Sf{x3, X2; B) Sf{x4, xi;B)] exp[-iS'pure qcd] • 

SF{xi, Xj;B) are the quark propagators in the external colour potential B^. The remaining functional 
integration is to be performed only over the gluon fields with the pure QCD action in the exponential, 
ie the one describing pure gluo dynamics. 

In leading order the determinant of the Dirac operator can be set to one. The two terms with 
propagators correspond to f-channel and u-channel exchange, respectively. Here we apply the eikonal 
approximation, ie we assume that the four momentum transfer is small compared to the total energy. 
Then the u-channel term can be neglected. 

After plugging (2.5) into (2.3), the propagators can be wrapped up in scattering matrix elements 
for "scattering" one quark in an external gluon field by following the above derivation in reverse. 

(PspT'IpiPi) ~ V / 2^65(^3, Pi ;B)5(p4,P2;B)exp[-iSpu,e qcd] (2.6) 

S{pi,pj;B) is the .S-matrix element for scattering a quark with momentum pj to one with momentum 
Pi in an external colour field Bf^. 



A quark in a background field 



Our next goal is to find an expression for the S'-matrix element S(pi,pj:B). One can show [Na91] 
that the quark scattering matrix elements in an external field can be expressed as a generalised WKB 
expression 

1 



S{pi,pj;B) = u{pi)^^u{pj)P exp 



-^9 j^'Bp 



dx" 



1 + 



Pi 



(2.7) 



where B is the gauge potential which takes values in sm(3). The path-ordered integral is taken along 
the classical path F. For high-energy scattering the relevant paths are the following light-like paths: 



Fi = (a;°,6/2,x^ = a;°) and F2 = (a;°, -6/2,a;^ = -a;°) . 



(2.8) 



They correspond to quarks moving in opposite directions along the x^ axis at the speed of light, h is 
the impact parameter. Let Vj be the phases picked up by the quarks along these paths: 



P exp 



-^9 



I B^(z) dz^ 



It can be derived [Na96] that the quark moving in positive x^ direction has the wave function 

i^i{x) 



Vi e-'f-' 'u(pi) + 



(2.9) 



(2.10) 



An analogous expression is valid for the other quark. From them the S-matrix clement (2.7) can be 
calculated. This can be plugged into the quark-quark scattering matrix element (2.6). The ensuing 
calculations lead to rather lengthy formulas which will not be given here. They are presented very 
well in [Na96]. The S'-matrix element (2.6) is simplified by performing the path-ordered integrations 
from (2.9). After again taking the high energy limit, the result reads: 



where {Express^on) = J VB Expre 



(2.11) 



ssion . 
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The non-Abelian phases Vj still depend on the impact parameter b and the gluon field B. The 
are the colour indices of the quarks and q — pi ~ ps is the transverse momentum transfer. (Small 

bold-face letters denote vectors in transverse space.) Even though the quark momenta are written 
explicitly here, the ^-matrix element depends only on s = {pi +^2)^ and t = {pi — Ps)^- 




Figure 2.3: The paths in a colour neutral three-quark cluster 



Nucleon-nucleon scattering 

In order to come to the nucleon-nucleon scattering amplitude we first consider the scattering of two 
groups of three quarks moving on parallel light-like world lines of the form 

Tl{xo,^ + x1,x^ =x% T^{xo,-^+x^,x^ = -x% a =1,2, 3. (2.12) 

The xf are the positions of the quarks in cluster i relative to its centre. In order to ensure that 
these quark clusters asymptotically form colour singlet states all colours are parallel-transported in 
the remote past and future to the centre of the cluster and contracted there antisymmetrically. This 
leads to the following S'-matrix element for the scattering of colour-neutral clusters [DFK94]: 

X (V})„„, (V?)^^, (V?)^^, e^^fs'ye,,. (V^)^^, (V^^)^^, (V^)^^, e,,^,.') • 

The lower Greek indices are the colour indices of the quarks running from 1 to 3. Because now the 
positions rather than the momenta of the quarks are given, the Fourier integral from (2.11) is missing. 
The non-Abelian phase factors V" are defined as in (2.9) with the U-shaped integration paths Fj in 
Figure 2.3. The Zi denote the wave function renormalisation for the clusters which in lowest order can 
be set to one. 

We will later use not the S'-matrix element but the reduced scattering amplitude J. They are 
related in the following way: 

J{b,{xn) = S{b,{xt})-l. (2.14) 



2.2. Calculation with a geometric model 
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(a) (b) (c) 



Figure 2.4: The three ways in which the odderon can couple to the proton in the lowest order. 

2.2.3 Application to perturbative odderon exchange 

In our perturbative calculation we expand the Vf in g. Up to order g^ Equation 2.9 becomes 

- ^ g'%KiBfA-''^''^''''U + 0{g') , (2.15) 

where 13^ ^ are the coefficients of the expansion of the path integral in the Gell-Mann matrices r^: 

I^Jz^ B^{z) = B^.t'^ . (2.16) 

In general wc would have to take into account path ordering in (2.15). However, the odderon 
contribution which we want to calculate has a symmetric colour structure. This is because the 7 
matrices which are the spinor space part of the vertices are antisymmetric. This already gives a minus 
sign under charge conjugation. Therefore the colour structure has to be symmetric to give overall 
oddness under charge conjugation. Due to the symmetric colour structure permuting the ^ changes 
nothing and we may discard the path ordering. 

The next step is to expand the S'-matrix clement (2.13) in g (using the expansion (2.15)) and to 
extract the terms which represent the lowest order odderon contribution. To this order the odderon is 
just composed of three gluons in a C = — 1 state. Therefore the relevant terms are of order g^. Three 
of the six factors g come from one hadron, three; from the other. Otherwise some of the gluons would 
contribute to the self-energy of a hadron rather than to the interaction between them. As we have 
already argued, the terms are symmetric in colour space. 

The expansion of (2.13) is not given explicitly here because it is too complex. Even the odderon 
involves more than a hundred terms so that they shall not be given in full. It helps a lot to classify 
them according to the ways the odderon couples to the proton at either end. This will be done in 
the next section on the basis of a Feynman graph description of the odderon exchange. This is the 
description which is easiest to implement in the numerical calculations which will be necessary. The 
derivation of an equivalent formulation from Equation 2.13 is described in Appendix A.l. It is rather 
involved and not of direct relevance for the results of our work. 

2.2.4 Colour structure and combinatorics 

Let us first consider one end of the odderon exchange, ie one half of the expression in the brackets in 
Equation 2.13. Three gluons can couple to the proton in the three ways displayed in Figure 2.4. Type 
(a) corresponds to the expression which contains the third-order term of one of the V" and the S^fs of 
the two others (cf. Equation 2.15). Type (b) represents one in which there is the second-order term, 
first-order term and 5 symbol of one V", respectively. In type (c) the linear terms in g of all three 
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Yf appear. Since three gluons couple to each hadron, the sum of the orders of the contributing terms 
from the expansion of the three V" has to be three. 

The colour tensors differ according to the type of the coupling but arc independent of the permu- 
tations, that is which gluon couples to which quark. They are calculated as follows: 

i^) = ^afff {t^t" T° )aa' <5/3/3' ^77' £a'l3'i' = £0/37 ea'/37 (r°T° r° 

= 25ac.' (tV^'t^")„„, = 2Tr(rV'='r^") = Tr([r'=, r'^V^" + {t',t''}t'") 

(b) C^^ = ea/37 (r'^r'^ )aa' T^0' ^77' ea'18'7' = ea/37 ea'/3'7 (r^'r'^ )aa' T^i3' 

= -Tr(T'=T'='T'=") = -\ if"'"" - \ d""'"" 

(c) Cf"" = ea/37 T^„/ Tjij, r^y ea/37 = - d""'"" (numerically) 

As has been mentioned before, only the symmetric colour structure constants d"" " are relevant for the 
odderon contribution because the structure in spinor space is antisymmetric. We obtain the following 
colour tensors for leading-order odderon exchange: 



(2.17) 



Having completed the investigation of the coupling types, we can now classify the Feynman graphs 
according to the types of the couplings at each hadron. Combining three types of coupling at either end 



(a) 


/~1CC c 

^a.. odd — 




(b) 


/ // 

r^cc c 

^b, odd — 




(c) 


r~icc' c" 

*^c, odd — 





Si 



-{E^H> •:^^3> -:£5a> 

(aa) (ab) (ac) 



fbbll mi) fbcl (ccl 



Figure 2.5: Seven exemplary Feynman graphs representing the seven graph types. 

yields six possibilities, ignoring which type occurs at which hadron. However, two (b)-type couplings 
can be combined in two ways: The pair of gluons coupling to the same quark may or may not be 
the same at both ends. Hence we obtain seven types of Feynman graphs. One representative of each 
type is displayed in Figure 2.5. This classification was already used by Riiter in non-perturbative 
calculations [Rue97]. 
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The colour structure of the odderon gives an overall colour factor of d""'"' d'^"'" = In addition, 
the prefactors of the symmetric structure constants in Equation 2.17 contribute different factors to 
different graph types. 

Another typc-dcpcndcnt factor comes from combinatorics: To facilitate our numerical calculations, 
we assume the gluons to be distinguishable. To compensate this overcounting, the result has to be 
divided by the number of permutations of gluons which couple to the same quarks. For instance, 
the contribution of the graphs of type (bbl) will be overcounted by a factor of two because for each 
graph an equivalent one will be calculated in which the two first gluons are exchanged (see Figure 2.5). 
By contrast, for type (bb2) the gluons are distinguished by which quarks they couple to and no 
overcounting occurs. 

Tabic 2.1 shows the colour and combinatorial factors for each graph type. The total type-specific 
prcfactor Ctype is arbitrarily defined as one for graph type (aa). That means that an overall factor of 
i i = ^ has to be multiplied with Ctypo to get the right factor for a specific graph type. Also given 
in Table 2.1 is the number of graphs belonging to each type. They add up to the number of different 
Feynman graphs, 165. 

A more rigorous derivation of the typc-dcpcndcnt prefactors is given in Section A.l in Appendix A. 
They arise from the coefficients in the expansion (2.15) and from summing up equivalent terms in the 
expansion of the six-gluon correlator (2.13). 



Type 


Colour factor 


Comb, factor 


Ctype 


# of graphs 


(aa) 


1 

4 


1 

6 


1 (by definition) 


9 


(ab) 


1 
8 


1 

2 


3 
2 


36 


(ac) 


1 

4 


1 


6 


6 


(bbl) 


1 

16 


1 

2 


3 
4 


36 


(bb2) 


1 

16 


1 


3 
2 


36 


(be) 


1 

8 


1 


-3 


36 


(cc) 


1 
4 


1 


6 


6 



Table 2.1: Prefactors of the different graph types and number of Feynman graphs belonging to each 
type. See text for explanation. 



2.2.5 The odderon-exchange scattering amplitude in the geometric model 

We can now put everything together and write down the reduced scattering amplitude for odderon 
exchange between two protons described by our geometric model. Since the quarks have a relative 
position fixed by the rigid star shape, we can discard the six position variables x f and use just two 
transverse-space vectors Ri, i = 1,2, which describe the radius and orientation of the star shape for 
each baryon. 

The reduced scattering amplitude J is defined as the scattering amplitude of two protons with fixed 

size and star angle, minus a total factor 2is. It contains the overall factors we have computed so far: 
The colour factor a factor ^ from Equation 2.13 and a conventional factor ^ extracted from the 
type-dependent prefactors. 

J sums up all Feynman graphs with the appropriate prefactors. It depends on the impact parameter 
b and the Ri which parametrise the size and orientation of each proton. 

J{b,Ri,R2) = — — —g^ Ctype ^ ^aibi^a2b2^a3b3 (2-18) 

types {ai}{bi} 
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The second sum runs only over the sets of indices which are compatible with a specific graph type. 
The indices denote the quarks of the first proton the gluons couple to, the 6j the quarks of the 
second proton. Xab is the gluon propagator in transverse space: 

Xab = :^ Ko{m\x-,-xl\). (2.19) 

It is derived in Appendix A. 2. to is a non-zero gluon mass which we have to introduce intermediately 
to compute X. In the gauge invariant expressions which arise after summing up all graphs it can be 
set to zero without causing problems. 

The scattering amplitude is obtained by integrating over the reduced scattering amplitude with a 
wave function (2.1) for each proton. Since we want to obtain the odderon contribution as a function of 
the squared four momentum transfer —t, we have to Fourier transform it back into momentum space. 

T{s,q) = 2is J d^be-''^^ J d^Ri j d^R2 |V(i?i,S)|2 \^iR2,S)\^ J{b,RuR2) (2.20) 

The polar angle in one of the three planar integrations is redundant since the forces between the 
hadrons depend only on their relative orientation. We can choose to measure the angles in the i?i and 
R2 planes relative to the vector b. Then J depends only on b modulus. Hence wc integrate out the 
angle in the b plane and get a factor 2tt and a Bessel function Jq instead of the exponential function 
from the Fourier transform. Using 191 = the result reads: 

T{s,t)=4nis j dbb j d^Ri j d^R2 Jo{V^b) \i;{Ri,S)\^ \'ip{R2,S)\^ J{b,Ri,R2) . (2.21) 

The remaining integrals can only be evaluated numerically. Wc do this using the Vegas integration 
routine from Numerical Recipes, an adapting Monte-Carlo method. In most cases, 6 iterations of the 
algorithm with 5 000000 evaluations of the integrand each were sufficient. For some values of the 
parameters, such as small angles in the star model, 8 iterations a 7000 000 evaluations were used. 

To determine the errors of the integration, it was executed several times with a difi'erent seed value 
for the random number generator on which the Monte-Carlo routine relies. The variation of the result 
proved a better indicator of the numerical errors than the built-in error estimation. The typical error 
in the scattering amplitude is 2 %, though it can be as much as 10 % for very small star angles or 
large t. 

2.3 Calculation of the odderon scattering amplitude from im- 
pact factors in transverse momentum space 

2.3.1 Introduction 

Impact factors are expressions describing the coupling of a Reggeon being exchanged to the scattering 
particles. By that token, the calculations described in the previous section use impact factors in 
transverse position space based on the geometric model. Impact factors in momentum space are 
however far more common. 

It can be shown that the scattering amplitude for a particular process factorises into two impact 
factors and a propagator. Hence, the odderon scattering amplitude is calculated by convoluting two 
odderon-proton impact factors with the propagators of the three gluons forming the odderon. 

= y I (4^ / '^'^^ '"'^^ - - sisiiA,-s,,-52,r ^'-''^ 

Sit,i = 1,2,3 are the transverse momenta of the gluons and Aj = 6u + ^21 + ^3t is the transverse 
momentum of the odderon. 40/3 is the colour factor (equal to d''" d'^'^ ) and 1/3! the combinatorial 
factor reflecting the fact that the gluons are indistinguishable. 



2.3. Calculation from momentum- space impact factors 
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In the high energy hmit A| = —t. This defines only the modulus of A^. Since the scattering 
amplitude depends only on t and not on the vector At, the integral has to be independent of its 
direction. This is true because both the propagators and the impact factors depend only on the 
relative orientation of the vectors 5^. Since 8xt and 62t are integrated over, the direction in which At, 
and therefore S^t, points is irrelevant. Consequently in our calculations we assume that At points in 
one particular direction (rather than averaging over the polar angle in the At plane, for instance). 

The two two-dimensional integrals in Equation 2.22 have to be evaluated numerically. We used the 
Vegas routine from Numerical Recipes with 5 iterations a 5 000 000 evaluations of the integrand. The 
errors in the scattering amplitude were smaller than 0.5 %. 



2.3.2 The general form of the odderon-proton impact factor 

It can be derived from the requirement of gauge invariance that the odderon-proton impact factor has 
to be of the following form: 

3 

% = [i^(At, 0, 0) - ^ F{6it, At - 6^, 0) + 2F{6u, 62t, ^st)] • (2.23) 

i=l 

The first expression in the brackets represents the diagram (a) in Figure 2.4, the second diagram (b), 
and the third corresponds to diagram (c). (Due to gauge invariance, all diagrams have to be there.) 
F(6it, S2t, ^3t) is the proton form factor which can be chosen in different ways. 

The form factor has to be constructed so that the impact factor vanishes when one of the Sa is zero. 
This is required by gauge invariance and has the additional benefit that the integrand in Equation 2.22 
is analytic. 



2.3.3 The form factor of Levin and Ryskin 

A paper on proton-proton scattering by Levin and Ryskin [LR90] uses a Gaussian form factor for the 
proton. It is constructed in such a way that it gives a significant contribution only when all the gluon 
momenta are small. 

F{du, S2U Sst) = exp {-Rl (Sl, + + Sl)) (2.24) 

There is some confusion about the proton radius Rp since it is given in the paper [LR90] as 
= 2.75 GeV^". The more reasonable value of Rp = 0.32723 fm is obtained when assuming the 
wrong sign to be in the exponent of the unit GeV. 



2.3.4 The form factor of Kwiecinski et al. 

In 1997, Kwiecinski et al. published calculations of the odderon exchange in the photo- and electro- 
production of the r]c meson [CKMS97] , based on earlier work by Fukugita and Kwiecinski [FK79] . To 

that end, they compute oddcron-photon and odderon-proton impact factors. For us only the odderon- 
proton impact factor is of interest. It is constructed according to Equation 2.23 from the following 
form factor: 

FiSlt, S2t, Sat) = -2 , 1 A \2 , fj; X \2 , {Ji r \21 " (2.25) 

Its width is proportional to A which was chosen by the authors to be half the rho mass, 384 MeV. 
Unlike the form factor of Levin and Ryskin, it can still give a significant contribution when the gluon 
momenta are large. This is the case when they are large but close to each other, ie when the gluons 
move in the same direction. 
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2.4 The Donnachie-Landshoff fit: A complete set of contribu- 
tions to elastic pp and pp scattering 

In order to fit experimental data of proton-proton scattering, one has to take into account other 
contributions besides the odderon. We decided to obtain those from the phenomenological fit by 
A. Donnachie and P. Landshoff pubUshed in 1984 [DL84]. It is well-known for reproducing the dip in 
the differential cross section in proton-proton scattering. Remarkably, the same parameter values fit the 
data for any centrc-of-mass energy, and the wandering of the dip with changing s results automatically. 

We use the Donnachie-Landshoff fit as a phenomenological input to our calculations. Since it 
depends only weakly on as, in most cases we leave as constant at 0.3, the value used in the original 
work. Whenever it is changed, this will be noted in the results in chapters 3 and 4. 

Donnachie and Landshoff take into account a number of contributions to the elastic proton-proton 
scattering amplitude. They will be described briefly in the following sections. For details see [DL84] . 

2.4.1 All contributions of the DL fit 
The pomeron 

The pomeron is described in the framework of Regge theory (see section 1.4.2). Donnachie and 
Landshoff's pomeron has intercept ap(0) = 1.08 and slope a'p = 0.25 GeV~^. 

To describe the pomeron-proton coupling, Donnachie and Landshoff use a proton form factor of 
the following form: 

^ 4m^-2.79^ 1 
' Arn^-t (l-t/0.71)2 " 

The complete expression for the scattering amplitude is as follows. (It differs by a factor of two from 
[DL84] since our definition of the scattering amplitude is different by that factor, cf. Equation 1.2.) 

Tv{s,t)=2is {3f3Fi{t)f 

Since the pomeron is even under charge conjugation, it contributes equally to proton-proton and 
proton-antiproton scattering. 

The Reggeon contribution 

The Reggeon contribution combines the effects of four different particle exchanges: p and w (which 
are odd under charge conjugation), and / and 02 (which are even). All four particles lie on the Regge 
trajectory with intercept 0.44 and slope 0.93 GeV~^. In addition to the form factor and constant 
factors of the pomeron contribution there is a complex factor. To account for the fact that two of the 
described particles are odd under charge conjugation, it has to be replaced by its complex conjugate 
when describing proton-antiproton scattering. 

Three-particle contributions 

Three different three-particle exchanges are taken into account: triple-gluon exchange, triple-pomeron 
exchange and simultaneous exchange of two gluons and one pomeron. For all of them, the scattering 
amplitude has the following form: 




T3{s,t) =2s / dxidX2dX36{l-Xi-X2-X3) {xiX2X3)'^V{xi,X2,X3y D{si,ti)D{s2,t2)D{s3,t3) . 

-t + to J 



2.4- The Donnachie-Landshoff fit 
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Figure 2.6: The cut-off gluon propagator used by Donnachie and Landshoff. 



Xi are the transverse momentum fractions of the three particles, and Si = xfs and ti = x^t. 
V{xi. X2, = AxiX2Xg, is the vertex function (wave function), and the D{s, t) arc the propagators of 
the particles. The constant A in the vertex function must be fitted so that the triple-gluon contribution 
approaches an asymptotic power law for the differential cross-section: da/dt = 0.09 mb GeV^^ • t~^. 
(See [DL79] for more on the power law.) 

The gluon propagator is ag/t for large \t\. For small \t\ it is cut off with a parabola through the 
origin (see Figure 2.6). It is independent of s. The following ansatz is made for the parabola: 



Depending on what particles are exchanged, three gluon propagators, three pomeron propagators 
or one pomeron and two gluon propagators must be used. The triple-gluon exchange contribution is 
odd under charge conjugation, the other two contributions are even. 

Even though the gluon propagator contains the coupling constant, the Donnachie-Landshoff fit is 
nearly independent of it. This is because the constant A which appears alongside the coupling constant 
is determined by fitting the ggg contribution to the asymptotic power law. Therefore any change in 
the coupling constant in the triple gluon contribution is compensated by a change in A. The other 
three-particle exchanges depend on A and hence change when A varies as a result of a change in Ug- 
However, this has little effect since they are rather less important. 

Reggeon-pomeron and double pomeron exchange 

To obtain an expression for simultaneous exchange of a Rcggcons and a pomeron, their trajectories are 
combined in the following way: aRp(t) = aR(0) -|- q;p(0) — 1 -|- 7'!"''/ t. In the case of double pomeron 

exchange, the analogous formula is apr{t) = 2ap(0) — 5 a'^t. The scattering amplitude is calculated 
from the trajectories of double pomeron and Reggeon-pomeron exchange in the following way: 



Dg{t) = Xt{t + t) for -t<ti. 

A and t are fixed by requiring continuity of Dg{t) and its derivative. This gives A = 2/t\ and t = — | t\. 
The pomeron propagator is derived from its Rcgge trajectory: 
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rr. , X -AiDf3^Fi{t) ,^ / . xaBp(t)-l 

Trp s, t) = / W {A + iB) ( — 

where {A + iB) is the same complex prefactor which occurs in the Reggeon contribution. 



How to produce the dip 

The dip in proton-proton scattering results from a three-way interference between the following con- 
tributions: pomcron exchange, double pomcron exchange and triple gluon (odderon) exchange. The 
constant D in the double pomeron contribution is adjusted so that the imaginary part of pomeron and 
double pomeron exchange cancel at the dip. Since both are nearly imaginary in the vicinity of the 
dip, only a small real part remains. The odderon contribution further diminishes the real part, which 
gives the dip. 

The other contributions are much smaller in size. They serve for fine-tuning the depth and position 
of the dip. 



2.4.2 The odderon contribution in the Donnachie-Landshoff fit 

This section will take a closer look at the triple gluon exchange contribution of the Donnachie-Landshoff 
fit. It evaluates a single Feynman diagram displayed in Figure 2.7. This diagram corresponds to two 
type (c) couplings according to the classification in Section 2.2.4. 



Figure 2.7; The diagram corresponding to Donnachie and Landshoff 's triple gluon contribution. 

This has two implications. First, since only a single diagram is included, the expression is not 
gauge invariant. This is the reason why the gluon propagator has to be cut off to prevent infra-red 
divergences. 

The second implication of choosing this specific diagram is its colour structure: It is the purely 
symmetric | d'^'^ (see Equation 2.17). Taking the antisymmetric structure in spinor space into 
account, it follows that the contribution is odd under charge conjugation. Hence it is an odderon even 
though the authors do not call it that. 

A final remark concerning the triple gluon exchange contribution has to be made and will be of 
great help to anybody reproducing the Donnachie-Landshoff fit. There appears to be a misprint in 
the paper describing the fit [DL84]. The gluon propagator cutoff ti and the scale which corresponds 
to the proton radius are given as "to = ti = 300 MeV^". This value results in a ggg contribution 
which is far too large, dominating the differential cross section for all t. In a later paper about a 
slightly modified fit [DL86] (described below), the authors mention ti as having previously had the 
value 0.3 GeV^. Indeed choosing this value produces a good fit. Therefore, we use this parameter 
value for both to and t\. 

2.4.3 The improved fit of 1986 

While the DL fit is fairly good at ISR energies, it does not fit higher-energy collider data. In Sec- 
tion 4.3.1, this will be shown in detail and attempts to improve it will be presented. 



2.4- The Donnachie-Landshoff fit 
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This state of affairs was already remarked upon by Gauron, Leader and Nicolescu in 1985 [GNL85]. 
In response to their criticism, Donnachie and Landshoff pubhshed an improved version of their fit which 
they claim both extends to collider energies and shows some improvement at ISR energies [DL86]. 
While it seems quite plausible that the Regge-plus-dynamics description of Donnachie and Landshoff 
should be able to accommodate high-energy data, I have not been able to reproduce the improved fit. 

In all probability the authors changed some aspects of the fit without mentioning it in [DL86] . It is 
unfortunate but was to be expected that today they cannot remember any details any more. I followed 
the prescription given in [DL86] but start to disagree with the paper at the point where according to 
the authors the double-pomeron contribution has to be increased significantly. In my reproduction of 
the fit, it in fact has to be made smaller (though very slightly) to make the imaginary part of the total 
amplitude vanish as prescribed in the paper. 



Chapter 3 

The results for pp scattering 



3.1 Experimental data 

Since we use the Donnachie-Landshoff fit as a framework for different odderon contributions, we use 

data for proton-proton scattering in the range where we know the fit is valid. This is the case for 
the energy range investigated by experiments at the Intersecting Storage Rings (ISR) at CERN. Data 
from the ISR were also used by Donnachie and Landshoff themselves. 

The differential elastic proton-proton cross section has been measured by a number of collabora- 
tions at the ISR. Measurements have been made at the centre-of-mass energies 23.5 GeV, 30.7 GeV, 
44.7 GeV, 52.8 GeV and 62.5 GeV. All the data we used were obtained from the Durham/RAL 
database [Dur]. For all energies except 52.8 GeV, we used the data sots as collected by Amaldi and 
Schubert [AS80]. Amaldi and Schubert did not do any measurements themselves for this paper; rather, 
they collected all ISR data and applied a normalisation to data which were not yet normalised. The 
data in the t range relevant for us were measured by Nagy and collaborators [Nag.. 79] and Bohm and 
collaborators [Bo. .74]. The latter data were normalised by Amaldi and Schubert. 

For the centre-of-mass energy 52.8 GeV, we used the data published by Nagy et al. [Nag. .79]. The 
data given by Amaldi and Schubert for that energy do not extend to the t region of interest to us. 
We complement the data from [Nag.. 79] with data measured by Breakstone and collaborators [Bre..85] 
even though they were measured at 53 GeV centre-of-mass energy. It can be seen in the plots that they 
do not always coincide, suggesting some experimental uncertainty. The s dependence of the differential 
cross section is too small to explain the discrepancy. We decided against using the data from Erhan 
and collaborators [Erh..85] since they were available only for relatively few values of t and tended to 
coincide with the Breakstone data. 



3.2 The fit of Donnachie and Landshoff 
3.2.1 Fixed coupling 

Initially we had some problems reproducing the Donnachie-Landshoff fit. This was due to a misprint 
in [DL84] (see Section 2.4.2). The constant to in the expression for the three-particle exchanges (which 
gives the scale corresponding to the proton radius) and the the gluon propagator cutoff ti is 0.3 GeV^ 
rather than the value given in the paper. The original fit reproduced with these parameter values is 
displayed in Figure 3.1. 
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Figure 3.1: The original Donnachie-Landshoff fit with constant coupling. Data for five centre-of-mass 
energies are plotted simultaneously. The scale on the left applies to the lowest curve only. Successive 
curves are shifted upwards by a factor of 100. The points with error bars are experimental data. See 
Section 3.1 for their sources. 



3.2.2 Different cutoffs of the gluon propagator 

To show that the Donnachie-Landshoff fit is independent of the precise nature of the gluon propagator 
cutoff, alternative ways of cutting off were tried. The two possibilities were a flat and a linearly rising 

propagator for small \t\, sec Figure 3.2. 

Changes in the integral over the propagator can be compensated by changes in the constants ti 
and A. The resulting curves are not shown here since they are practically identical to the original fit 
shown in Figure 3.1.^ 



3.2.3 Running coupling 

In a leading-order calculation, the coupling constant is just that, a constant. However, it is known that 
in reality the strong coupling constant depends on the scale. It is possible with limited extra effort to 
take this into account and see how it influences our results. We replaced the strong coupling constant 



-'-At the time this calculation was done it was not yq\. ck'a.r lio-w best to handle the misprint in the paper of Donnachie 
and Landshoff (see section 2.4.2). The scale to which corresponds to the proton radius was chosen as 300 McV'^. Since 
this calculation was quite peripheral to our main interest, investigating the effect of different models for the odderon- 
proton coupling, it was not repeated when we decided on 0.3 GeV'^ as the value most probably meant by the authors. 
Since this change in to did not in the least modify the results of other calculations, it can be safely assumed that this 
result is as valid. 
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Figure 3.2: The two alternative ways of cutting off the gluon propagator. 



with the following expression for the running coupling: 

o.(i) = , |„ 7 ^ (3.1) 

This is the one- loop expression for the running of modified to account for non-pcrturbativc effects. 
The constant 6 was added in the logarithm to avoid the divergence of the perturbative running coupling, 
which is called "freezing" of the coupling constant. This specific way of freezing ag has been used before 
by Ewerz [EwOO] . The number of colours is 3, and we have set the number of flavours Nf also to 3. 
Aqcd was chosen to be 250 MeV. 

The fit with the running coupling is as good as the original fit, see Figure 3.3. 



3.3 The fit with the geometric model odderon 
3.3.1 Fixed coupling 

The main aim of our calculations in the quark-diquark model was to determine the diquark size. For 
that purpose we substituted the diquark model odderon amplitude for the triple-gluon contribution in 
the Donnachie-Landshoff fit. The coupling constant of the DL contributions was left at 0.3 since that 
is what its authors chose. (In fact only the ¥gg and the PPP contributions depend on a^.) We use the 
DL fit as a framework in which several odderon contributions are substituted for the original one. For 
our geometric odderon contribution, we chose the coupling constant 0.4 since it is in the middle of the 
range 0.3 . . . 0.5 common in the literature at the scale given by —t at the dip. See Section 3.3.2 for a 
discussion of other values for the coupling constant. 

The diquark size (respectively the angle in Figure 2.2) was adjusted to give the best fit for a range 
of ccntrc-of-mass energies. The result can be seen in Figure 3.4. The optimal angle is 0.14 tt, or 
approximately 25 degrees (in transverse space) . This corresponds to a diquark size in transverse space 
of 0.22 fm according to Formula 2.2. This diquark size justifies a practice common in non-perturbative 
calculations. Since soft gluons cannot resolve a diquark as small as 0.3 fm, it is legitimate to treat the 
proton as a colour dipole, as in [DFK94, BN99, SS02], for example. 

Choosing the correct angle is vital for obtaining a good fit. Figure 3.5 shows that the fit gets 
considerably worse when the angle is maladjusted. The dependence of the odderon contribution on 
the diquark size will be analysed in greater detail in the next section. 



3.3. The fit with the geometric model odderon 
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Figure 3.3: The Donnachie-Landshoff fit with running coupling. 



3.3.2 The value of the coupUng constant and the significance of the fit 

The quahty of the fit depends strongly on the size of the odderon contribution in the dip region. That 
is why we were able to determine the diquark size (on which the size of the odderon depends) from 
the fit of da/dt. However, the odderon contribution also depends strongly on the coupling constant: 
it is proportional to q^? . Therefore the optimal value for the diquark size depends on the choice of the 
coupling constant (or the choice of Aqcd if cxs is running). 

The difficulty here lies in the fact that in a leading-order calculation the value of the coupling can 
only be estimated. Determining the scale relevant for the coupling and the "true" value of the coupling 
constant requires higher-order calculations. We have therefore determined the optimal diquark size for 
several values of a.. 





star angle 


diquark size [fm] 


0.3 


0.22 TT 


0.34 


0.4 


O.Mtt 


0.22 


0.5 


0.095 TT 


0.15 



Table 3.1: Optimal diquark size for the odderon in the geometric model for different coupling constants. 



The values for which are common in the literature at this scale in the momentum transfer lie 
in the range from 0.3 to 0.5. Table 3.1 shows how the optimal diquark size for several the coupling 
constants. Since the scattering amplitude rises both with the diquark size and with a^, the diquark 
size has to be larger for a smaller coupling constant. All resulting sizes for the diquark cluster in the 
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Figure 3.4: The fit with the odderon computed in the quark-diquark model, with constant coupling. 
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Figure 3.5: The fit with the quark-diquark model odderon for different star angles. The drawn-out 
curve is the best fit, 0.14 tt. The centre-of-mass energy is 44.7 GeV. 



proton arc < 0.35 fm, supporting the assumption of a reasonably small diquark unless the coupling 
constant is smaller than 0.3. 



3.3. The fit with the geometric model odderon 
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The other extreme, a rotationally symmetric "Mercedes star" proton, can be excluded with some 
confidence. This would imply an unrealistically small coupling constant of 0.17. 
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Figure 3.6: The dependence of the quark-diquark model odderon amplitude on the diquark size. The 
error bars reflect the numerical errors of up to 10%. The three points for each value of the diquark size 
correspond to three different transverse momentum transfers. The square of the momentum transfer 
-t is 0.6 GeV^, 1 GeV^ and 1.6 GeV^ from top to bottom. 



3.3.3 Properties of the geometric model odderon 
Dependence on the diquark size 

As can be seen in Figure 3.5, the quality of the fit depends strongly on the diquark size. The drawn-out 
curve in the figure is the overall best fit for all centre-of-mass energies. Though a better diquark radius 
could be found for one particular energy, it would give a worse fit for other energies. 

The odderon amplitude's dependence on the diquark size is plotted in Figure 3.6. Data for three 
different values for —t are in the plot. As one can see the three sets of points differ by a constant factor. 
This means that the suppression with vanishing diquark size is independent of the transverse four- 
momentum transfer. (The odderon is by construction completely independent of the centre-of-mass 
energy.) 

The odderon-proton coupling vanishes for zero diquark size [Za89, RD96], for the following rea- 
son: A proton with a point-like diquark is a colour dipolc. Since the odderon is odd under charge 
conjugation, it gives a contribution of opposite sign when the dipole is rotated by 180 degrees. After 
averaging over all orientations of the dipole the coupling has to be zero because contributions from 
opposite orientations cancel out. 

Figure 3.7 shows a comparison between the perturbative odderon in the quark-diquark model and 
the non-perturbative odderon computed by M. Riiter with the Model of the Stochastic Vacuum [Rue97]. 
The non-pcrturbativc odderon vanishes significantly faster than the perturbative one. This can be 
understood with physical arguments: The non-perturbative odderon which consists of soft gluons 
cannot resolve small diquarks and takes them to be point-like. 

To determine the power law according to which the odderon amplitudes vanish, a straight line was 
fitted to the computed amplitudes in a double logarithmic plot. The result can be seen in Figure 3.8. 
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Figure 3.7: Comparison between the non-perturbative and the perturbative odderon. The suppression 
function plotted here is defined as the scattering ampHtudc for a given diquark size divided by the 
scattering amphtude for the symmetric three-star with angle 2/3 tt. The non-perturbative odderon 
amplitude vanishes significantly faster than the perturbative one. The non-perturbative data were 
obtained from [Rue97]. 
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Figure 3.8: Double logarithmic plot of the data from Figure 3.7. The power laws which describe 
the vanishing of the amplitudes for small diquark radii were obtained by fitting a straight line. The 
two lines intersect at a diquark size of do = 0.868 fm which means a symmetric star and where the 
suppression function is 1. 



3.3. The fit with the geometric model odderon 
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The perturbative odderon vanishes as ^dii^ark' ^-^^ non-perturbative one as ^dii^uark* 
Relative importance of different graph types 

Since we classify the Feynman graphs in different types, it may seem easy to investigate the contri- 
butions of graphs belonging to one particular type. Unfortunately, the numerical integration over the 
reduced scattering amplitude is much more difficult when only one graph type is computed. Large 
cancellations between the different graph types are to be expected since contributions from single graph 
types are not gauge invariant. However, the reduced scattering amplitudes of single graph types seem 
to vary wildly, making numerical integration difficult. 

Nevertheless, calculations have been done to determine the relative size of the contributions of 
different graph types. 10 times 7 000 000 evaluations of the integrand were performed to evaluate 
the integral with the Vegas algorithm. The results carry large errors and are to be understood as 
indications of the order of magnitude rather than numbers of any precision. Table 3.2 shows the 
results. 



Graph type 


Scattering amplitude [fm/GeV] 


(aa) 


-70 


(ab) 


375 


(bbl) 


-190 


(bb2) 


-430 


(ac) 


-205 


(be) 


870 


(cc) 


-215 


total 


0.0127 



Table 3.2: The very approximate scattering amplitudes corresponding to each graph type. The trans- 
verse momentum transfer is 1 GeV^; the star angle in transverse space is 0.14 tt, the value found to be 
most realistic. The numerical errors are 30 to 50 %, except for the total. 

The graph type (be) gives the largest contribution and (aa) the smallest. The other graph types 

are clustered around 200 to 400 fm/GeV. The contributions add up to a number which is compatible 
with 0, considering the large errors. This is reassuring since the total scattering amplitude is much 
smaller than the individual contributions. 

Donnachie and Landshoff [DL84] argue that the type (cc) contribution is dominant at large t. They 
note that if the three gluons do not couple to three different quarks, the quarks' trajectories are not 
parallel immediately after the scattering. Therefore a further gluon has to be exchanged before the 
two triples of quarks can again form protons. This and the fact that one of the quarks has to be 
off-shell intermediately suppresses the couplings of types other than (c). Consequently, Donnachie and 
Landshoff 's triple gluon contribution contains only a graph corresponding to type (cc). 

In our residts, the (cc) contribution is not dominant. This indicates that the argument of Donnachie 
and Landshoff is not yet valid &it — l GeV^. That said, it has to be emphasised again that our results 
concerning the scattering amplitude by graph type are very imprecise. 

3.3.4 Running coupling 

The diquark model odderon can equally be inserted into the Donnachie-Landshoff fit with running 
coupling. The coupling constant in both the geometric model odderon and the Donnachie-Landshoff 
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Figure 3.9: The fit with the odderon computed in the quark-diquark model, with running couphng. 

contributions is substituted by the expression for the running coupHng (3.1). 

The result is nearly as good as the previous fits, though the dip tends to be too deep for the higher 
centre-of-mass energies. The plot can be seen in Figure 3.9. The best fit is achieved with the same 
star angle as for fixed coupling, 0.14 tt, corresponding to a diquark size of 0.22 fm. 

3.4 The fit with momentum-space impact factors 

3.4.1 The impact factor of Levin and Ryskin 

If one uses the impact factors with the Gaussian form factors as proposed by Levin and Ryskin 
in [LR90], the coupling constant has to be adapted slightly. To obtain the best fit, has to be chosen 
as 0.5. The curves for all centre-of-mass energies are displayed in Figure 3.10. 

3.4.2 The impact factor of Kwiecinski et al. 

In [CKMS97], Kwieciiiski et al. use the coupling constant = 1. Using this value for proton-proton 
scattering makes the odderon contribution far too large. It dominates over all other contributions. In 
Figure 3.11 on page 42 it is shown together with data for \/i = 44.7 GcV. 

A good fit can be obtained by changing the coupling constant.^ The best fit was obtained for 
as = 0.3 and is presented in Figure 3.12. That such a small coupling constant is required throws some 

^It was decided not to modify the constant A which determines the width of the form factor (2.25). It is half the 
mass of the p meson, the hghtest vector meson. Since the structure of the proton is measured in -yp scattering among 
other processes and the photon has to fluctuate into a vector meson to interact strongly with the proton, this mass is 
closely connected to the structure of the proton. 



3.5. Proton-proton summary 
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doubt on results for diffractive rj^ production through odderon exchange at HERA [CKMS97, EIKS98, 
BBCVOl]. The enormous dependence of the odderon on the couphng constant (oc a^) means that a 
factor of 0.3 in changes the cross section by orders of magnitude. Furthermore, while the situation 
in the process of r]c electroproduction could be said to be different from that in elastic pp scattering, 
the coupling constant would be expected to be smaller still due to the hard scale given by the rjc mass. 
Our results indicate that the amplitude for 77c production was overestimated by the three groups cited 
above. 



3.5 Proton-proton summary 

Figure 3.13 shows a summary of the four methods used to compute the odderon contribution: Don- 
nachie and Landshoff 's ggg contribution, our geometric model, and the impact factors with form factors 
by Kwiecihski et al. and Levin and Ryskin, respectively. The coupling is fixed in all cases. As can be 
seen, all methods yield good results after the coupling constant is adapted. 

Figure 3.13 shows clearly that the quality of the data is not sufficient to distinguish between 
different models for the odderon-proton coupling. However, the data can strongly restrict the range of 
parameter values for each individual model. The diquark size in the geometric model in position space 
could be fixed to 0.22 fm for a coupling constant of 0.4 and was found to be < 0.35 fm for all reasonable 
values of as . It could be shown that the coupling constant used by several groups [CKMS97, EIKS98, 
BBCVOl] would be far too large for the situation of proton-proton scattering and that they probably 
overestimated the cross section for rjc production through odderon exchange. For the impact factors 
computed from the Gaussian form factor of Levin and Ryskin, the coupling constant also had to be 
adapted, but less radically. 
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Figure 3.11: The fit with the odderon-proton impact factor by Kwiecinski et al. with coupUng constant 
ag = 1 and 0.3, respectively. 
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Figure 3.13: Summary of proton-proton scattering. All methods for computing the odderon contribu- 
tion are compared for all centre-of-mass energies for which experimental data are available. 



Chapter 4 

The results for pp scattering 



4.1 Experimental data 

For centre-of-mass energies in the range of our proton-proton data, the differential proton-antiproton 
cross section has been measured at only one energy, 53 GcV. The measurement was done at the CERN 
Intersecting Storage Rings (ISR) by Breakstone and collaborators [Bre..85]. 

At higher energies measurements have been made at the SppS collider, also at CERN. The proton- 
antiproton differential elastic cross section has been measured at the energy 546 GeV by Bozzo and 
collaborators [Boz..85] and at 630 GeV by Bernard and collaborators [Ber..86]. 

The measurement at Fermilab by Amos and collaborators at 1.8 TcV [Am. .90] is of no use to us, 
unfortunately. It does not extend to the values of the squared four momentum transfer —t at which 
the data show a kink, where the odderon begins to play a role. 



4.2 pp scattering at the ISR 

To treat proton-antiproton scattering, the sign of the odderon amplitude has to be reversed relative 
to proton-proton scattering. The contributions of the Donnachie-Landshoff fit which describe C-odd 
particles have likewise to be modified as described in Section 2.4.1 and in [DL84]. 

Figure 4.1 compares four calculated curves to the experimental data at the energy 53 GeV. The 
dotted line is the original Donnachie-Landshoff fit. The solid line shows the fit with the odderon 
computed in our geometric model for the odderon-proton coupling. The diquark size was the same as 
for proton-proton scattering, 0.22 fm. 

The two dashed lines show the fit with the odderon contribution computed from impact factors in 
momentum space, with form factors according to Kwiecihski et al. and Levin and Ryskin, respectively. 
The coupling constant was the same as the one used for proton-proton scattering: 0.3 for Kwiecinski's 
form factor and 0.5 for Levin's and Ryskin's. 

That the same parameter values as for proton-proton scattering describe proton-antiproton data 
as well adds some weight to our results in the previous chapter. It is also a tribute to the Donnachie- 
Landshoff fit. 

As was the case for proton-proton scattering, all models for the odderon-proton coupling fit the 
data equally well. 
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Figure 4.1: The computed proton-antiproton differential cross section compared to experimental data 
from [Bre..85]. The centre-of-mass energy is 53 GeV. 

4.3 Describing higher-energy data with an energy dependent 
odderon 

4.3.1 With the Donnachie-LandshofF fit as framework 

Since all the models for the odderon-proton coupling describe proton-antiproton scattering at ISR 
energies so well, it is interesting to see whether our calculations can be extrapolated to higher energies. 
Since all models led to equivalent results, we restrict this investigation to the Donnachie-Landshoff fit 
with the original triplc-gluon contribution as an odderon. 

As Figure 4.2 shows, a pure extrapolation is wide off the mark. This may be explained by the fact 
that Donnachie and LandshofF's ggg contribution (unlike other contributions) has no dependence on 
the centre-of-mass energy ^/s. In the relatively narrow range of energies measured at the ISR this is 
not necessary, but when extrapolating to much higher energies, an energy dependence for the odderon 
might improve the fit. 

We attempted to improve the fit for the proton-antiproton data at higher energies by adding an 
energy dependence to the odderon. This could in turn enable us to learn something about its s 
dependence, specifically its Rcgge intercept. 

We gave the odderon contribution an s dependence by multiplying its scattering amplitude with 
the s dependent factor known from Regge theory (see Section 1.4.2, Equation 1.6). As a first step, 
only the Regge intercept was used but not the slope: 



This amounts to taking into account the dependence of the odderon propagator on ,s but not on t. 
All the t dependence comes from the odderon-proton coupling contained in the original triple-gluon 
expression of Donnachie and Landshoff. Their expression corresponds to an intercept ag = 1. Since 
the intercept is not expected to deviate much from 1, neglecting the s dependence for the range of ISR 
energies makes sense. 




(4.1) 
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Figure 4.2: A pure extrapolation of the DL fit to higher energy pp scattering matches the data poorly. 
Data from [Ber..87]. 

Figures 4.3 and 4.4 show the fit with the energy dependent odderon for different intercepts. As a 
reference energy we chose ^/so = 53 GeV since the fit works well at that energy. Figure 4.3 suggests a 
Regge intercept of 1.4 for the odderon. This is the value that produces the best fit. However, Figure 4.4 
at -^s = 630 GeV which extends to larger |t| shows clearly that none of the curves matches. After 
the customary kink (which looks a bit like a dip at this energy) the data points start to fall off more 
rapidly again. The computed curves which are of the right magnitude do not reproduce this behaviour. 
Furthermore, such a high intercept already has an effect in the ISR range. In proton-proton scattering 
it makes the dip deeper for the higher ISR energies and flatter for the lower ones. This is the opposite 
of the behaviour shown by the data (see for instance Figure 3.1). 

This disappointing result is not in contradiction with the little which is known about the odderon 
intercept. The odderon intercept calculated perturbatively in the Leading Logarithmic Approximation 
is slightly below one.^ Given that the intercept of the BFKL pomeron was made smaller rather than 
larger by next-to- leading-logarithmic corrections, an intercept of 1.4 for the odderon seems unlikely. It 
is therefore no great surprise that calculations using this intercept fail to match the data, even if the 
fit was satisfactory at one energy. 

One can go one step further and use the complete Regge trajectory of the odderon, taking also the 

t dependence of the odderon propagator into account. That has the effect that the s dependence of the 
odderon contribution is different for each value of t (or vice versa). The expression for the scattering 
amplitude is the following: 



Figure 4.5 shows the result for some combinations of values for the intercept and the slope. As 
was to be expected, the already very high intercept had to be raised even more. The positive slope 
diminishes the odderon contribution relative to t = since t is negative. To retain the same size of the 
odderon contribution at —t = l GeV, for instance, the intercept has to be raised by 0.1 if the slope is 
set 0.1 GeV-2 higher. 

^For the new solution of the BKP equation found be Bartels, Lipatov and Vacca [BLVOO], it is exactly one. 




(4.2) 
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Figure 4.3: The fit with an s dependent odderon for proton-antiproton scattering at y/s = 546 GeV. 
The odderon intercept ranges from 0.8 to 1.5 in steps of 0.1. The reference energy was ,/so = 53 GeV. 
Only the Regge intercept is taken into account, not the slope. 




Figure 4.4: The fit with an s dependent odderon for proton-antiproton scattering at -^s = 630 GeV. 
The odderon intercept ranges from 0.8 to 1.5 in steps of 0.1. The reference energy was ^/so = 53 GeV. 
Only the Regge intercept is taken into account, not the slope. 
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Figure 4.5: Taking the whole Regge trajectory of the odderon into account for proton-antiproton 
scattering at \/s = 630 GeV. The plot displays the combinations of intercept and slope which came 
closest to matching the data. The reference energy was ^/so = 53 GeV. 
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Figure 4.6: Fit with the energy dependent odderon, taking into account its Roggc trajectory and the 
signature factor. The centre-of-mass energy is \/i = 630 GeV. The Regge intercept is 1.7, the slope 
0.2 GeV~^. The reference energy was again yfso = 53 GeV. 
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Figure 4.5 suggests a Regge slope of 0.2 GeV~^ and an intercept between 1.6 and 1.7. Such a large 
intercept seems extremely unlikely for the reasons stated above. Besides, the fit is far from perfect: 
The flattening off in the data points between ~t = 0.8 and 1.2 GeV~^ is not reproduced well. 

One final attempt was made to improve our description of higher-energy pp scattering. Wc have so 
far neglected the odderon's signature factor which changes its phase according to its Rcgge trajectory 
(see Section 1.4.2). This amounts to assuming that the trajectory stays close to one, ic the intercept 
is close to one and the slope close to zero. Since wc found quite different values for both this is no 
longer consistent. The expression of the odderon amplitude with the signature factor reads: 



The result is disappointing, as Figure 4.6 shows. Far from improving the fit, the added signature 
factor makes it worse. The curve falls off even more rapidly in the region between 0.8 and 1.2 GeV^. 

It can be concluded that for extending the Donnachie-Landshoff fit to higher energies, modifying 
the odderon contribution alone is not enough. It seems that some of the other contributions are no 
longer correct at energies above the ISR range and also have to be modified. 

4.3.2 With the fit by Gauron, Leader and Nicolescu 

Another attempt was made to determine the odderon intercept and slope. A different fit for da/dt was 
used as a framework. It was created by Gauron, Leader and Nicolescu in 1989 [GLN90]. It provides 
a much better fit than the one due to Donnachie and Landshoff [DL84] for both proton-antiproton 
and proton-antiproton at the cost of a large number of free parameters.^ In particular, it also fits 
higher-energy proton-antiproton data. 

Gauron, Loader and Nicolescu take into account a rather large number of contributions to the 
scattering amplitude. First, they include two "asymptotic" contributions which unlike the others 
remain significant in the limit s — > oo. One of them is a C = +1 "froissaron" , the other a C = — 1 
maximal odderon.'^ Besides there are contributions which play a role at lower (ISR) energies: a pomeron 
pole, a double pomeron cut, an odderon pole, an odderon-pomeron cut, three types of Reggeons and 
two Reggeon-pomeron cuts. 

To siibstitute our odderon contribution, wc leave out from the fit the three contributions related 
to the odderon — the maximal odderon, the odderon pole and the odderon-pomeron cut — and add 
our geometric odderon. However, far from being able to study the s dependence, we failed even to 
produce a satisfactory fit for a single energy. Figures 4.7 and 4.8 show curves with a range of different 
coupling constants for our odderon. 

A detailed analysis reveals that the maximal odderon is by far the most important of the three 
odderon contributions in the GLN fit. Due to interference between the terms it contains, it already 
shows a dip structure itself. Figure 4.9 shows a comparison between the modulus of the maximal 
odderon amplitude and our odderon (which is real- valued). Apparently the odderon plays an even 
greater role in the GLN fit than it does in the DL fit. With an odderon with a much simpler structure 
such as ours, we cannot reproduce their fit. 

4.3.3 Consequences for the odderon 

As Figure 4.9 shows, the maximal odderon used by Gauron et al. has a much richer structure than 
our geometric odderon and, by implication, the odderon in the Donnachie-Landshoff fit. It seems 
implausible that the GLN odderon should work in the framework of the DL fit and vice versa. 

^Donnachie and Landshoff have improved their fit in response to criticism by Nicolescu et al. [DL86]. However, I 
could not reproduce the improved fit (sec Section 2.4.3). 

^It is maximal in the sense that it is oc (logs)'^, the maximum growth allowed by asymptotic theorems. 
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Figure 4.8: The GLN fit with our geometric odderon substituted for the three odderon-related contri- 
butions. Proton-antiproton scattering at ^/s = 630 GeV. 
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Figure 4.9: Modulus of the GLN fit's maximal odderon contribution compared to our geometric odd- 
eron for = 0.4. The maximal odderon has a much more complex structure. The centre-of-mass 
energy is ^/s = 44.7 GeV. 

The fits indeed do not work with the wrong odderon contributions: Figure 4.10 shows the fits 
by Donnachie and Landshoff and Gauron, Leader and Nicolescu, respectively, with the other group's 
odderon substituted, together with the unchanged fits. The fits with the wrong odderon match the 
data quite badly. 

This has disturbing implications for the odderon itself. While many people agree that there is 

an odderon, there is no consensus on what it should look like. This comparison shows that we do 
not know the universal odderon given by the location of the odderon pole in the complex angular 
momentum plane. Though the odderon plays an important role in both fits, the form of the odderon 
contribution is so different as to be incompatible. While in principle C-evcn and C-odd contributions 
can be determined independently of each other, in practice the lack of more data and the freedom in 
the C-even terms mean that it cannot. The only explanation for the two very different odderons in 
the two fits is that they make a different separation between the even and odd contributions to the 
scattering amplitude. (Their non-odderon odd contributions are similar.) Since both fit the data well, 
the data are not yet able to decide between the two conflicting views. 
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Figure 4.10: Comparison of the unchanged fits by Donnachie and Landshoff (DL) and Gauron, Leader 
and Nicolcscu (GLN) with fits in which the other group's odderon has been substituted. This is 
proton-proton scattering at -^i = 44.7 GeV. 



Four-pomeron vertices from 
unitarity corrections 



This part concerns vertices occurring in unitarity corrections to the BFKL equa- 
tion. A 1 — > 3 pomeron vertex is computed from a 2 — > 6 reggeised gluon vertex 
occurring in the six-gluon ampHtude in the Extended Generalised Leading Logarith- 
mic Approximation (EGLLA), part of the gluon-ladder approach to higher-order 
BFKL. This gluon vertex is part of the integral equation for the six-reggeised-gluon 
amplitude in EGLLA. It is an irreducible element in the conjectured conformal field 
theory of unitarity corrections. 

As a preparation, five basic functions of which the I n pomeron vertices are 
ultimately composed are presented in both momentum and configuration space. 
A further function G, which combines them and in terms of which the pomeron 
vertices can be expressed, is shown to be conformally covariant. 

The irreducible 1 — > 3 pomeron vertex is computed in terms of integrals over 
conformal eigcnfunctions. It is cast into the form of a conformal four-point func- 
tion. The freedom left by this form, a function of an anharmonic ratio of the four 
pomeron coordinates, is computed for all terms of the vertex and expressed in 
terms of three integrals. The partly reggcising part of the six-gluon amplitude also 
gives rise to 1 ^ 3 pomeron vertices, which dominate in the limit Nc oo. They 
are computed in the last sections of this work. 

These vertices arc of great interest because 1 ^ n pomeron vertices for n > 2 
were shown not to exist in the dipole approach to higher-order BFKL. The findings 
of this work, which were calculated in the gluon-ladder approach, demonstrate a 
discrepancy between the two major approaches to higher-order BFKL. 



Chapter 5 

Fundamental functions 



5.1 Introduction 

This chapter collects the preparatory work which has to be done before we start on our derivation of 
pomeron vertices in the following chapter. Section 5.2 contains an introduction to conformal trans- 
formations, relevant elements of conformal field theory and the BFKL equation. Sections 5.3 and 5.4 
present five functions of which reggeised gluon vertices are composed, and derive their form in config- 
uration space. Section 5.5 deals with the function G, a combination of these five functions in terms of 
which reggeised gluon vertices are most easily expressed, and determines its properties under confor- 
mal transformations. None of this is truly original work. Notably the derivation of the transformation 
properties of G to some extent repeats the proof of the conformal invariance of the 2^4 reggeised 
gluon vertex by Bartels, Lipatov and Wusthoff [BLW96] , which is a sum of nine terms of G functions, 
and the sometimes sketchy treatment of G by Braun and Vacca [BV99] . Likewise, the five functions of 
which G is composed have bccin known in momentum space for some time. All the same, this is to my 
knowledge the first systematic treatment of the five functions and G in both momentum and position 
space and will hopefully provide a clear and reproducible guide for future students of the subject. 



5.2 Prerequisites 

5.2.1 Conformal transformations 

The topic investigated in this and the following chapter is the behaviour of expressions under conformal 
transformations. Conformal transformations (or Mobius transformations) are those transformations 
in the complex plane which preserve both angles and their orientation. Rotations and dilatations 
(rescalings) are conformal transformations. So arc translations. 

The complex conjugation preserves angles but inverts their orientation and therefore is no conformal 
transformation. The inversion, z ^ \/z, however, conserves orientations as well as angles. It is not 
quite the geometric inversion with respect to the unit circle. The inversion with respect to any circle 
changes the orientation of angles, just like the reflection with respect to any straight line. The geometric 
inversion with respect to the unit circle can be written bs z ^ 1/z m. complex coordinates. Hence 
the complex inversion is (graphically speaking) the concatenation of the geometric inversion and the 
complex conjugation (reflection with respect to the real axis), which is consistent with the fact that it 
conserves orientations. 

There is a general form for conformal transformations: 

^^az + b ad-bc^O. (5.1) 

cz + d ^ ' 
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However, it would not be easy to determine how an expression transforms under such a relatively 
complicated transformation law with four parameters. 

In fact, all Mobius transformations can be generated from the four we have already mentioned: 
rotations, dilatations, translations and the inversion. Therefore it is sufficient to investigate their effect 
separately. Their formulaic expressions are as follows: 

z^e'"^ z ip€R 
z^cz ceM>o 
z 1-^ z + a a G C 
1 

2 1-^- 
z 

There is one point which should be mentioned to avoid confusion with other work concerning con- 
formal transformations. What we regard as conformal transformations here differs from the definition 
of conformal transformations in the theory of complex functions and in string theory. There, conformal 
transformations are defined as functions which arc locally conformal, ie whose Jacobi map preserves 
angles and their orientation. This holds for any analytic function with a non-vanishing derivative. 
Here, however, we are interested only in global conformal transformations. 



5.2.2 Conventions 

Before starting out in earnest, let me explain a few fundamental points of my notation. Variables with 
two indices are customarily defined to be the difference between two coordinate vectors, ie the distance 
vectors between two points. 

Pi2 := Pi - P2 

The derivative with respect to a coordinate is abbreviated by a differentiation symbol with the 
same index, ie 

dpi 

The same convention will be used for Laplacians and other derivative operators. 

I adopt the following convention for the constant factors in Fourier transformations (for an arbitrary 
function / depending on a two-dimensional vector): 

f{p) = /(fe) e'''" ; /(fe) = / rfV f{p) e-"'" . (5.3) 

The function and its Fourier transform will not be denoted by a different symbol; the distinction will 
be clear from the type of its argument. 



5.2.3 The vectorial and the complex formulation of expressions in two- 
dimensional space 

We will describe aspects of high-energy scattering in two-dimensional transverse space. There are 
two ways to describe a location in the plane: as a two- vector or a complex number. While the 
formulation in terms of vectors is mostly the simpler and more intuitive, the complex notation is much 
more convenient when applying conformal transformations. In the following the two notations will be 
distinguished by writing vectors as bold face characters and complex numbers as plain ones. 

Vectors can be transformed into complex numbers and vice versa. The two components of two- 
vectors arc; identified with the real and imaginary part of the corresponding complex number. In the 
context of conformal field theories, however, the real and imaginary part are not customarily seen as the 
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two degrees of freedom of a complex number. Instead, the complex number and its complex conjugate 
are viewed as independent degrees of freedom. While this approach is hard to grasp intuitively, it is 
practical when dealing with conformally invariant expressions. After performing the calculations, one 
imposes again the physical condition that p* be the complex conjugate of p. The conversion formulas 
to and from the complex notation are: 

p = ei - p + ie2 - p, p = ei - p-ie2 - p, where ei = , 62 = 



If p + p* 



P 2\-i{p-p*) 

The scalar product of vectors can easily be transcribed into the complex notation by using (5.4). 
To distinguish the scalar product from the ordinary product of complex numbers, we write it with 
parentheses. 

Pi - P2 = \ {pip*2 + pIP2) =■■ (pi, P2) (5.5) 

In addition to points in the plane, we need differential operators, notably V and A. The derivation 
of the first-order differential operators will not be performed explicitly here. It is straightforward; 
one can obtain it for instance by requiring the first-order term in a Taylor expansion of an arbitrary 
function to be equal in both notations. The Laplacian is obtained as the scalar product of two Vs. 
The results are: 



v = 



d + d* \ (5.6) 
i{d-d*)) 



A = 4dd* 

Note that the conversion formula is different from the one for the coordinates. In rare cases, I will use 
nabla operators in the complex notation to indicate gradients. They are converted from the vector 
nablas with the formula for coordinates (5.4), and consequently V = 29* and V* = 2 9. 

Finally, there are integrals. The Jacobi determinant of the transformation (5.4) is —2i, hence the 
integral measure in both coordinate systems differs by a factor of 2. Knowing that, we can relate 
vectorial integrals and delta functions to integrals over and delta functions of the complex coordinates: 

I dPpfip) = I <fpfip,p*) = ^ J dp J dp* f{p,p*) 

5\p)^5\p)^2 5{p)5{p*) 

Two-dimensional integrals or delta functions of complex coordinates are customarily defined to include 
the constant factors and are therefore equivalent to the vectorial ones, as shown in the equations. 

5.2.4 Conformal eigenfunctions 

This and the following section will introduce the few elements of conformal field theories we will need, 
conformal eigenfunctions and n-point functions. For a thorough introduction to conformal field theory, 
see [Dot88, Gi91, Sch96]. 

Eigenstatcs of operators which are invariant under conformal transformations arc severely restricted 
in their form. They are called conformal eigenfunctions. Conformal eigenfunctions of integral operators 
with two arguments (like the BFKL kernel, see below) have the following form (except for constant 
factors independent of the ps): 

l+n , ■ 1-n , ■ 

PlaP2aJ \P\aP2aJ 
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They are parametrised with a real number u and an integer n. v is called conformal dimension, 
n conformal spin. The exponents in Equation 5.8 are called conformal weights and denoted with the 
symbols h and h. So we have: 

1 + • 7 1 - ' i 

n = — Viu n = — ^ 11^ = 1 — n 

_ 1 _ (5-9) 

n = h — h i' = — {h + h—l) 

It is well known from the theory of complex functions that an arbitrary powc;r of a complex base to 
a complex exponent is not well defined. The value of the power .x^ := exp(y log./;) depends on which 
branch of the multi- valued logarithm one chooses. So there is a certain ambiguity in the definition (5.8) 
which is never remarked upon. The most natural choice is to choose the principal value logarithm 
denoted by Log, the imaginary part of which lies in the range (— tt, tt]. 

Having defined the expression (5.8), one can work out how the transform under exchange of 

the first two arguments. This transformation amounts to changing the sign of the base of the power, 
ie the argument of the logarithm. The exact result is: 

r (-1)" for r G C\M, 

(P2, Pi; Pa) = £;("'"Hpi,P2;Pa) • < -exp(-27ri/) for r e R>o , (5.10) 

[ — exp(27r!^) for r € ]R<o , 

where r = p\2l{p\apia)- For real r, Logr and Logr* transform in the same way {j^i is added to or 
subtracted from both) . This has the consequence that not the ^ part but the rest of the exponents is 
retained and determines the transformation property. 

It can be safely assumed that the transformation law intended by the authors of the conformal 
eigenfunctions was the first line of (5.10) only.^ Reluctantly in line with physicists' usual nonchalance 
towards mathematical rigour, I will assume in the following that the eigenfunctions transform with 
a factor (—1)", ie are symmetric for even n and antisymmetric for odd n. A limited justification for 
this is that the deviation from that rule is restricted to a set of measure zero, which can never be 
measured experimentally and will furthermore be irrelevant under an integral. However, this remains 
an abuse of mathematics. The desired transformation law could only be obtained be choosing two 
diflPerent branches of the logarithm for the non- conjugated power and the conjugated one, so that 
log*'"'^-'(— r) — log'-"'^-' r = — (log'-^-'(— r*) — log*-^-* r*) for all r. This is obviously mathematical nonsense. 

In the later chapters, we will restrict ourselves to the ground state with zero conformal spin. Here 
is its explicit form, which depends only on the moduli of vectors, not their direction: 




(5.11) 



The third coordinate argument of the eigenfunctions represents an external coordinate of the cor- 
responding fields. We will often suppress it in the calculations, but it is always implied. 



5.2.5 Conformal n-point functions 

In a conformal field theory there are special fields called "primary fields" . They are distinguished 
by satisfying the following relation under conformal transformations p'-? 

^Kh{p.P*) = {^)'^n,H{p',p'l- (5.12) 

^The original work [Li86, Li89] does not remarlj on the symmetry of the nor on the ambiguity of complex 

exponentials. 

■^Most authors formulate (5.12) as a transformation law rather than an equation. But since for physical reasons 
$'(p',p'*) = $(p, p*) is required, the above relation has to hold as an equation for primary fields. 



58 



Chapter 5. Fundamental functions 



It is easily checked that the conformal cigenfunctions (5.8) fulfil this requirement. In this case, h and 
h are related via (5.9), even though in general they can be independent variables. 

As a consequence of this transformation property, n-point functions of primary fields are constrained 
in their form. The two- and three-point functions are fixed up to a constant: 

Pl2 Pl2 

„ _ _ _ [„hi+h2-h3 h2+h3-hi hi+h3-h2 hi+h2-h3 * h2+h3-hi * hi+h3-h2'\ ~^ 

— (^hiMMMMM lPi2 P23 P13 P12 P23 P13 J 

The conformal cigenfunctions introduced in the previous section can themselves be imdcrstood as 
conformal three-point functions. In the context of the BFKL equation, two of the participating fields 
can be interpreted as rcggciscd gluons, one as a pomeron. The conformal cigenfunctions, interpreted 
as three-point functions, describe the pomeron in terms of reggeised gluons. The explicit formula is 
the following, where the correlation between n, h and h is as in (5.9): 

£('^'")(pi,P2;Pa) = {^0,0{PUP*l)^0,0{p2,P*2)^h:h{Pa,P:))- (5.15) 

The four-point function has some more freedom. It is fixed only up to a function which may depend 
on an anharmonic ratio of the four coordinates, though not on the individual coordinates. 

{^h,M (^1' Pi) ^h2M iP2,P*2) ^h3,h3 (/'3, Ps) %„-hAp4, PD) = 

i<j 

P12P34 
X = 

P13P24 

At this point some words about anharmonic ratios are helpful. In principle, one can construct six 
different anharmonic ratios from four coordinates: Both numerator and denominator can contain three 
different pairs of differences of coordinates. If the fraction is not to be constant = 1, numerator and 
denominator must not be equal and six possibilities remain. These are pairwise inverses of each other. 
They all can be expressed in terms of one anharmonic ratio, eg x. Here are three of them (the other 
three being their inverses): 



P12P34 1 P14P23 X P12P34 

^ ^3 — ^ 

P13P24 ' P13P24 ' 1- X P14P23 

For example, here is the proof for the last identity: 



(5.17) 



P12P34 

P14P23 



/ Pl4p23\ ^ ^ / P1P2 - PiP2 - Plp3 + P4P3 \ ^ ^ 
VP12P34/ V P12P34 / 

/ 1^ _l_ -P1P2 + P3P2 + PiPA - PsPi + P1P2 - PiP2 - Plp3 + P4P3 \ 

\X P12P34 / 

/ 1- _l_ P3P2 + PlP4 - P4P2 - Plp3 \ ~^ = i _ 

\x P12P34 / 



5.2.6 The BFKL equation 



The BFKL equation describes the evolution in rapidity of the system of two reggeised gluons which 
makes up pomeron exchange in QCD. It is an integral equation of Bethe-Salpeter type. The unknown 
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in the equation is a four-point function of reggeised gluons, or equivalently the pomeron propagator. 
The BFKL equation in momentum space reads as fohows: 



^2 

2^"rf2fe'2' 1 1 



(5.18) 



(27r)3 fef k 



772 7^772 Kbfkl (^i j ^2 ! ^i) ^2) 4'<^ ('^i) ^2; fei , ^2) • 



It takes the form of an eigenvalue equation for the amphtude (f)^ with the eigenvalue w, the complex 
angular momentum. The BFKL kernel has the form: 



KBFKL{k'i, k'i; k[,k'^) = -N,g^5\k[ +k'^- k'l + k'i 



h,'2 ul/2 



"-l "-2 



m = 



(2^)3 fef fcf <52(fc; - fc'/) <52(fci - k'i) + , (5.19) 

(27r)3 P{l-kY ' 



a(fc^) = 1 + /3(A;) is the Regge trajectory of the gluon. For that reason /3 is also called, somewhat 
inaccurately, the trajectory function. 

The configuration-space representation of the BFKL equation is rather more complicated. The 
convolution of the BFKL kernel with a pomeron wave function can be expressed in terms of three of 
the five fundamental functions we will discuss later (see Section 5.3 and Equation 5.46). 



k'j k'2 



k'j k'2 






Figure 5.1: Graphical representation of the BFKL equation. 

The BFKL equation can also be presented in a graphical way, see Figure 5.1. The vertical curly 
lines with slanted lines through them represent reggeised gluons. The BFKL kernel in the second 
term on the right represents a superposition of exchanges of ordinary QCD gluons between the two 
reggeised ones. If one imagines solving the BFKL equation through iteration, it becomes obvious that 
this adds up to a superposition of ladders with an arbitrary number of rungs — the pomeron. While 
the s-channcl gluon which forms the rung is an ordinary QCD object, the vertices which couple it to 
the reggeised gluons are not. The BFKL kernel contains several terms for each vertex, which has the 
effect that Feynman graphs in which rungs cross are also included in the final superposition of ladders. 

The; solution of the BFKL equation was initially found only in the forward direction, ie for q ~ 0. 
Only when it was realised that the BFKL kernel is invariant under global conformal transformations 
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in position space, the non-forward solution was found. The consequence of this state of affairs is that 
the conformal eigenfunctions (5.8) are eigenfunctions of the kernel: 



dV'i Kbfkl{puP2;p[,p'2) S('^'")(p;,p^) = xM E^'''"\puP2) , (5.20) 
with the eigenvalue 



xiy-, n) = 



2^(l)-^(^+-V/''"'"' 



11^ 



(5.21) 



where ij^ is the logarithmic derivative of the F function. 

The solution of the BFKL equation can be written in terms of the conformal eigenfunctions: 



oo ^ 

1>uipi, P2; p'l, p'2) = / 



27r [Av^ + {n + If] [4j/2 + (n - 1)2] 

(5.22) 

/ d^poE('^^^\puP2;Po)E('^'^^*{p'^,p'^-po). 

This is in effect the propagator of the pomeron, more precisely of two reggeised gluons in a colour 
singlet state. 

When pomeron exchange occurs in a scattering process mediated by photons, the pomeron couples 
to a quark loop into which the photon fluctuates. It can be shown that the impact factor (which 
describes the coupling) of any number of reggeised gluons in a colour singlet state coupling to a quark 
loop can be written as a superposition of two-gluon state impact factors.^ Therefore it is sufficient to 
consider the coupling of two reggeised gluons to the quark loop. 

The amplitude of the pomeron can be obtained by convoluting an impact factor with the pomeron 
propagator: 

£>2(fei,fe2) = J d^k[d''k'2 4'°{k[,k'2)Mkuk2;k[,k'2). (5.23) 

The impact factors <p^ are usually symmetric under exchange of their arguments. Prom this. Equa- 
tion 5.22 and the properties of the conformal eigenfunctions, it follows that D2 is also symmetric. 
Furthermore, it vanishes when its configuration-space arguments are set to be equal, since the confor- 
mal eigenfunctions then vanish. Besides, D2 vanishes when one momentum argument is set to zero 
or one configuration-space coordinate is integrated out. To sum up all relations in momentum and 
position space: 

£>2(fel, fe2) = -D2(fe2, fel) r>2(pi, P2) = -D2(P2, Pi) 

d'^lD2{l,k-l) = D2{p,p)=0 (5.24) 



L'2(fe,0) = = ^2(0, fc) J d^P2D2{pi,P2)=0 = j d^PlD2{puP2) 

Strictly speaking, D2 is called the "amputated" amplitude because QCD gluon propagators have 

to be added to attach it to anything. The corresponding amplitude with gluon propagators is denoted 
by <l>2. There is a third type of amplitude called C2 which differs from D2 in that it does not contain 
the Reggeon propagator (w — /3(fei) — /3(fe2))~^, which is implicit in the BFKL equation. 



^Strictly speaking, one has to use a "cut" quark loop, in which quark lines are forced to be on shell. See for instance 
[Ew98] or [BE99]. 
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5.2.7 Useful formulas 

This section contains the derivation of some identities which wiU be used frequently later on. The first 
of these is the Fourier transformation of the gluon propagator in (two-dimensional) transverse space. 
The two-dimensional gluon propagator has to be regularised with a gluon mass in order to perform 
the Fourier integral: 

oo 27r 

1 / d^k—^e'''P = -^ I dk ^ , / d^e^'^l"!-^^ 
(27r)2 J k^ + rn^ (27r)2 J k^+rn^ J ^ 



oo 



This integration is performed in somewhat more detail at the end of Appendix A. 2. The integration 
formula by which the Bessel function Kq is obtained is cited there, in Equation A. 16.^ 

Since the argument of Kq contains the gluon mass m which will be set to zero in physical expressions, 
we are interested in the expansion of Ko for small arguments: 

Kq{z) ^\n2 + -4)(l)- hi z + 0{z^). (5.26) 

Here is the logarithmic derivative of the Gamma function and tp{\) — —7 « —0.5772 is the negative 
of Euler's constant. Since all mass-regularised expressions imply the physical limit of m — » 0, KQ{m\p\) 
will often be substituted by [In 2 + — Inm — In [p|] and vice versa. 

Many of the expression we will handle contain Laplace operators, some together with the regularised 
configuration-space propagator we just calculated. Therefore the Laplacian of ln|p| will be needed 
occasionally. In fact the logarithm is the Green's function of the Poisson equation in two dimensions: 

A ln|p| =27r^2(p). (5.27) 

To prove this, we look at the function ln(|p| -|- e). We calculate 

1 1 
IpI + e - IpI {\p\^) " M+ ~ \p\ (|p| + 



Aln(|p| + e) = V • = o^^^ - JTZTT^ = 1 „in „i , ,,2 ' (5-28) 



where is the unit vector in radial direction. The third expression was obtained by product dif- 
ferentiation and by using V • = j^. In the limit e ^ 0, the final expression diverges if p = 
and approaches zero otherwise. These are the values the Dirac delta function takes. To obtain the 
multiplicative constant in Equation 5.27, we integrate over a disc with radius R around the origin: 

R R 

J dVAln(|p| +e)=2w J drrAln{r + e) = '^^^ J * ^^T^y = ^tt ^- -^-^ + 1^ ^ 2w . 

disc 

This accounts for the constant 27r in (5.27). 

Equation 5.27 has an important implication for the derivatives of complex coordinates. To see it, 
we express the Laplacian in complex derivatives (see (5.6)). 

2w6'^{p) = A\n\p\=4dd*\n\p\=4dd*^ \n{pp* ) = 2d\ 

2 p 

By performing the derivation with respect to p first, we obtain 2d*^. Hence we have: 

d^ = d*^=nS\p). (5.29) 



*An observant reader will have iiotieed that the result is finite only for p ^ 0, ic it is still (logarithmically) ultraviolet 
divergent. This requires another regularisation if the propagator is integrated over. We will perform it later. 
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Put into words, this means that p and p* are not independent variables at the origin. Note that 
this relation holds only if the total power of the complex conjugated variable in the expression to be 
differentiated is —1. It is illegal to split off a factor 1/p* from an expression and use (5.29) in the 
product rule of differentiation. 



5.3 The five component functions in momentum space 

Having presented a selection of basic formulas, we will now start to investigate expressions closely 
related to the problems we intend to discuss. The main aim of this and the following chapter is to 
investigate a vertex describing the transition from two to six reggeised gluons. It was shown in [BE99] 
that such vertices can be expressed in terms of five fundamental functions. The starting point of 
these vertices is always a two-gluon amplitude. Therefore the amplitude D2 is already included in the 
definition of the five functions. Their momentum-space representation is the following: 

aik„ kM=l^ (f - k,nt fc, - fc,)2 ^^(^' ki + k. + k^- I) 



S{kuk„k,) = j (|^ f2(f\)2 J2(fcr+fc2,fc3) 



The function a is symmetric in its last two arguments. One can see this by substituting I ki+k2 + 
k^ I and using the fact that D2 is symmetric in its two arguments. All the functions vanish when 
their first argument vanishes, s and t also vanish when their last argument is zero, because the D2 
amplitude then vanishes according to the third formula of (5.24). 

There are some relations between the functions a, b and c and between s and t and the trajectory 
function /3. 

6(fei, fez) = a(fei, 0, fez) = a(fei, fez, 0) 
c(fe) = 6(fe,0) = a(fe,0,0) 

s(fei, fez, fea) = - P{ki)D2{ki + fez, fes) (^-31) 

i(fei, fez) = s(fei, 0, fez) = - /3(fei)£)z(fei, fez) 

The trajectory function /3, which was already introduced in Section 5.2.6, has the following form in 
momentum space: 

This integral can be performed after applying a regularisation scheme. The two regularisation schemes 

which have been used for this expression arc mass regularisation and dimensional regularisation using 
the dimension d = 2 + e. The derivations are lengthy and will not be performed here. The results are 
as follows: 

Mass regularisation: /?(fe) = lim ^ — ttIu — ^ (5.33) 

m—^o 2 m 
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Figure 5.2: Diagrams representing the five functions a, b, c, s and t with location of configuration 
space coordinates and momenta. Coordinates or momenta not appearing as arguments are integration 
variables. 

Dimensional regularisation: /3(fe) = lim — 27rA^c5^(fc^)' - (5.34) 

e— >o e 

The fimctions a, b and c represent real corrections, that is terms composed of them describe the 
production of on-shell gluons. s and t describe virtual corrections, self-energy corrections to the 
reggeised gluon propagator. 

Figure 5.2 shows a diagrammatic representation of the five functions. With a view to the following 
sections, it shows both the momenta of the gluon lines and the coordinates of the vertices. The momenta 
of the unlabelled lines are fixed by momentum conservation at the vertices. The loop momentum I and 
the coordinate po in the function t arc integration variables and therefore do not appear as arguments. 
Note that the middle momentum or spatial coordinate in the diagram for a is conventionally taken as 
the first argument of the function a. 

The formulas 5.30 can be obtained from the diagrams with a simple prescription [Ew98, BE99]: 

• Find the lowest vertex in the diagram. 

• Write down propagators for the lines attached to it from above. 

• Put the square of the sum of the momenta attached from below in the numerator. 

• Multiply with the amplitude D2 with suitable arguments. 

• Integrate over the loop momentum I and divide by (27r)^. 

Even though this is the simplest method, it is unfortunately somewhat asymmetric: Not all vertices 
and lines are treated equally. This has the consequence that when attaching a wave function to the 
lower legs, the full amplitude $2 has to be used even though the one attached to the upper legs is an 
amputated amplitude, D2. 

There is a more symmetric prescription for obtaining integrals from such graphs [BLW96]. Here 
all the lines get propagators and all vertices a square of the total momentum attached from above or 
below. As a consequence, amputated amplitudes D2 have to be attached both from above and below. 
However, this does not lead to the same functions a, b, c, s and t which are defined in the asymmetric 
fashion described above. 
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5.4 The five component functions in configuration space 
5.4.1 The real corrections a, b and c 

We will be interested in the properties of gluon vertices under the conformal transformations in trans- 
verse position space. For that purpose wc need to express thcni in position space. So wc have to 
Fourier transform the five functions a, b, c, s and t of which they are composed. Before we start with 
the function a, a remark concerning the wave functions is required. We will see that the configura- 
tion space functions arc expressed more conveniently in terms of the full amplitude $2 instead of the 
amputated amplitude D2. The two are related in the following way: 

Mkuk2) = ^^D2{ki,k2) (5.35) 
D2{ki,k2) = kl kl $2(fci, fc2) =klklj Spi d?p2 $2(pi, P2) e-'^'^'P'+'^-P-^ 
= J rfVi d^P2 $2(pi,P2)(-Ai)(-A2)e-^('=^''^+'==''=) 

d^pi (fp2 (Ai A2 $2(pi, P2)) e-^^'^'P'+'^^P^^ (5.36) 



The final result was arrived at by two-fold partial integration in each of the integration variables. 

Put into words, it means that the configuration space D2 amplitude is the double Laplacian of the 
configuration space $2- We will make frequent use of this identity since the Laplacians will enable 
us to perform the Fourier integrations easily and to derive the equivalents of the relations (5.31) in 
configuration space. 

Let us now turn to the Fourier transformation of the function a. As a first step, we substitute the 
expression (5.36) for D2: 

a(pi,p2,P3) = ^ j d^kicFk2cFk3 a(fei,fe2,fe3) e'^'''P'+'''P'+'''P''> 

r fc^ 

cPkid^k2d^k2.cPl d^p'^d^p'^ ^ 



• (Ai/ A2' $2(p'i,P2)) e'^''^''^"'"''^''^"'"''^''^"''''!"^*'^"'"*'^"'"*'^"'^''^) . 

Now we express the numerator k\ as a Laplacian which we can pull out of the integral since the variable 
pi appears nowhere but in the Fourier exponential. Besides we can perform the integration over fes, 
which gives a delta function, and then integrate over 



^ (-Ai) j d'^ki d''k2 d'ka dH j d^p\ d'p'^ 



'2 „' 



(27r)9 ^ ^' ] ' ' " » y ^ Kl". K2 _ fe2)2 _ fc^ _ fe2)2 

• (Ai' A2' <l?2(p'i,P2)) e*^'°^''^+*'^''^+'°^''='~''''i~^'°^+'°^+'°^~'^''2^ 

^jd-k,d-k2dH jd-p\ (,_fc^)2(/_fc^_fc^)2- 

• (Ai, A3$2(pi,P3)) e*('"i''i+''2''^"''''i"(''i+''^"'^''=*) 

• (Ai, A3 $2 (pi, Pa)) e^C^iP^+fesPsa+iPai') 
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After shifting the integration in ki and fe2 (fei := ki + k2 — I, k2 '■= ^2 — I): the momentum-space 
integral factorises into three separate integrals. Two of them are Fourier transformations of the two- 
dimensional gluon propagator which we regularise with a gluon mass (5.25), the third gives a delta 
function. 



(27r)7 \J kl 



= ^™ - 7;^TT^o(m|pi2|)-ft'o(m|pi3|) A2 A3$2(P2,P3) 

= lim - [ln2 - Inm - In |pi2|] [ln2 + - Inm - In |pi3|] A2 A3 $2(^2, Ps) 

m— »0 (iTTj'' 

This is the configuration space representation of the function a. Which of the last two expression to 
use is a matter of taste and convenience. 

The; transformation of the function h works quite similarly and will therefore not be shown in such 
detail. The steps are quite analogous up to the point where we shifted the momentum integration 
variables. In the case of the h function the integral does not factorise completely because there is one 
Fourier integration less. 

Kp1'/^2) = ... = ^/ <fk,dH I ^Vi ^2(f!fc^)2 (Ai^A2Mp;,P2)) e^C^^-^+^P-) 

Here we substituted fei := ki — I. We proceed by performing the Fourier integrations over the 
propagators with mass regularisation and perform a double partial integration in p'l, shifting the 
Laplace operator to one of the Kq Bessel functions. Here it proves useful that we use the amplitude 
$2 instead of D2. The surface terms vanish since the amplitude $2 vanishes at infinity. 

...= hm f d^p[Ko{m\pi2\)KQ{m\piv\)ArA2Mp'i,P2) 

= lim --^ / ciVi^o(m|pi2|)(Ai.[ln2 + V(l)-lnm-ln|pn-|]) A2$2(p'i,P2) 

m—^O [ZTT) J 

We have already expanded the Bessel function. As we saw in Section 5.2.7, the Laplacian of the 
logarithm is proportional to the two-dimensional delta function. The constants vanish when performing 
the diflferentiation. Hence we can finally perform the integral in p[ . 



1 

"0 (2"7)2 



6(pi, P2) = jim^ Aii4:o(m|pi2|) A2 $2(pi, P2) 



= lim --i-^ Ai[ln2-FV(l) -lnrn-ln|pi2|] A2$2(pi,P2) 

In the case of the function c, two configuration space integrations remain to be done after trans- 
forming the propagators. They are solved in the same way as for the b function, by shifting the Laplace 
operators from the amplitude to the Bessel functions and integrating out the resulting delta functions. 
So Laplacians act on both Bessel functions and turn them into delta functions which are integrated 
away. The result is just a Laplacian of $2 multiplied with a constant. Here are the position space 
representations of all three functions, with the Bessel functions expanded: 
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a{pi,p2,Pi) = - -^^^^ Ai[ln2 + '^(1) - Inm - In |pi2|] [ln2 + '^(l) - Inm - In |pi3|]A2A3$2(p2, Pa) 
b{pi,P2) = T-^Tg Ai[ln2 + V(l) -lnm-ln|pi2|] A2 4'2(pi,P2) (5.37) 
c(p) = -^A$2(p,p) 

The limit to ^ is not written explicitly any more but is always implied. The symmetry of the 
function a in its latter two arguments is evident in position space. 

The relations between the three functions (5.31) have analogues in configuration space. While in 
momentum space b could be obtained from a, and c from b, by setting one of the momenta to zero, in 
configuration space the corresponding spatial coordinate has to be integrated out. 



b{pi,P2) = J d'^p' a{pi,p',p2) = j d^ p' a{pi, p2, p' 
c{p) = J d'p'b{p,p') 



(5.38) 



This can be checked explicitly by performing the integrations analogously to the last steps in the 
Fourier transformation of b (or c). One Laplacian is shifted from the amplitude $2 to one of the Bessel 
functions. This results in a delta function which allows performing the integration. The factor 27r 
in (5.27) accounts for the different prefactors in a, b and c, and the minus sign in front of the logarithm 
in the expansion of Kq accounts for the alternating sign. Here is the explicit calculation for deriving 
c from 6: 

J d''p'b{p,p') = -^ J dVA[ln2 + ^(l)-lnm-ln|p-p'|]A'$2(p,p') 

= ^A J d^p' (A'[ln2 + ^(l)-lnTO-ln|p-p'|])$2(p,p') 
= ^A I d'p'{-2n6'{p-p'))Mp,p') 
= -^A$2(p,p) = c(p). 

5.4.2 The virtual corrections s and t and the trajectory function /3 

It is convenient to transform the expressions for s and t in terms of /? first and to perform the Fourier 
transformation of f3 separately. To this end, we insert the Fourier transform of the configuration-space 
trajectory function for the momentum-space (3 contained in s and t. Besides, we replace D2 by the 
double Laplacian of $2 as we did in the previous section. We get the expression for s: 

S(pi, P2, P3) = ^ ^ J d^ki d^k2 d^ks /3(fei) £)2(fel + fe2, fes) e'('=i''i+'==''=+'=3''3) 

(27r)6 AT,/ / d^kid^k2d^k3 J d^ p' d^ p[ d^ p'2 /3(p') (Ai, A2- $2(pi, P^)) • 



1 -2 

i(felpl+fe2p2+fe3p3-felp'-(fel+fe2)pi-fe3P2) 



• e 

Here we have already inserted the Fourier transform of the configuration-space amplitude $2- AH the 
momentum integrals give delta functions, which makes the spatial integrations trivial. The result is: 

2 

s(pi, P2, Pa) = - /3(pi2) A2 A3 $2(P2, Pa) • (5.39) 
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In the case of the t function, there are only two Fourier integrals and consequently one integral 
over a spatial coordinate remains in the final expression. Repeating the steps done for s, we get the 
following expression: 

t{p,,p2) = ... = ^ J d'ki j d'pd^po dip) Ao A2 $2(P0, P2) e'^^^P^-P-P'^^ 

= - J d'^po fi{pw) Ao A2 $2(po, P2) = y d^Po s(Pi, Po, P2) • (5.40) 



The relation between s and t (sec (5.31)) is already obvious from comparison with Equation 5.39. 

The main task remains to be performed, namely the Fourier transformation of the trajectory 
function (3. The techniques employed here are very similar to those used for the real corrections a, b 
and c. Wc express the in the numerator as a Laplace operator and factoriso the two-fold momentum- 
space integral through a shift in the integration variable. Both factors give us a two-dimensional gluon 
propagator in position space, a Kq Bessel function (in the mass regularisation scheme). 

= \imJ^^A{K,{m\p\)f = \\rn N,g^ j^[K,{m\p\)AK,{^^^^^ 
After expanding the Bessel functions we can perform the differentiations. Using (5.27) we obtain: 

/3(p) = Um ATeff^ ^ (-277 5^ (p) [ln2 + V(l) - In m - In |p|] + ^) (5.41) 

This formula contains two terms which are ultraviolet divergent under an integral and will later 
have to be regularised: 5{p) In p and l/\p\^. However, we will ignore this issue in most of this chapter. 
Wc will write down integrals over (3 as formal expressions. In section 5.5.3, the problems relating to 
regularisation will be discussed. 



5.4.3 Summary of Results 

This section contains a summary of our results so far: the configuration space representations of the 
five component functions. Figure 5.2 shows how the coordinates correspond to the location of the 
vertices in the diagrams. 

a(pi,p2,P3) = - -^^-^ Ai[\a2 + - Inm - ln|pi2|] [\a2 + ^{l) - Inm - In |pi3|] A2A3$2(P2, Ps) 

6(pi,P2) = 7;r^ Ai[ln2 + V(l) -lnm-ln|pi2|] A2 $2(pi,P2) 

{2-kY 

c(p) = -i-A$2(p,p) (5.42) 

s(Pi,P2,P3) = (27r52(pi2) [ln2-FV'(l) - Inm -In 1/912 1] - f— 72 ) A2 A3 $2(p2, Ps) 



i(pi,p2) = j d^po (^27r<52(pio)[ln2 + V(l)-lnm-ln|pio|]- 1^^^ AoA2$2(po,p2 

The limit m ^ is implied in all expressions even though it is not written there. 
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The relations between the five functions and the trajectory function have the following form in 
position space: 



b{pi,p2)= j d^p' a{pi,p',p2) = j d^p' a{pi,p2,p') 

c{p)= J d''p'b{p,p')= J d^p'd^p"a{p,p',p") 
2 

P2, P3) = - -Tf-^ P{P12) A2 A3 $2(P2, Ps) 



(5.43) 



t{pi,P2)= j d^ po s{pi, Pq, P2) = - j /3(pio) Ao A2$2(P0,P2) 

5.5 The building block of vertices: the function G 

5.5.1 The function G in momentum and configuration space 

It has already been mentioned that 2 ^ n reggciscd gluon vertices, one of which we will investigate 
in the following chapter, can be expressed in terms of the five functions which were the topic of the 
previous section. However, the vertices do not seem to be arbitrary combinations of them. The two 
known vertices (V2^4 and V2^e) arc sums of terms of one function which can in turn be expressed 
through the five fundamental ones [BW95, Va98, BV99, EwOl]. It is conjectured [EwOl] that higher 
vertices also can be expressed in terms of this function. The function is denoted by G and has the 
following form: 

G{ku k2, fes) = Y [2 c(fei + fe2 + fes) - 2 b{ki + fes, fes) - 2 6(fe2 + fes, fei) + 2 a(fe2, fei, fes) 

+ t{ki + k2, ka) + t{k2 + fes, fei) - s{k2, fei, fes) - s{k2, fes, fei)] • (5.44) 
To avoid writing so many sums of momenta, one can express this in a short-hand notation: 

G(l, 2, 3) = ^ [2 c(123) - 2 &(12, 3) - 2 &(23, 1) + 2 a(2, 1, 3) 
2 

+ i(12,3)+i(23,l)-s(2,l,3)-s(2,3,l)] . ' 

Each number stands for the momentum of the corresponding reggciscd gluon, that is 1 for fei, 2 for fc2 
and so on. A set of several numbers, for instance 123, stands for the sum of the momenta, fci + fc2 + fca. 
An empty argument, which will be denoted by "— ", stands for zero momentum. This notation has 
the advantage that it can be easily transformed into a configuration-space expression as well. In 
configuration space, a number stands for the corresponding coordinate. A set of numbers means that 
the coordinates are identified with a delta function with which the expression has to be multiplied. 
This follows from the identity 

^ J d2fc^rf2^^^(fc^^^^-)gi(fe,p.+fe.p.)^_l_ J d2fe/^2^^_^(fe,)gi(feip,+fe.(p.-pO)=j2(p^^-)_^(^^) 

which holds for an arbitrary function /. An empty argument "— " translates as an integral over the 
corresponding coordinate argument in position space. 

It can be easily seen from (5.45) that G is symmetric imdcr exchange of its first and third argument. 
The symmetry does not extend to the second argument. However, it has a different interesting property: 
When the second momentum argument is set to 0, the a and s functions vanish and G reduces to a 
derivative of the BFKL kernel multiplied with a constant. 

(5.46) 

G(l, 3) = y [2 c(13) - 2 6(1, 3) - 2 6(3, 1) + t{l, 3) + 1{3, 1)] = — Ai A3 (Kbfkl O $2)(1, 3) 
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Besides it can be proven that G vanishes when its first or last momentum argument is zero. This 
is because then the component functions either vanish or turn into a different component function 
with fewer arguments (according to (5.31)). They then cancel among each other. Here is the explicit 
calculation if the first argument vanishes: 

G(-, 2, 3) = y [2 c(23) - 2 6(2, 3) - 2 6(23, -) + 2 a(2, -, 3) 
+ i(2, 3) + 1(23, -) - s(2, -, 3) - s(2, 3, -)] 

= y [2 c(23) - 2 6(2, 3) - 2 c(23) + 2 6(2, 3) 
+ i(2,3) + 0-i(2,3)-0] 

= 0. 

For our investigation of 2 ^ n gluon vertices, we want to know how G transforms under conformal 
transformations in transverse position space. Therefore we need the configuration space representation 
of G. It can easily be obtained from the coordinate-free symbolic notation as described above. The 
configuration-space expression for G is: 

G{pi,P2,P3) = y [2(5^(P12)^^(P23)C(P2) - 2,5^(pi2)6(/9i,P3) - 2 S'^{p23)b{p3, Pi) + 2 a{p2, pi, ps) 

+ 6'^{pi2)t{pi,p3) + (5^(p23)i(P3, Pi) - s(P2, Pi, Pa) - s(p2, Pa, Pi)] . (5.48) 

5.5.2 The conformal transformation properties of the function G 
Initial simplifications and properties under translations, rotations and dilatations 

The function G is a component of expressions for 2 — > n reggeised gluon vertices which occur in unitarity 

corrections to BFKL. Since the BFKL kernel was found to be conformally invariant in transverse 
position space, it is a legitimate question to ask whether these vertices, too, are conformally invariant. 
This has been proven for the 2 — > 4 gluon vertex [BLW96]. As the higher 2 — > n gluon vertices contain 
a significantly larger number of terms of the five functions a, 6, c, s and t, to approach this question 
directly for them would be all but impossible. But since the largest known vertex, ¥2^5, and quite 
possibly higher vertices, are composed of G functions, their conformal transformation properties can 
be inferred from the properties of G. Therefore investigating G's conformal transformation properties 
is of great importance. 

Conformal n-point functions are defined as expectation values of a product of n conformal fields. 

By contrast, G and the gluon vertices are defined to contain only one conformal field $2 (which is 
a superposition of E^'^'^'^s). The expression we will identify with a conformal n-point function is a 
convolution of the 2 — > n gluon vertices with several conformal fields. Since the convolution integrals 
themselves transform under dilatations and the inversion, we will not expect the gluon vertices, and by 
implication G, to be strictly invariant. Rather, G should transform in such a way that it is invariant 
together with convolution integrals over its arguments. In the following, we will call this property 
conformal covariance. 

The amplitude $2 defined into G has to be exempted from our investigation. The transforma- 
tion properties of the conformal eigenfunctions are given exclusively by their external coordinates. 
This (or the corresponding transformation properties of general conformal fields) is what leads to the 
specific form for conformal n-point functions (5.13), (5.14) and (5.16). Hence we need not take the 
transformation properties of $2 into account when investigating G. 

To investigate how G transforms under conformal transformations, we start out from Equation 5.48. 
However, expressing G in terms of the five functions a, 6, c, s and t is of limited use for determining its 
conformal transformation properties. Of these functions, only c is conformally covariant on its own. 
Therefore we will rewrite G in a way which may seem far more complicated at first sight, but which 
will allow us to find conformally covariant groups of terms. 



(5.47) 
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Let us first group the terms in (5.48) somewhat differently: 



+ 5\pi2) {-2b + Pa) + S^ip23) (-26 + i)(p3, Pi) 
+ (2a(p2, Pi, Pa) - s(p2, pi, Pa) - s(p2, Pa, Pi))] ■ 



(5.49) 



It is well known that the combination —26 + t and the expression containing a and ,s in the last 
line are infrared finite. We will now simplify these two expressions. After integrating out the delta 
function in the first part of t, we obtain for —26 + t: 



(-26 + t)(pi,p3) = 



(27r)3 



(27r)3 



- 2(27r)Ai[ln2 + V(l) - Inm - In |pi3|] A3 $2(pi,P3) 

+ 2 (27r) [ln2 + ^(1) - Inm] Ai A3 $2(pi, P3) 

1 



! / d'po ( 



277 5 (pio)ln|pio| 



27r Ai In |pi3| A3 $2(pi, Pa) 



iPiol 



Ao A3 $2(P0,P3) 



27r(5^(pio)ln|pio| + -. 

\Pw\ 



(27r)3 



{2TTfS\pi3) A3 + 47r ^1 ^ ' A3 + 2n In jpial A1A3 
iPiar 



Ao A3 $2(P0,P3) 

^2(P1,P3) 



J rfVo (|27r52(pio)ln|pio| + 1^^^ AoA3$2(po,P3)^ . (5. 



50) 



The last line in (5.49) can also be simplified. Since the first Laplacian in the a function affects only 
the Kq Bessel functions, we can expand them and perform the differentiations with the product rule. 
The resulting terms can be combined with the terms from the s functions to form a concise result. 

2a(p2,Pl,P3) - s(P2,Pl,P3) - s(p2,P3,Pl) = 

2 



(27r)3 



- A2-fCo(m|pi2])-fl'o(m|p23l) - 27r5 (P12) -K'o(m|pi2]) - -j pj- 

IP12I 



27r 6^{p23) Ko{m\p23\ 



\P23\ 



(27r)3 



-2Tr6^{pi2)Ko{m\p23\) + 2 



Pl2nP23P 



Ai A3 $2(Pl,P3) 

2nS^{p23)Koim\pi2\) 



- 27r^2(pi2) Ko{m\pi2\) - 2Tr6^{p23) Ko{m\p23\) + 



1 



1 



|Pl2p |P23P 



Ai A3 $2(P1,P3) 



{2^r 



- 27r{S^pi2) - 5\p23)) {Ko{m\p^2\) - Ko{m\p23\)) 



+ 



P12 _|_ P23 



lPl2l^ ' |P23'^ 



Ai A3 $2(pi,P3) 



(27r)3 



27r(,52(pi2)-<52(p23)) In 



IP12I 

|P23| 



P12 



P23 



IP12P IP23P 



2n 



(5.51) 



Ai A3 $2(P1,P3) 



5.5. The building block of vertices: the function G 



71 



In the last step we have again used the expansion (5.26) of Kq for small arguments, assuming m to be 
small. 

We can now put everything together and write down the whole G function: 



G(P1,P2,P3) = 



(27r)3 



(27r)2(52(pi2)(52(p23)(-A2 $2(P2, P2) + (Ai + A3) $2(Pl, Pa)) 



P23 • Vi 



P12 • V3 



+ 47r(5^(pi2) ^f;^ 12 ' A3 $2(P1,P3) -47r(52(p23) ^^J^-j^ Ai$2(pi,P3) 



IP23I 



IP12I 



-5'^{pi2) j d^Po (^27r<52(pio)ln|pio| + 1^;^^ Aq A3 $2(po, P3) 

-'^^(p23) J d^po (^27r52(p3o)ln|p3o| + 1^^^ Aq Ai $2(po, Pi) 

+ 27r((5^(pi2) In |pi2| + (5^(p23) In |p23|) Ai A3 $2(pi, P3) 
2 



+ 



P12 _|_ P23 



IP12I 



IP23I 



Ai A3 $2(P1,P3) 



(5.52) 



The second term in the first line comes from the first term of (5.50). (Note that there are two terms 
containing —26 + t in (5.49).) The following three lines are from (5.50), too; the minus sign in the 
second line arises because wc have replaced P21 with pi2. The third term in (5.50) cancels part of the 
first term in (5.51), leaving the last two lines. 

At this point we can already say that G is invariant under translations and rotations (ignoring the 
amplitude $2)- This is obvious because G contains only relative vectors and differential operators. 
Both are invariant under translations. Furthermore, all expressions are rotationally invariant since the 
vectors occur only as moduli or scalar products. 

We can also already determine G"s transformation properties under dilatations. Let us first look 
at how elementary expressions of which G is composed change under dilatations: 



/ 



P 
V 

A ■ 

d^p 



5\P) 
In IpI 



Ap 

A-^A 

A^y d^p 

In IpI +lnA 



(5.53) 



Using those rules, it is easy to sec from (5.52) that all terms of G get a factor A^'' under dilatations 
(ignoring the transformation of <1>2). In addition, the terms containing a logarithm of a coordinate give 
an additional term In A. Taken together, these extra terms cancel: 

-S'^ipii) j d^po 27r5^(pio)lnAAo A3 $2(po, Ps) - '^^(P23) j d'^po 27r<5^(p3o) In A Aq Ai $2(po, Pi) 

+ 2n{S^{pi2) InA + S^{p23)ln\) Ai A3 $2(pi,p3) 

= 277 lnA( - (^^(pia) Ai A3 $2(pi,P3) - (52(p23) A3 Ai $2(P3,Pi) 

+ 5^(pi2) Ai A3 $2(P1, P3) + ^^(P23) Ai A3 $2(P1, P3)) 

= 0. (5.54) 
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This result is what we expected. It means that when we convolute G with conformal fields, the 
transformation properties of the resulting expression are given exclusively by the attached fields. The 
three convolution integrals provide a factor under dilatations which cancels the we obtained 
from G. Therefore the transformation properties of the whole expression under dilatations are given 
by moduli of the external coordinates of the fields (including the $2 contained in G), raised to a power 
given by the conformal weights. This is a property of conformal n-point functions. One can conclude 
that vertices constructed from G in the way presented in the following section are conformal four-point 
functions. It only remains to be shown that the same applies under inversion, which will be done in 
the following section. 



The complex representation and properties under inversion 



Unlike the other basic conformal transformations, the inversion cannot easily be expressed in the 
vectorial notation. Therefore we will rewrite G in terms of complex variables. The result looks very 
similar to (5.52). 



G{pi,p2,P3) 



(27r)3 



(27r)2,52^Pl2)(52(p23)(-A2 $2(P2,P2) + (Al + A3) <^2{pi, Pa)) 

Ai$2(/0l,P3) 



+ &TTd {P12)—. 7^ A3^2[Pi,P3) - ond {P23)—, TT- 

\P23\ |Pl2r 



-S^{pi2) J d^Pn (^27r<52(pio)ln|pio|+ 1^^^ Aq A3 $2(^0, Ps) 

-^^(P23) J d^po (^27rS^{p3o)in\p3o\ + Aq Ai $2(^0, Pi) 

+ 27r((5^(pi2) ln|pi2| + ^^(^23) In |p23|) Ai A3 $2(^1, P3) 
1 ' 



+ 

P12 P23 



Al A3 $2(P1,P3) 



(5.55) 



Now we can investigate how G changes under the inversion transformation. As we did above for 
dilatations, we first write down how simple expressions occurring in G transform; 



P 
P12 

d 
A 



_ P12 
P1P2 

H'A 

d'p \p\-' 



(5.56) 



5\P12) 

In 1/912 1 



\pi\'5\p,2) 

\n\p12\ - ln|piP2| 



Let us now look at the terms in (5.55) line by line. The terms in the first line all get a factor of a 
coordinate modulus to the power of twelve. Eight powers of coordinates come from the delta functions, 
the remaining four from the Laplacians. Since the delta functions identify all three coordinates, we 
can express this factor in terms of any of them. For reasons which will become clear shortly, we will 
choose to write it as |pi|^|p2|^|p3|^- 
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For the second line, we will make use of the explicit form of the derivative of the conformal 
cigcnfunction E^'^'"'^ . Since $2 is a linear combination of these eigenfunctions, it is sufficient to treat 
inversion of the expression for $2 = E'-'''"-\ The derivatives of E^"'""^ are calculated in Appendix C. 
For the first term in the second line of (5.55), we have: 

0* A3 i;(-.")(pi, p3) = 1 V1A3 £;('^'")(pi, p3) = {-l-n + 2iu) A3 E^-''^\puP3) ■ 

2 P1P13 

Leaving out the E^"'^^ and constant factors from the first term of the second line, we obtain the 
expression we have to investigate. Taking into account the definition of the scalar product (5.5) and 
the delta function, it can be simplified as follows: 



* * / I 12 ^3 = <> (P12) X h — -. 1 

PlPhJ \P23? 2 pU P23 



5\pi2) ( P23, r-^ A3 = 5\p,,) - ( ^ + ^ ) ^ A3 . (5.57) 



Under inversion, this becomes: 

I |4 x2/- \ 1 /^^l I Pl\ IP2p|P3p I |4 A I |2| |4| 16 r2/ \^ ( P'^ , Pl\ ^ A 

IP2I 5 (P12) ^ — + I 12 m\ A3 = \pl\ \P2\ \P3\ 5 (P12) ^ — + -T T— 12 A3 

2 VP3 P3/ lP23r 2 VP3 ^3/ |p23r 

|4l„ l4i |4 r2/ \ 1 /^^3 , 1 



2 vPi Pi/ lP23r 

The last expression is the same as (5.57) with the additional factor |pi |^|/92 I^IpsI^- It has been obtained 
by pulling a factor |pip/|p3p out of the scalar product. The second term in the second line of (5.55) 
can be treated analogously. 

Now come the most difficult terms, the integral terms in the third and fourth lino. Wc will start 
with the second half of the term in the third line, the term l/|/9iop in the parentheses. Under 
inversion, the delta function yields a factor the integral Ipol"^) the fraction |/9ip|poP) and the 

Laplacians IpoHpsI^- From our experience with the previous terms of the G function, wc expect to 
obtain |/9i|^|p2|^|p3|'' multiplied with the original term. Therefore, we rewrite the factors with the help 
of the delta function as: |pip|p2|^|p3i^ / • • • ipop • • •• What remains to be done is to convert the |poP 
in the integral to a To that end, we have to prove that 

/ dVodpiT- IpoH ^^AoA3$2(PO,P3) =0. 

J \Pw\ 
We write the difference down in a different way: 

\P\? - \po\^ = PipI - PoPo = PiPio + PwPo ■ 
For the first term, we obtain by writing out the first Laplacian and with partial integration: 

[ cFp^ PL 4ao5* A3 *2(po, P3) = - / rfVo fa* -^"j 4^0 A3 $2(P0, P3) ■ 
J pw J \ PwJ 

V ' 

=0 

Note that the derivative in parentheses is zero even though there is a po contained in the denominator 
of the fraction. Since partial integration is based on the differentiation of a product, relation (5.29) 
does not apply. The second term can be shown analogously to vanish under the integral. This tells us 
that we can replace the |po|^ under the integral by a Hence we get the same factor as for the 

previous lines, |/3i|''|p2|''|p3|^- 

The fraction term in the fourth line of (5.55) can be treated analogously. The logarithm terms in 
the third, fourth and fifth line also get a factor |pi |'*|/92|^|p3|^ (taking into account the delta functions). 
In addition, the logarithms give rise to additional terms according to (5.56). Because of the delta 
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functions, they can all be written as logarithms of |p2p- As was the case under dilatations, they cancel 
among each other (cf. (5.54)). 

That leaves the last term in G, the sixth line in (5.55). The square modulus takes the following 
form under inversion: 
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Together with the factors from the two Laplacians, this again gives a factor |/9i|^|p2nP3|'*- 

When one convolutes G with conformal fields, the factor |/0i|^|p2|^|p3|^ is exactly cancelled by a 



factor 



Pt 



from each convolution integral. It follows that G is conformally covariant in the sense 



defined above, and gluon vertices composed of G-functions are conformal n-point functions. 



5.5.3 Reguleirising G 

In the previous section, we have proven that G is conformally covariant, that is its conformal trans- 
formation properties are precisely the opposite of the properties of convolution integrals over its argu- 
ments. However, we derived this only for the unregularised, formal expression of G. This is sufficient 
for proving that the vertices constructed from it are conformal n-point functions. Since the vertices are 
physical quantities, any regularisation terms must cancel. Therefore the transformation properties of 
the vertices can be derived from the formal non-regularised expression for G that we have investigated. 
For instance, the proof of the conformal invariance of the (regularised) 2 — > 4 reggeised gluon vertex 
in [BLW96] shows explicitly that the extra terms from regularisation cancel. 

As a matter of interest, we will now look into whether and how regularisation changes G's con- 
formal transformation properties. In this section, we will always imply convolution integrals over G's 
arguments. We will see that this is necessary to treat all terms equally. 

As we have remarked in Section 5.4.2, the reggeised gluon trajectory function /3 contains two terms 
which are ultraviolet divergent under the integral of the function t. We regularise the fraction 
with a theta function 9{\p\ — e), ie a UV cutoff. The term (5^(p)ln|p| also constitutes a logarithmic 
divergence which we regularise by replacing In |p| with Ine. This causes these two divergences to cancel 
on Fourier transformation back into momentum space. f3 then takes the form: 

/3(p) = |^(-27r<52(p)[ln2 + V(l)-lnm-lne] + ^M_£)^ • (5.58) 

In the G function, the terms (5^(p)lne cancel, that is the first terms in the third and fourth line 
of (5.55) cancel with the fifth. The result is: 

{2TTfd^{pi2)S^{p23){-A2 <^2{p2, P2) + (Ai + A3) $2(^1, Ps)) 

+ 87r52(pi2) ^7^^ A3 $2(P1,P3) - &1^6\p23) ^1 *2(P1, Ps) 

\P23r \Pi2r 

-5\p,2) I rfVo^^^^ AoA3$2(po,P3) 

-S'{p23) /d^o ^^'f"',;'^ AoAi$2(po,Pl) 

J \P30\ 

^ f 9i\pi2\-e) 2 ei\p23\-e) \ .... . 

+ \ ^ 1 T ^ \ Ai A3$2(pi,P3) 

V IP12P IP12P23I IP23P / 
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(5.59) 
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The limit e — » is implied in this expression, even though we do not write it down explicitly. We split 
up the expression in the last line since not all of it requires regularisation. 

Let us now look at what changes under conformal transformations. The invariance under trans- 
lations and rotations is unchanged. Under dilatations, however, the arguments of the thcta functions 
change. The theta function in the third hne becomes 6'(|A/9io| — e) = ^'(|pio| — ^1^)- This gives rise to 
an additional logarithmic term according to the following identity: 

e 

= 2Trfipi) lim / d|pio|^ = 27rlnA/(pi)=27rlnA / d^ po d^p.o) f (po) . (5.60) 
e^oj Ipiol J 

e/\ 

We have used the abbreviation f(pa) = Aq A3 $2(/00: Ps); which is an analytic function in the vicinity 
of the point po = Pi- The domain of integration is limited by the theta functions to a ring around pi 
between the radii e/A and e. Since both are small quantities, we can approximate / by its value at pi. 

So we get two logarithmic terms from the third and fourth line in (5.59): — 27r lnA(5^(pi2) + 
i5^(P23))Ai A3 $2(pi, Pa)- If we imply an integral over the arguments of G (which we have whenever 
we use G), the derivation (5.60) applies also to the two regularised terms in the last line. Using the last 
expression in (5.60), we obtain the terms 2tt In A {S^{pi2) + S^{p23))Ai A3 $2(pi, Ps)- That is exactly 
the same with the opposite sign. 

The situation midcr inversion is more complicated. The |pio| in the thcta function (third line 
in (5.59)) transforms as |/Oio|/(|pi | |Po|)- The \pi\ is a constant with respect to the integration over po, 
so we can use relation (5.60) for A = l/|pi|. That we can do the same thing with the \po\ is not so 
clear. But we can argue that it can be replaced by a \pi\. 

The theta function 0(|pio|/|po| ~ e) restricts the domain of integration to the outside of the curve 
on which \piq\ = \po\e. Before the inversion, the corresponding curve was given by \pi — po\ = e. The 
Po on this curve differs only infinitesimally (of order e) from pi . Since the inversion is a continuous 
mapping, this still has to be the case after the inversion. Therefore 

^(|pio|/|po| - e) = ^(Ipiol - e Ipol) = ^(|pio| - e |pi| + 0{e^)) . 

This means that replacing \pq\ by \pi\ in the theta function only introduces an error of order e^, which 
we can ignore. We can then pull the pi out of the theta function as above, using formula (5.60).^ 

Hence we get an additional term 5^{pi2)4:Tr ln|pi| from the third line in (5.59), which came from 
a t function. The fourth line (also from a t function) gives an analogous term. Implying integrations 
over the arguments of G, the fifth line (which came from an s function) also gives two logarithms, but 
with a minus sign. As was the case for dilatations, the extra terms from s and t functions cancel. So 
we obtain the result that the regularised G function is conformally covariant. 

This is somewhat surprising given the fact that the regularised G contains a fixed scale, the UV 
cutoff e. One would have expected this to destroy covariance under dilatations and the inversion, 
unless the terms containing e themselves cancel. We argued at the beginning of this section that the 
terms containing Ine do indeed cancel. However, the fraction terms do not cancel. One can see this 
by writing the terms from the third and fourth line of (5.59) in a form similar to the fifth line. This 
requires convoluting the three arguments of G with conformal fields, which we will collectively denote 
by $3. We can then integrate out the delta function and rename the integration variable po ~* Pi- 



^The authors of [BLW96] use an analogous procedure, without remarking on why it is legitimate. 
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- [ <fpi(fp2d'^p3^3{pl,p2,p3)S^{pl2) [ rfVo^%^2-^ Aq A3 $2(po, Ps) = 

J J \Pw\ 

= - [ d^p2d^P3<^3{p2,p2,P3) [ rfVo^^^o-^ Aq A3 $2 (po, Ps) 
J J \P20\ 

= - f d'^pid^p2d^P3^3{p2,p2,P3) ^^^^^"^^^ Ai A3$2(P1,P3) 

This is almost but not quite the negative of the first term from the fifth fine, convoluted with $3. The 
fifth-line term would lead to the same integral, but with the first argument of $3 being pi, not p2- 

So why does the presence of the fixed scale e not destroy scale invariance? The answer is that, 
despite all appearances to the contrary, e does not constitute a fixed scale. Changing e gives rise to 
the same type of logarithmic terms as a dilatation transformation, according to (5.60). These cancel 
out among the difl[erent terms of G, as we have seen. So the regularised expression for G is in fact 
independent of the choice of e. It just cannot be written down without it. 



Chapter 6 

The 1 — ^3 pomeron vertices 



6.1 The irreducible 2^6 reggeised gluon vertex 

As explained in Section 1.6, the Extended Generalised Leading Logarithmic Approximation to BFKL 
involves amplitudes of more than two reggeised gluons in a colour singlet state. The six-gluon amplitude 
is a sum of a reggeising part, a partly reggeising part and a new irreducible part. The sum of the latter 
two parts fulfils an integral equation which was first derived in [BE99]: 



^ - E '^(^) I (^e'"" + D'/) «^«^«3a4aaaa 2, 3, 4, 5, 6) = 



= (y2'!\."6'"'°'"'"'-D2)(l, 2, 3; 4, 5, 6) 

+ ^ ^ /aia203 f ataman 2, 3; 4, 5, 6) 

+ E rf"^"^"^'"* ^aaaa ^(1. 2, 3, 4; 5, 6) (6.1) 
+ J(l,2,3,4;5,6) 

I R II 

Dg denotes the partly reggeising six-gluon amplitude, Dg' the completely irreducible part. The 

arguments of the spatial functions arc written as the indices of the reggeised gluons instead of the 
momenta. So 1 stands for ki and so on. A set of several indices (eg 134) for one argument would mean 
a sum of momenta. This notation has the advantage that it can easily be transcribed into configuration 
space. There an index stands for the corresponding coordinate, and a set of indices means that the 
coordinates are identified with a delta function. This correspondence can be derived by performing 
the Fourier transform of a function taking a sum of momenta as its argument (see Section 5.5.1). An 
argument without any index, denoted by a hyphen "— ", would mean a zero argument in momentum 
space and an integral over the respective coordinate in position space. 

The first term on the right hand side of (6.1) is a 2 ^ 6 (reggeised) gluon vertex which is local in 
rapidity. Since it cannot be expressed in terms of smaller vertices, it constitutes a new element of the 
conjectured conformal field theory of unitarity corrections. It has the following form: 

^ya^a,a,a,a,a,j^^^^^^ 2, 3; 4, 5, 6) = da,a,a,da,a,a, (W^i^2)(l, 2, 3; 4, 5, 6) 

a2 04 
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The colour structure is given by the symmetric structure constants dabc of 8u{2>). The sum runs over 
all ten possibilities of dividing six momenta and corresponding colour indices into two groups of three. 
The function {WD2) is a convolution of a function W and a pomeron amplitude D2. The amplitude 
D2 represents the pomeron attached to the vertex from above. Since this vertex occurs in unitarity 
corrections to BFKL, there is always a pomeron coming from above, and we will always consider 
the combination {WD2) rather than W alone. ^ {WD2) can be expressed in terms of the function 
G introduced in Section 5.5 or in (quite many) terms of the functions a, 6, c, s and t. Here is its 
representation in terms of G: 

(W£)2)(l,2,3;4,5,6) = |- ^ (-1)*^ G(123\M,M,456\M) . (6.3) 

MeP({i,...,6}) 

This sum runs over all subsets M of the set of indices. The notation "123\M" means the indices 1, 
2 and 3 except those contained in M. #M is the number of indices contained in M. This sum and 
the symmetry of the function G under exchange of its first and third argument cause the function 
(VFD2) to be symmetric under all permutations of its first three and its last three arguments, and 
under exchange of the two groups of arguments. Together with the sum in (6.2), this makes the 2^6 
vertex symmetric under arbitrary permutations of the gluons. 

The representation (6.3) is easily generalised to a different (even) number of reggeised gluons. In 
fact the 2^4 reggeised gluon vertex can be written as a sum over functions (FD2) which are defined 
analogously to (T4^r>2) [BW95]. The sum runs over all (three) possibilities of dividing four indices into 
two groups of two. It is conjectured that there exist higher 2 ^ 2n reggeised gluon vertices which arc 
also new elements of the field theory of unitarity corrections and which are of analogous form [EwOl]. 

6.2 The 1 — ^3 pomeron vertex from the irreducible 2 — ^6 
gluon vertex 

6.2.1 Projecting the 2 — > 6 reggeised gluon vertex 

Since the pomeron is a bound state of two reggeised gluons, a 1 3 pomeron vertex may be obtainable 
by projection from the 2^6 reggeised gluon vertex. This is analogous to the projection of the 2^4 
reggeised gluon vertex to obtain a 1 — > 2 pomeron vertex [Lot96, Ko99]. It is not clear a priori 
that the projection is not zero, in which case no pomeron vertex can be derived from the 2^6 
gluon vertex. The projection is performed by convoluting the spatial part of the gluon vertex with 
three pomeron amplitudes and by contracting the colour tensors with three Kronecker deltas with 
corresponding indices. Because of the symmetry of the vertex, it does not matter which pomeron has 
which arguments resp. colour indices. We will choose one pomeron to have the arguments 1 and 2, the 
second 3 and 4 and the last 5 and 6. 

These three pomerons shall have the external coordinates pa, Pb and pc, respectively. The fourth 
pomeron, whose amplitude has been defined into the function G and its component functions, shall 
have external coordinate pa- 

We will restrict ourselves to the ground state of the pomeron. This amounts to setting the full 
amplitude $2 = E^'^'^\ both for the three pomerons attached to the vertex from below and the one 
contained in G. We will not write down the external coordinates as arguments to save space; which 
amplitude is which will be clear from the conformal dimensions Va, ^6, and Ud. 

The pomeron with index d is special in that a corresponding amplitude is already defined into the 
vertex V2^q and its component functions. This amplitude is a D2 which contains a pomeron impact 
factor (see Equation 5.23). However, the principle of Regge factorisation states that impact factors are 

^So, strictly speaking, this is not a 2 — » 6 gluon vertex but a 1 pomeron 6 gluon vertex. However, we will follow 
the nomenclature in [Ew98, BE99]. 
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independent of the evolution of the amplitude rapidity. This allows us to discard the impact factor and 
attach an amplitude E^^'^^ , thereby obtaining a pomeron vertex which is independent of the impact 
factor. This E^'^''^^ corresponds not to D2 but to the non- amputated amplitude $2- 
Hence our projection of the vertex can be written in the following way: 



— >3P — ^aia2 ^asa^ ^a^ae 



(yai^a2a3a4a506£,2)(l, 2, 3; 4, 5, 6) 



<S>2//E'."d.o) 



E 



6.2.2 The colour structure 

The colour structure of the pomeron is 5ab- Because the symmetric structure constants date vanish 
when two of their indices are contracted, this has the consequence that less than half of the terms in 
the sum in (6.2) contribute to the 1^3 pomeron vertex. The only remaining terms are those where 
the two colour indices of each pomeron are contracted with different d tensors. Since we chose as 
the pomerons' arguments (12), (34) and (56), this leaves the following four permutations of (colour, 
momentum or coordinate) indices: 

(1,3, 5; 2, 4, 6), (1, 3, 6; 2, 4, 5), (1, 4, 5; 2, 3, 6), and (1, 4, 6; 2, 3, 5) . (6.4) 

Each of these permutations gives rise to the same colour factor because of the symmetry of date- 
Since each date tensor is contracted with one index of each of the ^s, the result is the contraction of 
two d tensors. Explicitly for the first permutation: 

40 

daia3a5da2aias^aia2^a3ai^ac,ae — daiaaa^daiaaa^ — • 



6.2.3 The spatial part 

Introduction 

The lion's share of the work to be done for the projection is sorting out all the terms of the spatial part 
of V2^6- This work is complicated by the fact that the 1 — > 3 pomeron vertex cannot be expressed 
in terms of G. Because E'^'^'^^ (p, p) — 0, some of the component functions of G vanish on projection, 
others don't. 

The brute-force way of extracting those that remain would be to write out {WD2) in terms of 
the functions a, 6, c, s and t and cancel those that vanish. This is indeed how it was first done, but 
there is a more systematic way. First, we observe that the four permutations (6.4) remaining after the 
contraction of colour tensors all lead to the same term. The reason for this is that they differ only by 
the exchange of two coordinates of the same pomeron wave function, which is symmetric. Hence it is 
sufficient to consider only one of the permutations and multiply the result by four. 

The second important observation is that the conformal eigenfunctions E^"'^") vanish if their two 
coordinate arguments are equal. Therefore terms in which both argument indices of the same pomeron 
occur in the same argument of G vanish under projection. This is the key to finding out which terms 
remain after projection. We now have a look at the first permutation in (6.4). The spatial part of the 
vertex for this permutation has the following form: 

(W^£>2)(1,3,5;2,4,6) = ^ ^ (-1)#^ G(135\M, M, 246\M) . (6.5) 

Mep({i,...,6}) 

We can see immediately that a term from this sum vanishes if two coordinates of the same pomeron are 
contained in the middle argument of G. The first and last argument of G contain only coordinates or 
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momenta belonging to different pomerons and hence cannot cause the term to vanish under projection. 
Therefore all the terms for which M does not contain both indices of a pomeron remain: 

(W£>2)p^3p(l,3,5;2,4,6) = |- ^ G(135\M, M, 246\M) . (6.6) 

M6-P({1 6}) 

MJ){1,2}AM21{3,4}AM2{5,6} 

To say it in words: M may contain cither 1 or 2 or none of these two indices, and 3 or 4 or neither of 
them, and 5 or 6 or neither. This yields 3 • 3 • 3 = 27 terms. However, two of these vanish because of a 
property of the G function: It vanishes when its first or third argument is empty, see Equation 5.47. 
This kills the terms for M = {1, 3, 5} and M = {2, 4, 6}. Hence 25 terms remain. 

The function G is composed of functions with one, two and with three arguments. The function 
with one argument (c) does not contribute at all to the pomeron vertex. Since its argument always 
contains all indices, c does not survive the projection. The functions with two arguments, b and t, 
are also components of the BFKL kernel. We will discuss them in the next part of this section. The 
functions with three arguments, a and s, will be dealt with after that. 

The BFKL term 

The functions h and t survive only if they are convoluted with pomerons in one specific way, namely 
if cac;h pomeron's arguments arc convoluted with different arguments of b or t. This means that both 
arguments of 6 resp. t contain three indices, one from each pomeron. First of all, this is the case for 
the term in (6.3) where M = 0. 

Besides, there are terms in the sum (6.6) where M contains only indices from the first or only 
from the second group of arguments of iWD^). Since one argument of b and t contains two arguments 
of G taken together (see (5.45)), one of the two b and t functions would have momenta of different 
pomerons in each argument and hence would contribute. However, by replacing M by its complement 
with respect to the group of arguments of {WD2) of which it is a subset, an identical term with 
opposite sign can be obtained (unless M is identical to the whole group of arguments, but these terms 
vanish, see above). Therefore all such terms cancel out. 

To illustrate this, let us look at an example. We again consider the first of the permutations in (6.4) 
with the spatial function (6.5). For example, M = {1,3} contains only indices from the first three 
arguments of {WD2). The corresponding term in the sum contains the following G function: 

G(5, 13, 246) = Y [2 c(123456) - 2 6(135, 246) - 2 6(12346, 5) + 2 a(13, 5, 246) 
+ 1(135, 246) + 1(12346, 5) - s(13, 5, 246) - s(13, 246, 5)] . 

The terms with more than three indices in one argument, such as c and the second 6 and t functions, 
vanish because this identifies two coordinates of a pomeron. For now, we ignore the functions a and s; 
we will deal with them in the next section. The remaining terms, —2 6(135, 246) and f(135, 246) do not 
vanish. However, the same terms occur in the term for the complement of M, M' = {1,3, 5}\M = {5}. 
The corresponding G function is: 

G(13, 5, 246) = y [2 c(123456) - 2 6(135, 246) - 2 6(2456, 13) + 2 a(5, 13, 246) 
+ 1(135, 246) + 1(2456, 13) - s(5, 13, 246) - s(5, 246, 13)] . 

Again, —2 6(135,246) and i(135,246) remain. But this term has the opposite sign because of the 
prefactor (—1)*^ in (6.6). So the t and 6 functions for M ^ cancel out.^ 

^The deeper reason for this is that 6/2 = 3 is odd. It can be conjectured that the same mechanism of cancellation 
occurs in the projection of higher 1 — > n pomeron vertices from 2 — > 2n reggeised gluon vertices for odd n. 
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Figure 6.1: Graphical representation of the part of the 1 — > 3 pomeron vertex which contains the 
BFKL kernel. The half-circular lines represent pomeron wave functions. 



Only for M = the functions b and t contribute to the 1^3 pomeron vertex. Knowing that 
G becomes a derivative of the BFKL kernel when its middle momentum argument vanishes (Equa- 
tion 5.46); wc can immediately write out this term. Since there is one such term for each of the four 
permutations whose colour structure does not vanish, we have to multiply the result by 4. We obtain 
the following BFKL term: 

• ^ Ai A2 {Kbfkl ® ^('^<^'°))(pi,P2) 



3 16 (27r)3 ''^ ^ ^ 



\P12\ 



\P12\ 



\Pla\\p2a\ 



\Pl2\ 



\Plb\\P2b\ 



\Pl2\ 



\Plc\\P2c\ 



\P12\ 



l/'ld||P2<i| 



(6.7) 

l+2ivi 



We have evaluated the convolution with the BFKL kernel using the BFKL eigenvalue (5.21) and put 
in the derivative of the wave function (C.6). 

In the last step, we have rewritten the BFKL eigenvalue x ^^s a function ^ which was already used 
in [Lot96]. It differs from x only in that it docs not contain factors of Nc and the coupling constant: 

Stt^ 



^{v) = 27r 



2^(1) -V 



IV 



Ncg- 



■X(2^,0). 



(6.8) 



Figure 6.1 shows a graphical representation of this part of the vertex. 



The OL terms 

What remains now arc the functions with three arguments, a and s. Since they have the same properties 
(both vanish only when their first argument, ic the second argument of G. is empty), they always occur 
together, and the result of the projection can be expressed in terms of the function 

a(2, 1, 3) := 2a(2, 1, 3) - s(2, 1, 3) - s(2, 3, 1) . (6.9) 

a is symmetric in its last two arguments. Figure 6.2 shows a graphical representation of this equation. 
The middle leg of a, which corresponds to its first argument, is marked to remind you of the fact that 
it is symmetric in the other two arguments, but not this one. 
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123 123 123 123 



Figure 6.2: Graphical representation of Equation 6.9. The arcs represent pomeron wave functions 
defined into the functions a, s and a. They were drawn as eUipses labelled "D2" in Figure 5.2 but we 
will drop that here since we will have many more of them. 



r\ r\ r\ r\ 




a(l,246,35) a(13,5,246) a(13,25,46) a(136,24,5) 



Figure 6.3: Graphical representation of the classes of terms with the a function. One representative 
term is printed underneath each graph (assuming the pomerons (12), (34) and (56)). The half-circular 
lines represent pomeron wave functions. 

We will now divide the remaining terms from (6.6) into groups that lead to the same type of integral 
after convolution with the wave functions. There are 24 terms since the case M = has already been 
dealt with above. The most obvious classification criterion is the number of elements M contains, ie 
the number of momenta or coordinates in the first argument of a. There are six terms for which M 
has one element, twelve for which it has two and six with three elements. (Normally there would be 
eight in the last class, but for two of them M contains a whole group of arguments and G vanishes, see 
above.) It makes sense that the numbers of terms in all classes arc divisible by six. Six is the number 
of permutations of three objects. Because of the high symmetry of the 2^6 reggeised gluon vertex, 
there are terms with every pomeron coupling to every pair of arguments of the function G. Terms 
which differ only by permutations of the labels (= external coordinates) of the pomerons belong to the 
same class since they lead to the same type of integral. The terms from different classes contribute 
with different signs due to the factor (— 1)**'^ in (6.5). The class in which M contains two elements 
has positive sign, the others negative sign. 

The class in which M contains two elements can be subdivided further. It contains those terms for 
which the first argument of a contains two indices. It makes a difference whether the second and third 
arguments also contain two or whether one has three and the other only one. With this subdivision, we 
have four classes in total, with six terms each. There is an intuitive interpretation of these classes which 
is presented graphically in Figure 6.3. The classes differ in how many pomerons are attached to each 
leg of the a function (ie how many are convoluted with the corresponding coordinate in position space; 
or how many there are whose momentum occurs in the corresponding sum of momenta in momentum 
space). 

Bearing in mind that a is symmetric in its last two arguments but not its first, it is easy to see 
that the four classes are indeed disjoint. What is more, they are also complete. It is impossible to 
construct a term which does not belong to one of the classes (taking the properties of a and the 
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pomeron amplitude into account). 



A closer look at the ce terms 



Before writing down the integrals containing a which make up the largest part of the 1 — > 3 pomeron 
vertex, let us recall the configuration space form of a which was already calculated in the previous 
chapter (Equation 5.51). In complex notation, it can be written as: 



a(P2,Pi,P3) 



(27r)3 



27r(j2(pi2)-(52(p23))ln 



|P12| 
|P23| 



Pl3 



P12P23 



Ai A3$2(P1,P3). (6.10) 



The two delta functions identify the second or third coordinate argument of G with its first. When 
this function is convoluted with further pomeron wave functions as illustrated in Figure 6.3, one or both 
of the delta function terms vanish after integration because the delta functions identify two coordinates 
of a pomeron. For the latter two groups both delta function terms vanish (since there are pomerons 
between the middle and both outer legs of a), for the others only one. 

Several terms of the function a are potentially divergent under an integral. The terms 
S^{pi2) ln|pi2| are obviously dangerous, but also the fraction may cause a logarithmic divergence 
with respect to the integration over pi2 or p23- The pomeron wave function for zero conformal spin 
E^'^'^\pi, P2) contains one power of \pi2\ and can hence regularise these divergences. Eliminating the 
delta functions requires that a pomeron is attached between the middle and one of the outer legs of a. 
Therefore the terms from the last two groups in Figure 6.3 are finite while in the others one divergence 
remains and has to be regularised. 

We will now write down the integrals resulting from each of the four groups of terms. We will 
denote them (and the groups) with an upper index equal to #M. the number of elements in M. The 
two groups with #M = 2 will be distinguished by the minimal number of coordinates in the other 
arguments of the G function and denoted by "(2,1)" and "(2,2)". The divergent terms are regularised 
in the same way as the function G in Section 5.5.3. The logarithm next to the delta function is replaced 
by Ine, and the fraction is regularised with a theta function. For the first group (#M = 1), we get: 



(1) 



P— >3P 
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3 16 (27r)3 
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27r<52(pi3) In^ - + 



P12I IP12PI/O13P 



^(|pi3|-e) 



A2A3£;("<^'°)(P2,P3). 



The fraction 40/3 is the colour factor, the factor 4 the combinatorial factor from the four permuta- 
tions (6.4). The factor 5^/16 is the product of the prefactors of {WD2) and G. We have rewritten 
the sum over the six terms of this group as a sum over the permutations of the pomeron indices which 
is equivalent to it by way of renaming the integration variables. It is denoted by Xia^b^c- ^^"^ 
now insert the explicit form of the pomeron wave functions. This allows us to evaluate the double 
Laplacian which was calculated in Appendix C (Equation C.6). 
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(6.11) 
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We have converted a logarithm into a factor inside the theta function. We had done the same thing 
in Section 5.5.3 when deriving the conformal transformation properties of the regularised G function. 
Here we have applied relation (5.60) with respect to the integration over pi, for A = I/IP12I and 
/(pi) = iJ^"" °) (pi,p2). As we argued in Section 5.5.3, (5.60) can be applied even when the factor A 

contains the integration variable. 

It is easily shown that V^^^3p alone has the transformation properties of a conformal four-point 
function. Invariance under translations and rotations is trivial. Invariance of the integral operator 

under dilatation with a factor A is not hard to prove: The three integrals give a factor A^, the 
denominator IP23I'' a factor A^**. Both the delta function and the fraction in the rectangular brackets 
give a factor A"^.^ 

There remains inversion. The first two integrals together with the denominator are invariant. So 
is the pi integral together with the delta function and the fraction. However, the theta function gives 
rise to an extra logarithmic term inside the brackets. Using relation (5.60) with respect to the pi 
integration for A = I/92//33I, it is computed to |p23p/|/5i2p • 27r(5(pi3) In |/32//53|. After performing the 
pi integration, everything except this logarithm is symmetric under the exchange p2 P3- Hence the 
additional term is overall antisymmetric and vanishes under the integration over p2 and ps. 

The integrals resulting from the other groups of terms are derived analogously. In all cases, the 
remaining integration variables are pi , p2 and p4 , which we will rename p3 . The second group, for which 
#M = 2 and one of the arguments of G contains only one coordinate, also requires regularisation. We 
deal with that as shown for Vpj23p, including absorbing the logarithm into the theta function. Like 

^^-'i3P) the vertices belonging to the other groups are also conformal four-point functions on their own. 
Here they are: 
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(6.12) 
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(6.13) 



■^Thc theta function is invariant under dilatations. This was tlic main reason for absorbing the logarithm into the 
theta function. In our derivation of the transformation properties of the regularised G function, we did not do that 
because the logarithm terms had cancelled with others from the component functions b and t. We aimed at the simplest 
possible expression for G. 
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(6.14) 



6.2.4 Simplifying the spatial peirt and the function * 
The 1 — >■ 3 pomeron vertex as a conformal four-point function 

In the previous section we have learnt that each of the five terms of the 1 — > 3 pomeron vertex is 
a conformal four-point function. This means that they are of the form (5.16). One small difference 
arises from the fact that we used complex conjugated wave functions for the pomerons attached to the 
vertex from below. This is equivalent to replacing u —u and n —n. Besides, the formula (5.16) 
can be simplified somewhat because we use only ground state wave functions with n = 0, and therefore 
h = h = iv. The resulting formula for our four-point function reads: 

E('^»'°)*(p„) E('''"°)*(p(,) E(''-°)*(pe) E^^^^'^Xpa)) = <l>{x,x*; {ui}) JJ K'^^-^^+i >>.) , 
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where i/i is — for i = a,h,c and Vd for i = d. The single argument of the conformal eigenfunctions 
is the external coordinate of the pomeron state. We have omitted the coordinates of the reggeised 
gluons contained in the pomerons since they are integrated over to obtain the expectation value on 
the left-hand side. By writing down all the factors, one can obtain the following explicit form of the 
four-point function: 
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Now we rewrite the ratios in brackets as follows (shown here for the first expression): 
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The first of the anharmonic ratios appearing here is what we defined as x in Equation 5.17, the second 
is x/{l — x). For rewriting the next two factors in (6.16) we use (6.17) with the indices a, b and c 
permutated in order: a ^ b ^ c. This has the consequence that the constant exponents — i of the 
first three anharmonic ratios cancel. From the last factor, we pull out the two anharmonic ratios which 
have padpbc in the numerator, so that the variable pa does not occur any more in the remaining part. 
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Pcd\ ) \ \Pad\ ) \ \Pbd\ J \ \Pb. 

To obtain the \E' function corresponding to our four-point function, we solve this equation for ^. We 
get the formula: 

\Pac\\Pad\ \ / \Pab\\Pbd\ \ ( \Pbc\\Pcd\ \ ( \Pbd\\Pcd\ \ ^ 

\Pcd\ ) V \Pad\ ) \ \pbd\ ) \ \pbc\ ) ■ ^ - ) 

From now on we will not write the arguments of ^ explicitly any more. It depends on the Vi as well 

as on X and x*. 

The BFKL term 

Since each of the integral expressions derived in Section 6.2.3 is a conformal four-point function, the 
function of the vertex can be written as a sum over five 4* functions associated with the BFKL term 
and the four a terms. They can be computed by inserting the integrals into Equation 6.18. We will 
start by calculating the ^ function of the BFKL term, ^bfkl- Inserting (6.7) into the formula for ^, 
we get: 
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Note that all the factors in (6.19) are conformally invariant by themselves. The function S is defined 
as the integral multiplied with the powers of the anharmonic ratios. It was defined with the arguments 
—fa, —Vb a-nd —Uc since in our case the first three wave functions arc complex conjugated. This is not 
the case in a general four-point function, so this notation facilitates comparison and generalisation of 
our result. 

The function S is symmetric under simultaneous exchange of the pomeron coordinates and the i) s 
{^a/h/c = ~i^a/b/cTVd = i^d)- This foUows from the fact that both the part of the vertex without the 
^ function (see (6.15)) and the integral of the BFKL term (6.7) have this symmetry. Note that a 
permutation of the coordinates changes the anharmonic ratio which is the last argument of S. Which 
anharmonic ratio results can be seen with the help of Equation 5.17. 
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The a terms 



In this section we will calculate the ^ functions corresponding to the terms containing the function a. 
We will see that they can be written in a form quite similar to "^bfkl- 

We will start with the term which arose from the set M having one element, V^^\ It contained a 
sub-integral with an expression which had to be regularised. The same integral has already occurred 
in a term of the 1—^2 pomeron vertex. The associated integral operator is conformally invariant and 
was found to have the eigenvalue 1/2 ^{u) [Lot96]. 
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Inserting this relation, we get for this part of the vertex: 
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Since the integrand except for the ^ function is symmetric in the pomeron indices a, h and c, we could 
remove the sum over the permutations and insert a sum over ^s instead. The result closely resembles 
the BFKL term (6.7), except for the sign and the argument of the ^ function. Therefore we can 
immediately write down the ^ function of this part of the vertex. 
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(6.22) 



The next group of a terms, V'^^'-'^) contains a new integral. We will see later that this integral 
occurs in none of the other terms. Therefore we will just write down ^'(2,1) without defining a function 
corresponding to S. Plugging (6.12) into (6.18), we get: 

^^2,1) = '^^Aj^{Avl + lf |a;|-3 + ii(-2''a+'^!.+^c->'<i) |i_^|i + §iK-2!/6+i/c+2,.rf) . 
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There remain the integrals y^^'^) and V'^^\ Looking at the formulas (6.13) and (6.14), one can see 
that they are very similar. One has to be aware that the only distinction between different terms is 
the number of pomerons attached to a given pair of coordinates, (12), (13) or (23). Which of the three 
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pomerons a, b and c is attached to which does not matter since all permutations are summed up. Also, 
the indices of the integration variables are immaterial because renaming them changes nothing. Only 
one pair of indices is distinguished from the others by the fact that the upper pomeron with index d 

is attached to it. 

Taking all this into account, the only difference between jj^j^^j ^(3) ^j-^^, following: In y(2,2)^ 

each of the three pomerons attached to the vertex from below has a different pair of coordinates as 
their arguments. In V'^'^\ none of the lower pomerons has the same pair of coordinate arguments as 
the upper (d) pomeron. Instead, two of them have the same pair of arguments. (This can also be 
seen in the graphical representation in Figure 6.3 when one realises that the upper two lines have the 
same coordinates as the outer lines below.) This implies that V'^^^ can be obtained from y(2,2) j-^y 
exchanging pd with the external coordinate of one of the two lower pomerons which have the same 
arguments {ph or pc in (6.14)), and simultaneously exchanging ^ —Vb or Vd ^ —i^c, respectively. 
Then the pomeron b resp. c has the arguments which previously none of the lower pomerons had, 
and pomeron d has the same arguments as one of the lower pomerons. This is the situation in F^^'^^ 
(modulo permutation of the lower pomerons and/or renaming of the integration variables). 

y(3) y(2,2) 24^ 

Therefore it can already be stated without any calculation that 1/(2'^) and V^^^ lead to the same type 
of integral. 

In the following, we will calculate ^'(3). ^(2,2) can then be obtained by exchanging the coordinates 
and conformal dimension according to (6.24). Inserting (6.14) into the formula for ^ (6.18), we get: 
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Plc||P3c| \pbd\ J \\p2d\\P3d\ \Pbc\ 

=:-y 1^4— ^(4zyJ + l)2 X] '^{-n,-i^c,-i'a,i^d;x,x*). 

' a<->6<-»c 

The step from the first to the second expression is not trivial. It is legal only because the vertex 
without the function is invariant under permutations of the pomerons a, b and c, as can be seen 
from (6.15). (In fact it is symmetric under all permutations of conformal fields if one takes care of 
the different sign of i^d-) It is this part of the vertex which we have written into the sum over the 
permutations, even though the symmetry is not obvious any more. The sum over the permutations 
also extends to x in the sense that it is transformed into a different anharmonic ratio by a permutation 
of the coordinates. In fact, all six anharmonic ratios occur in the sum, some with a minus sign. 
T is defined as follows: 
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We have defined it so that the conformal dimensions of the two pomerons with the same coordi- 
nate arguments come first. It is symmetric in these first two arguments provided the corresponding 
coordinates are exchanged as well, which leads to a different anharmonic ratio, Xjx instead of x. 

T{~i^b, -J^c, -J^a, i^d; X, X*) = T{-i^c, -n, -i^a, ^d] - , ^) (6.27) 

X X* 

The conformally invariant function ^'(2,2) corresponding to V^^'^^ can be obtained from ^"(3) by way 
of (6.24). The only slight difficulty is the anharmonic ratio x. Exchanging the coordinates pd pb 
also changes the anharmonic ratio on which T still depends: 

^^PM_^PM^PadPbc^^_^^ according to (5.17). 

PacPbd PacPdb PacPbd 

So we get for ^'(2,2)^ 

40 2 

*(2,2) = y Y6^(2;0^^^''^^^^' ^ T(-zvd,-i.e,-i^a,J^6;l-a;,l-ar*). (6.28) 

6.3 Summary of the 1 — > 3 pomeron vertex from the irre- 
ducible 2^6 gluon vertex 

This section will summarise the results of the calculations of the previous section. We have seen that 
the 1^3 pomeron vertex has the form of a conformal four-point function: 

14^3p({Pi}; {'^i}) = x*;{'^i}) n |pi,f + ^ , (6.29) 

i<j 

where x = (pabPcd)/ (pacPbd), h = -i/j for i = a,b,c and I'd = i^d- 

The freedom which remains in the otherwise fixed form of the vertex, the function ^, can in this 
case be written as a sum of five terms: 

^{x, X*; {ui}) = ^bfkl{x, X*; {vi}) + {x, x*; {vi}) 

(6.30) 

+ *(2,i)(a;, x*;{vi}) + *(2,2)(a;, x*; {vi}) + *(3)(a;, x*; {vi}) . 



The five terms of the function ^' have the following form: 

40 ^ 2 
y 16 (27r)3 

40 2 
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^BFKL{x,X*;{Ui)) = — J^^J^^Z^^i^d) (4z^d + 1)^ 'B.{-Ua,-V>b,-l^c,yd;X,X*) 

40 5^ , 2 
1 

*(2,1)(X,X*; {^,}) = ^ f!4 (4:.^ + 1)2 \^\-k + iK-^^^+'^r,+^c-u.) |1 _ ^|i + |iK-2..+..+2..) . 
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The sums over T in *(2,2) and ^"(3) run over all six simultaneous permutations of the conformal 
dimensions Ua/b/c and the corresponding coordinates Pa/b/c- The latter entails that x is replaced by a 
different cross ratio according to the permutation of the p s. The two integral expressions S and T are 
defined as follows: 
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6.4 Conclusions from the existence of the vertex 

In the previous sections, we have derived the explicit form of a 1 — > 3 pomeron vertex. It was 

obtained from a 2 ^ 6 rcggciscd gluon vertex from the integral equation for the six-gluon amplitude 
in Extended Generalised Leading Logarithmic corrections to BFKL. This gluon vertex, and hence the 
1^3 pomeron vertex obtained by projection, is local in rapidity. 

This has some; implications. Leading Logarithmic Order BFKL is known to be equivalent to leading- 
order calculations in the dipole approach proposed by Mueller [Mue94, Mue95] . In the dipole approach, 
scattering colour dipoles are assumed to split up repeatedly into smaller ones, which then interact by 
exchanging gluons. This splitting up of the dipoles was fomid to be equivalent mathematically to the 
gluon ladders the BFKL equation was originally derived from. 

How far this correspondence extends to higher orders is not yet known. The state of the art 
regarding this question is rather complex. Peschanski has derived a general \ ^ n pomeron vertex 
from Mueller's dipole approach [Pes97]. However, there are two qualifications to this statement: For 
one thing, the higher vertices were obtained by an extrapolation of the 1^2 pomeron vertex. The 
other point of caution concerns the expression for the 1^2 pomeron vertex which forms the basis 
for the interpolation. Though Peschanski starts out from the dipole formalism, the extrapolated 
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expression is very close to the form of the 1^2 pomeron vertex obtained by Lotter [Lot96] in the 
gluon ladder approach. Hence Peschanski's results constitute no proof that these vertices exist in the 
dipole approach. Braun and Vacca [BV99] have since argued, starting from the generating functional 
of the n-fold dipole densities, that the vertices of 1 — > n pomerons for n > 2 do not exist in the dipole 
approach. 

In contrast to that, the results of the previous sections show that the 1 — > 3 pomeron vortex does 
exist in the gluon ladder approach. This constitutes a discrepancy between the dipole and the gluon 
ladder approach and implies that the correspondence of the two approaches does not extend beyond 
Leading Logarithmic Order. 

6.5 The 1^3 pomeron vertices from reducible 2^6 gluon 
transitions 

6.5.1 Reducible 2 — > 6 gluon transitions 

As was already mentioned at the beginning of Section 6.1, the 2^6 reggeised gluon vertex treated 
in the first part of this chapter comes from the integral equation (6.1) for a six- gluon amplitude which 
occurs in unitarity corrections. Besides V2^6, this equation contains other components, some of which 
are also 2^6 gluon transitions local in rapidity. Unlike ¥2^6, however, these terms do not represent 
irreducible 2 — > 6 gluon vertices, and we will therefore not call them "vertices" but "transitions" . They 
can be expressed in terms of the function {VD2) of which the 2^4 gluon vertex is composed. They 
are therefore related to the irreducible part of the four-gluon amplitude. 

We will not consider the completely reggeising part of the six-gluon amplitude here. It is expected 
to yield a 1 ^ 3 pomeron vertex, just as the reggeising part of the four-gluon amplitude gives rise to 
a 1 ^ 2 pomeron vertex. The latter vertex was discussed in [BRV02] and found to be subleading in 
the large limit. 

We will derive the 1^3 pomeron vertices contained in the partly reggeising components in the 
second, third and fourth line of (6.1): 

• • • + ^ ] faia^as fa^a^ae -^(1) 2, 3; 4, 5, 6) 

+ Yl 'P'"'""'' Sa,ae (1, 2, 3, 4; 5, 6) (6.33) 

+ 'i"^"^"-'"* Sa,ae J{1, 2, 3, 4; 5, 6) + . . . 

The sums run over simultaneous permutations of colour indices and momenta/coordinates. However, 
the symmetry of both colour tensors and spatial functions is limited. These terms arc not totally 
symmetric in the indices, as ^2^6 was. This forces us to make a choice for the symmetries of the 
three-pomeron state we use for the projection. If, for instance, one of the pomerons is distinguishable 
from the other two, the result of the projection may depend on which pairs of arguments of the 2^6 
gluon transition term it couples to. Which pomerons are distinguishable depends on where there are 
t-channel cuts, which in turn depends on the physical process considered. The existence of cuts also 
places restrictions on which arguments of the transition terms each pomeron can couple to. Since we 
have no particular process in mind, we will for the larger part of this section choose a superposition of 
three-pomeron states in which the pomerons' six arguments have all possible permutations: 

i. J2 ^('^<"°)*(7ri,7r2)£(''^'°)*(7r3,7r4)^(''-°^*(7r5,7r6) 

■ 7ren(l,2,3,4,5,6) 

This is a completely symmetric state of six reggeised gluons, with the additional requirement that 
they form three pomerons. It implies that we assume the pomerons to be indistinguishable. This 
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(6.36) 



has the advantage that in our projection we will encounter all the terms which can possibly result 
from a projection onto an arbitrary three-pomeron state, provided we refrain from making use of the 
symmetry to find cancellations. Later, in Section 6.5.6, we will also consider the case of distinguishable 
pomerons. 

Two of the terms in (6.33) contain the colour tensor (^"ia203a4^ jg defined as follows: 

Here are the generators of SU{Nc), half the Gell-Mann matrices for = 3. d^'^o-^o-aai obeys the 
following relation: 

^ 1 2 3 4 _ '^0102 '^0304 + ^ {dkaia2^ka3a4 ./'iaia2 /ia3a4) • (6.35) 

fabc are the structure constants and date the symmetric structure constants of the algebra. 
The three spatial functions have the following form: 

L(l, 2, 3; 4, 5,6) = ^ [{VD2){12, 3; 45, 6) - {VD2){12, 3; 46, 5) + {VD2)il2, 3; 4, 56) 

- {VD2){U, 2; 45, 6) + (yr'2)(13, 2; 46, 5) - {VD2){13, 2; 4, 56) 
+ {VD2){1, 23; 45, 6) - {VD2){1, 23; 46, 5) + {VD2){1, 23; 4, 56)] 

7(1, 2, 3, 4; 5, 6) = g^[-{VD2){l, 234; 5, 6) - (yr'2)(123, 4; 5, 6) + {VD2){U, 23; 5, 6)] 

J(l, 2, 3, 4; 5, 6) = g^[-{VD2)il3A, 2; 5, 6) - {VD2){124, 3; 5, 6) 

+ (F£>2)(12, 34; 5, 6) + (Fr>2)(13, 24; 5, 6)] 
The function {VD2) can be expressed in terms of the function G. 

{VD2){1, 2; 3, 4) = ^ [G(l, 23, 4) + G(2, 13, 4) + G(l, 24, 3) + G(2, 14, 3) (6.37) 
- G(12, 3, 4) - G(12, 4, 3) - G(l, 2, 34) - G(2, 1, 34) + G(12, 34)] 

This formula is in fact analogous to (6.3), but since there are only nine non-zero terms, the sum 
has been written out. {VD2) is symmetric under exchange of the arguments 1 and 2, and 3 and 4, 

respectively, and under exchange of its two pairs of arguments. This can be easily derived from the 
symmetry of the G function under exchange of its first and third argument. 

{VD2){1, 2; 3, 4) = iVD2){2, 1; 3, 4) = (m)(l, 2; 4, 3) 

(6.38) 

(yi?2)(l,2;3,4) = (T/i?2)(3,4;l,2) 

We will shortly derive the 1^3 pomeron vertices contained in the three reducible 2^6 gluon 
transitions. We will see that they are composed of the same terms as the pomeron vertex obtained 
from V2^6- To determine which vertex will dominate in the large Nc limit, wc will calculate the colour 
constants for general N^- Before that, however, we will rewrite the irreducible vertex computed in 
Section 6.2 for general Nc and in a concise graphical notation for the convolutions of the spatial part. 

6.5.2 The irreducible vertex for general Nf. 

The colour constant of the irreducible vertex is datcdabc- For Nc = 3 this is equal to 40/3. For general 
Nc, it can be derived from the following relation: 

N^ - 4 

dabcdabd = Scd 

dabcdabc = ^{Nl-A){Nl-l). 
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Collecting all the prefactors, the vertex takes the following form: 
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(6.39) 



The prefactors j of the a terms come from a factor 4 for the four permutations (6.4), from the factor 
I in (6.3) and the factor i in (5.45). The prefactor of the BFKL term differs because Kbfkl already 
contains a factor Ncg^, lacks a factor 2 relative to a, and because there is only one BFKL term for six 
of the a terms. 

The functions represented graphically correspond to the integrals we obtained in Section 6.2. They 
depend on the conformal dimensions of the four pomerons and on the cross ratio of their coordinates. 
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Note that the graphical symbols for the alpha terms, as defined here, are not symmetric with respect 
to exchange of the pomeron indices. This does not constitute a problem here, since we have decided to 
use a three-pomoron state which is a superposition of the six possible permutations. In Section 6.5.6, 
where we consider distinguishable pomerons, we will write indices next to the pomeron lines. In 
the immediately following sections, however, the identities in the graphical notation are always to be 
understood as modulo permutations of the indices of the three lower pomerons a, b and c. We will 
take care always to average over these permutations. 
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6.5.3 The colour-antisymmetric L terms 

The first set of reducible terms we will now deal with are the ones containing the spatial function L. 
They have a colour structure given by the antisymmetric structure constants. 

^ ] faia2a3 f ataman L{\, 2, 3; 4, 5, 6) = /oia2a3 fata^ae -^(1) 2, 3; 4, 5, 6) 

(6.41) 

"I" faia^ai fasa^ae 2, 4; 3, 5, 6) + . . . 

The sum runs over the same permutations as for the irreducible vertex, that is the ten possibilities 
of dividing a set of six indices into two groups of three. The indices of the / tensor are always in 
ascending order. 

Unlike V2->6; this expression is not completely symmetric in the reggeised gluon indices. The / 
tensors are symmetric only under cyclic permutations and antisymmetric under anti-cyclic ones. The 
spatial function L is even less symmetric. It can be worked out from its form (6.36) in terms of the 
function G that it is invariant under the exchange of its first group of three arguments with its second 
group of three arguments and of the first and third argument in each group. 

Because of the complete symmetry of the three-pomeron state we project onto, it is sufficient 
to perform the projection of one of the terms in the sum (6.41). Therefore the sum over the ten 
permutations in (6.41) can then be discarded and replaced by a factor ten. To find out which different 
terms result from the projection, we have a look at the permutations of the six arguments of three 
pomerons. There are 6! = 720 permutations of six arguments. A factor 2^ = 8 can be divided out for 
the exchange of each pomeron's two arguments. Since they are symmetric, no averaging over these 
possibilities is necessary. The remaining 90 possibilities decompose into the 15 possibilities in which 
three identical pomerons can couple to six positions and six permutations of the three pomerons. We 
will deal with them separately. 

Because the colour tensor / vanishes when contracted with a Kronecker delta, each pomeron has 
to have one argument in each of the two groups of three. This is the case for 6 of the 15 possible ways 
of coupling thcm.^ So only six of the fifteen possibilities of the pomerons' coupling contribute. Since 
the function L is not symmetric under exchange of the middle argument in each group with one of the 
others (of the same group), a pomeron coupled to this argument is distinguishable from the others. 
Two different terms arise depending on whether the same pomeron couples to the middle argument in 
each group or not. They can be represented graphically. 




TT 

Since the other two arguments within each group are exchangeable, these are the only two distinguish- 
able terms. Two of the six terms are of the first type, four of the second. 

However, the antisymmetry of the / tensors causes these terms to cancel out in pairs. A certain 
way of coupling three pomerons to two groups of three arguments can be written as a permutation of 
the numbers one to three. The permutation maps 1 to 1, 2 or 3 depending on which argument of the 
second group is attached to the same pomeron as the first argument of the first group, and so on. The 
colour constant corresponding to the permutation (7ri7r27r3) is 

/oi02a3./a,ia,2a,3 = X{t^) Nc{N^ - 1), 

where x(7r) is the character of the permutation, +1 for even and —1 for odd permutations. Exchanging 
the first and third argument in one argument group of L leaves the spatial expression unchanged 
(because of L's symmetry) but reverses the permutation and hence the sign of the colour factor. 

^One pomeron has to be attached to the first index of the first / tensor, or equivalently the first argument of L. Its 
other argument can be attached to any argument in the second group, which gives three possibilities. The pomeron 
attached to the second argument in the first group has two possibilities, the third just one. 
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Therefore, all terms come in pairs which cancel out. The terms with the L function do not give rise 
to a 1 — > 3 pomeron vertex. 

We could have obtained this result more quickly by invoking symmetry arguments. However, our 
aim was to present all possible terms which might result from a projection onto an arbitrary three- 
pomeron state. In a projection onto a three-pomeron state which is not completely symmetric, the 
above cancellations do not occur and non-zero terms remain, see Section 6.5.6. 



6.5.4 The terms with the function / 

The sum of the terms with the function / runs over all fifteen ways of dividing six arguments into one 
group of four and one group of two. The arguments are always in ascending order of the indices: 

V d-i-^-3a4 ^^^^^ 2, 3, 4; 5, 6) = d'^^'^^'^'''' Sa,a, /(I, 2, 3, 4; 5, 6) 

(6.42) 

^^aia.a3a5^_^^_^^ j(1^2,3,5;4,6) + ... 

The terms containing the spatial function / have very little symmetry. Both the colour structure 
and the function I are invariant under exchange of the last two arguments. Besides, the colour tensor 
^01020304 jg invariant under cyclic permutations and under inversion of the order of the indices. The 
function / is symmetric under exchange of its first with its fourth and of its second with its third 
argument, respectively. This implies invariance under reversal of the order of the first four arguments, 
a symmetry shared by the colour tensor. This and the exchange of the last arguments is their only 
common symmetry. 

Since the three-pomeron state we project onto is completely symmetric, we average over the permu- 
tations of the three pomerons attached below the 2^6 gluon transition. We can then again discard 
the sum over the fifteen I terms and replace them with a factor 15. (This cancels a factor from the 
average over the pomerons, so we have in effect replaced the sum over the terms by a sum over the 
ways the pomerons couple.) 



The colour factors 

We will first calculate the possible colour factors which can occur in the terms with the / function. We 
will see that all of them actually occur. To make it easier to match them to the spatial expressions, 
we will write them in a graphical notation. The colour structure of the terms is: 

Half-circles denote Kronecker deltas. We will now calculate the possible colour factors by attaching 
three of them from below in all possible ways. For instance, one pomeron (which has a singlet colour 
structure described by a Kronecker delta) could couple to the first two legs, another one to the third 
and fourth, and the third to the last two. This expression factorises, and we can treat each factor 
separately. 

^S? n =^A? • n 

The second factor is the easier to evaluate. It is the contraction of two Kronecker deltas, ie the trace 
of the unit matrix in the adjoint representation: 

n = SabSab = - I . (6.44) 



The other factor can be computed with the help of relation (6.35). 



^abScd + ^ {dkabdkcd - flab fled) 



^ab^cd 



^ SabSed SabScd = ^ {N^ - if (6.45) 
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Both the three-index d and / tensors vanish when contracted with a Kronecker delta in this way. So 
we obtain the first colour factor: 

fS? n =^(^c-l)'- (6-46) 

www -^J'c 

The second possibility is for a pomeron to couple to the first two legs (again) and for the other 
two to link the rf"'"^'' tensor to the Kronecker delta. Due to the symmetry of the delta symbol, there 
is only one way of doing this: 

w v.;^^^ 

It is easy to see that this "linking" has the same result as attaching a Kronecker delta to the two last 
indices of the (f^'^'^ tensor: 

n = ^cf^de^ef = Scd = (6.47) 

\wy w 

So we get the same result as for a d"^'''^ tensor contracted with two Kronecker deltas (6.45). 

^S? n = = ^ (iVe^ - 1)^ (6.48) 



Because the d'''"^'^ tensor is not completely symmetric in its indices, terms in which the delta symbols 
attached to it cross are different from the ones we had already. We obtain for the rf"'"^'' tensor alone: 



= d"'^'^'^5acSbd = 



SabScd + ^ {dkabdkcd — flabflcd) 



Sachd 



= Saa + ^ {dkabdkab - flab flab) 

= 2^^ ~ + i ~ ~^^~\ ~ = ~ ^ ~ • ^^-^^^ 

As for the uncrossed diagram, there are two different colour factors related to this one. In the first, 
the Kronecker delta from the transition term is contracted with a pomeron separately from the 
tensor; in the second, they are linked. The factorising diagram gives the result: 

="^^^'"'^'- ^^-^^^ 

Because of (6.47), the non- factorising diagram gives the same result as the contraction of the d"'^'''^ 
tensor alone. 

There arc no further difFcrcnt colour factors. All other diagrams are equal to one of those already 
presented by way of the symmetry of the colour tensor. For instance, 

(6.52) 

yvjy w w 

because of the tensor's invariance under cyclic permutations. 
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The spatial part 



The spatial part of this transition term is given by the function I: 



/(1, 2, 3, 4; 5,6) = g^[- (F£)2)(l, 234; 5, 6) - (yr'2)(123, 4; 5, 6) + {VD2){U, 23; 5, 6)] . (6.53) 



We will now look at all possible ways of convoluting three pomeron wave functions with it. 

Let us first have a look at the terms in which one pomeron couples to the last two arguments of I. 
There are three possibilities of coupling the remaining two pomerons to four arguments. They can 
be represented graphically. A half-circle denotes a pomeron amplitude which is convoluted with the 
respective arguments of I. We get three diagrams: 



The legs of the / function have been marked to indicate its symmetries: Exchanging the pair of 
arguments marked in the same way, or the unmarked ones, leaves it invariant. However, the colour 
factor changes if one exchanges only one of the two pairs of marked legs, since d"'"^'' is invariant under 
reversal of its indices but not under arbitrary permutations. So the latter two diagrams represent the 
same spatial convolutions but have different colour factors. 

Let us deal next with the terms in which the last two arguments of / are connected by pomerons 
to two arguments which also form a pair (ie under the exchange of which I also is invariant). There 
are two possibilities here: The other pair can either be the second and third argument, or the first and 
fourth argument: 



These terms have the same colour factor (6.52), taking into account (6.47). From each of them a new 
diagram can be obtained by exchanging the two rightmost legs. These diagrams are equal in both 
colour structure and spatial part since this exchange is a common symmetry of both the colour tensor 
and /. So they can be accounted for by multiplying the displayed terms by two. 

There remain eight terms. They can be obtained from each other by applying the exchanges of 
arguments with respect to which the function / is invariant. Therefore their spatial convolutions are 
all equal. However, two different colour factors occur, with four terms having the same factor. The 
reason for this is that d"''"^'^ is invariant only under reversal of its indices, which amounts to exchanging 
the first with the fourth and the second with the third index simultaneously. Graphs which differ by 
only one of these exchanges have different colour factors. Here is one representative of each group, the 
spatial part of which is equal: 



We will now put the spatial diagrams and the colour factors together and sum up all the terms. In 

addition, we will write down the average over the permutations of the pomerons, the only symmetry 
of the three-pomeron state which we have not yet taken into account. The result is the following sum 
of the five different terms with the given prefactors: 






a<-»o<->c L y w www 
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- E 



6 V 



w www 



(6.54) 



What now remains to be done is to express / in terms of {VD2), {VD2) in terms of G, and G 
in terms of BFKL kernels and alpha functions. We will see that a G function convoluted with three 
pomerons always leads to one or the other. Using (6.53) and bearing in mind the symmetries (6.38) 
of V and the fact that a pomeron wave function with two equal coordinate arguments vanishes, we 
obtain: 







www I I I I 




(6.55) 



w 

The sum of all terms including the correct prefactors now reads: 



5 E 

a<-^b<-*^c 



2Nr. 




{Nl - 2){N! - 1) CO - 4 (2iV2 - 3) CO + 4 {N^ - 2) 




"1^ 



(6.56) 



Now the three different terms with V have to be expressed in terms of G and then as a or BFKL 
terms. Before we do that, we will show that a G function contracted with pomeron wave functions 
always is a sum of the two, ie that the functions b and t cannot lead to anything other than the 
well-known BFKL term. Leaving aside the functions a and s which lead to a, G is composed of three 
functions: 

G(l, 2, 3) = y [2 c(123) - 2 6(12, 3) - 2 6(23, 1) + 1(12, 3) + 1(23, 1) + . . . ] . 

The function c does not play a role when one projects onto pomerons, since it always identifies both 
coordinates of the pomerons. The functions b and t may also identify two coordinates of a pomeron 
and thus vanish under projection. But if, for instance, no pomeron couples to the pair of arguments 2 
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and 3, one of the b and t functions remains. In the simpler case with only two pomerons, this looks 
like this: 

G(l, 2, 3) (g> 2) 3) = 

2 2 

= |- (-2 6(23, 1) + 1{23, 1)) ® 2) £;("'"0)*(1, 3) + y a(2, 1, 3) ® 2) £;(''''>o)*(l, 3) . 

Since the pomeron wave functions are symmetric, only the symmetric part of b and t survives the 

projection. Therefore we are allowed to exchange their arguments. Wc symmetrise the expression in 
these arguments and add a c function (whose projection is zero anyway) to obtain a complete BFKL 
kernel: 

^ (-2 6(23, 1) + t{23, 1)) O £(''-o)*(l, 2) ^(''''>o)*(l, 3) 

2 

= |- (-2 6(2,l) + i(2,l)) ®£;("»'0)*(l,2)E('''"0)*(l,2) 

= (^(12) - 2 6(1, 2) - 2 6(2, 1) + t(l, 2) + t(2, 1)) (5§ £;('^-o)*(l, 2) .E(''''>o)*(1, 2) 
= i ^ [Ai A2(i^sFKL $2) (1, 2)] ® i;(^-o)*(l, 2) i;('^^>o)*(l, 2) . (6.57) 

Note that the prefactor of the BFKL kernel is half that obtained when the middle argument of G is 
left empty, sec (5.46). The Laplacians only serve to cancel propagators defined into the BFKL kernel. 
The wave function $2 is the amplitude attached to the kernel from above, which will in our case be 
replaced by an as in Equation 6.7. 

Now we can evaluate the three terms containing the 2 — > 4 vertex function V. Taking into account 
the symmetry of G, one arrives at the result: 




(6.58) 



In the second V term, the first and the last G term both contain a BFKL term. They partly cancel 
each other. 

We now plug the formulas (6.58) into (6.56) to compute our final result.^ 



^In principle it would be sufficient now to average over the cyclic permutations of the pomerons since due to the 
symmetry of a only one pomeron can be distinguished from the others. But we will leave the sum over all permutation 
to make the comparison with (6.39) easier. 
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One can sec that this vertex dominates over the one derived from the irreducible 2^6 gluon vertex 
in the large A^c limit: Here the first term is of order N^, while (6.39) is of order A^^ . This differs from 
the situation in the case of the 1 — > 2 pomeron vertex, where the irreducible vertex was A^c-leading 
(see [BRV02]). 

6.5.5 The terms with the function J 

The terms with the J function are very similar to the I terms treated in the previous section. The 

sum 

^ d---^-^-- Sa,ae J{^, 2, 3, 4; 5, 6) = d-=^-i-3a4 ^^^^^ 2, 3, 4; 5, 6) 

(6.60) 

^^a^aiaaaa^^^^^ j(l 2,3,5;4,6) + ... 

runs over the same fifteen terms obtained by dividing six arguments into two groups of four and two. 
The only difference is in the colour structure: While the indices of the d"'^"'^ tensor were in ascending 
order in the case of the / function, here the first two indices are interchanged. 

The first steps of the derivation of the J terms are exactly analogous to the terms with /. The same 
colour factors occur. The functions / and J have the same symmetry, so the graphs which led to equal 
terms with the I function are also equal for the J function. The colour tensor is again symmetric under 
simultaneous exchange of the arguments/colour indices 1 with 4, and 2 with 3. Due to the difference in 
the colour tensor, this is brought about not by its invariance under reversal of the order of the indices 
(as for the I terms) but by invariance under cyclic permutations (by two positions, in this case). 

What is different, however, is the matching of colour factors to spatial convolutions. One can see 
this graphically by drawing a tensor d''"^'^'^ as a d"''^'^ with crossed legs. In the previous section we 
derived two different contractions of the d"'^"'^ alone and found a third to be equal to one of them. 
One of these colour factors remains the same for the "crossed" d"''"^'^ tensor, but the other two are 
interchanged: 
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The additional Kronecker delta affects the colour factors in the usual way: If it is linked to the d'**"^'' 
tensor (as in (6.48)), the result is the same as for the rf"^'^'' tensor alone, otherwise there is an additional 
factor (A^2 _ i) 

We will give the colour factors for each term below, together with its decomposition in terms of 
the 2 — > 4 vertex function V. In contrast to (6.55) we give all terms whose spatial part is equal but 

which might differ in the colour factor separately. As it happens, the two pairs of equal terms also 
have equal colour structure. The multiplicity of the seven terms is given in the first column. 



W W 



2g'Cv 





V 




V 



_(iV2-l)2 



www www ^iyr. 



w www 



_(iV2_ 1)3 





2N, 



{Nl - 1) 



2Nr^ " ' 



4- 



Adding them all up, we obtain: 



1 ^ 



2Af. 



WW 2N^ ' > 



Putting in (6.58) we get the result: 

1 




2 {Nl - 2){Nl J - 4 {2Nl - 3) V J + 8 (7V^^ - 2) ( V 





12 



(2iV4 - 7) - i (7V,2 - 2) ^ i (2A^,2 - 3)(^T^ 




+ ^(2^c^-3)(j:«:])-^(2Ar,2-3)^(|?^ 



(6.61) 
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As was the case for the / terms, the first term is of order and therefore dominates in the large 
Nc limit. The other terms of the reducible vertices and all of the irreducible one (6.39) were only of 
order N^. 

6.5.6 Distinguishable pomerons 

Having computed the projection of reducible 2 — > 6 gluon transitions onto a completely symmetric 
tliree-pomeron state, wc will now have a look at the opposite extreme: Three distinguishable pomerons 
which couple to specific arguments of the transition terms. We will choose as an example the three- 
pomeron state in which pomeron a couples to arguments 1 and 2, pomeron h to arguments 3 and 4, 
and pomeron c to arguments 5 and 6. This gives the following thrcc-pomeron state: 

£(^-°)*(l, 2) £;('^'-°)*(3, 4) i;(''-°)*(5, 6) 

The calculation of the projection is analogous to the projection onto a symmetric thrcc-pomcron 
state performed above. The one additional difficulty is that the terms are not only distinguished by 
whether a pomeron couples to a given pair of arguments, but also which pomeron couples to them. 
The pomerons in the graphical notation have to be labelled a, h or c. The identities (6.55) and (6.58) 
for the / terms and analogous ones for L and J have to be rewritten to take the distinctiveness of the 
pomerons into account. Other than that, the derivation is quite analogous to the previous sections. 
Wc will not present it in detail here, but give just the results. In all diagrams, two of the pomerons are 
interchangeable, either because they couple to the same pair of arguments or because of the symmetry 
of a. Where necessary, the third pomeron is labelled in the diagrams. 

The result for the L terms is: 

5a,a, Sa,a, S^.a, f-^-^as fa,a,a, L{1, 2, 3; 4, 5, 6)) £(''»'0)*(1, 2) £(''^'0)*(3, 4) £;K.0)*(5, 6) = 




a b c 



One can see that the coefficients of some terms in which the pomeron with index h is special differ 
from those with indices a and c. The reason for this is that the pomeron b couples to the arguments 
of L which are distinguishable from the others, see Section 6.5.3. 
The projection of the terms with the functions I and J is: 

<5a,a. <5a3a. ^a.aa (E d^^'^^'^^'^' 6a,ae 2, 3, 4; 5, 6)) ® i?(''-°)*(l, 2) i?(''-0)*(3, 4) £;(--0)*(5, 6) = 
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The respective iVc-leading term is symmetric in the pomerons. The other a terms are stih symmetric 
in the pomerons a and c. This was not obvious a priori. It is a consequence of the symmetries of I 
and J and of the fact that their two groups of arguments arc ordered with ascending indices. 

All expressions are of the same order in Nc as in the result of the projection onto the completely 
symmetric three-pomeron state. This applies not only to the sums of terms with L, I and J, but also 
to each type of diagram separately, even though the exact prefactors have changed. In particular, the 
same a diagram of the / and J terms dominates in the large Ng. limit. 
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A.l Obtaining the graph classification from the correlator of 
path integrals 

This section closes the gap left between Sections 2.2.3 and 2.2.4. In Section 2.2.4 we switched to the 
Feynman graph formulation of our problem and explained the graph classification and the different 
prefactors on its basis. It is interesting to see why this is allowed and how the same factors follow from 
the continuation of our derivation on the basis of the path integral formalism described in Section 2.2.2. 

Wc shall now start out from the expression for the S'-matrix and the expansion of the gluon potential 
path integrals V". See Section 2.2.2, Equation 2.9 for a definition of the V". 

S{b,{xt}) = ^(e„^, (V0_, (V?)^^, (V?)^^, e^,0,ye,,^ (V^)^^, (V^)^^, (Vf)^^, ^.v) (A.l) 

(^tU = - igKi^ii - \ a^Ki^U^'^^'U - ^ a^KKKi^-^'^'^^'^'U + o{g^) (A.2) 

Bl.r'^ = I^Jz^ B,{z) (A.3) 

The expansion of the S'-matrix in g (or, cquivalcntly, in the J up to 0{g^) contains a huge number 
of terms. Wc have to extract those which represent perturbative odderon exchange. Since they are of 
the order g^, they all contain a correlator of six ^. 

In a perturbative calculation this correlator can be expanded in correlators of two B^ j . The 
reason for this is that even though the B^ ^ are expansion coefficients of path integrals over the gauge 
potentials, they will ultimately lead to correlators of gluon potentials themselves. (This is shown in 
Section A.2 for the two-correlator.) These correlators in turn can be expanded in gluon propagators 
in perturbative QCD. 

We do not want terms in the expansion which contain correlators of two B^ ^ from the same baryon. 

They represent a gluon exchange between two quarks of the same baryon, ic a self-energy correction. 
We are interested only in the lowest-order odderon exchange in which three gluons are exchanged 
between two baryons represented by three quarks. Therefore we can restrict ourselves to those terms 
with three B^ ^ from each baryon: The expansion of a six-correlator in which a different number of 

B"^ ^ comes from each hadron always contains self-energy terms. We obtain six terms of the expansion 
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which are relevant for triple-gluon exchange: 

\^a,l -Oa',1 -Oa",l ^h,2 ^b" ,2/ 



\-Oa,l ^b,2) \^a',l ^b',2l \^a",l ^b" ,2/ 
+ {Ba,l -^6,2) By', 2) i^a" .1 By ,2) 

+ {Bi, Bii,) {Bi,^, Bi^) {bs,^ Bi:.2) 

+ (-^a,! By,.^ {Ba'.l -B^,2> {Ba" ,1 By ,2) 

+ (4,1 4<2) (4'm 4< 2) (4'M 4,2) 

+ self-energy terms. (A-4) 



The procedure now is to list every six-correlator containing three ^ from each baryon and expand 
them according to (A. 4). We will see that in many cases there appear equivalent terms which can be 
summed up by calculating one of them and multiplying them with a constant factor. The different 
terms correspond one-to-one to the different Feynman graphs. 

The fact that some terms in the expansion of the S'-matrix arc equivalent is one source of constant 
prefactors. Other factors arise from the coefficients in the expansion of V" (A. 2) and from the colour 
structure. The colour tensors and the corresponding factors have been treated in Section 2.2.4. They 
arc listed again in Table A.l. The factors from the exponential series of Vf , which arc straightforward, 
and the factors from equivalent terms, which are quite tricky, will be discussed in the following. 

For coping with the multitude of terms, we again employ the classification of the odderon-proton 
coupling described in Section 2.2.4. In the language of the above formulas, the three types of couplings 
are distinguished by how many B^ ^ come from each of the three V" of a baryon: 

a) The third-order term of one of the V" and the Kronecker delta of the others contributes. 

b) One second-order, one linear and one zeroth-order term come from the three expanded V°. 

c) Each of the V° contributes the linear term from its expansion. 

One can see immediately from Equation A. 2 that a type (a) coupling entails a factor of ^, (b) a 
factor of ^ and (c) a factor of one. This gives the factors from the expansion of V" listed in Table A.l. 

The derivation of the factors due to terms which represent the same Feynman graph and are 
numerically equal is more involved. Practically the only way to obtain them is to write down at least 
a couple of terms for each graph type and their expansion according to (A. 4). This will be done in the 
following for two examples, (aa) and (cc). 

Two (a)-type couplings yield nine different six-correlators, depending on which quark the three 
gluons couple to (or, to put it differently, which of the V" contributes its third-order term). Here are 
some of them: 

(aa): (BJ B{^^ B\ -^ B^ ^ ^1,2 -^1,2) + (4,i 4,i 4,1 -^1,2 ^1,2 -^^,2) 

+ • • • + (-83,1 -^3,1 Bss 4,2 4,2 4,2) • (A.5) 

They represent a three-gluon exchange between the first quark of the first baryon and the first quark 
of the second baryon, between the second quark of the first baryon and the first quark of the second 
baryon, and so on to the third quark of the first baryon and the third quark of the second baryon. 
Let us have a closer look at the first of these six-correlators. It can be expanded in two-correlators 
according to (A. 4). 

(aa), term: {B^, b(, b(, B^, Bt[, bQ = {Bl, bI,) {b(, b(,) {b(, bQ 

+ (-^1,1 -^1,2) {Bi^i 4,2) (-^1,1 4,2) 

+ (4,i4,2)(^G^'2>(4;'i4,2) 
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+ {bi^bC){b(,b(,){b(,bI,) 

+ self-energy terms 

+ (J^''" (^'^''^ (5'="''' + S"'^" S"'"^' S""'^] 

+ self-energy terms (A. 6) 

Here we have used the fact that the colour part of the gluon propagator is a Kronecker delta. The 
configuration space part of the propagator has been wrapped up in Xab which will be calculated in the 
following section and is defined as follows: 

(^a,l4^,2) = -^'''^a6. (A.7) 

It should now be obvious that all six terms in Equation A. 6 are equal: When the symmetric colour 

tensors of the oddcron are contracted with the delta symbols, it does not play a role which index 
of d'^" " is contracted with which index of d'^'^ . Therefore the first term for graph type (aa) is equal 
to -6 Xfi. 

Lot us now have a look at the graph type (cc). It contains only one term from the expansion of 
the S'-matrix since there is only one way to couple three gluons to three different quarks: one to each. 
Its expansion in two-correlators is the following: 

(cc): ^2,1 B^^i Bf 2 -62,2 -^3,2) = 

+ 

+ 

+ 

+ 

+ 
+ 

This time the six terms are all different: They differ not only in the colour indices but also in which pairs 
of quarks exchange a gluon. They correspond one-to-one to the six Feynman graphs one obtains when 
assuming that the three quarks of a baryon are distinguishable. This preposition is quite necessary 
since the quarks differ in their position according to our geometric model for the proton. 

The graph types (aa) and (cc) illustrate two extreme cases: All terms may be equal because of the 
symmetric colour structure of the oddcron, which leads to a factor 6 for the expression obtained in the 
case (aa), or all may be different, which leads to a factor 1 for (cc). Graph types (ab) and (ac) also fit 
into this simple scheme, resulting in six equal terms and therefore a factor of 6 like (aa). For the type 
(be), the expansion (A. 4) contains three pairs of equal terms, giving a factor 2 for each different term. 

Two (b)-type couplings are the most interesting case. The six terms of the expansion in two- 
correlators split into a set of two and a set of four equivalent terms. They correspond to the graph 
types (bbl) and (bb2) whose prefactors differ by a factor of 2. The two equivalent terms are the 
ones with two X functions with the same set of indices. They get a prefactor 2. They correspond to 
Feynman graphs in which the pair of gluons coupling to the same quark is the same at both ends, ie 
type (bbl). The four equivalent terms (bb2) contain X functions whose pairs of indices are all different. 
They get a factor 4. 



/ f>d \ I f}C jjd \ I rye f}d \ 

K^l.l ^1,2) \^2.l ^2,2) {^3.1 ^3,2) 

{Bl.i B12) {B2.1 Ma) (-^3,1 -^2,2) 

K^l.l ^2,21 \^2,l ^1,21 \^3,1 ^3,2/ 

/ f}d' \ I f}c' f}d" \ I f}c" f}d \ 

^2.2) \^2,1 ^3,2) \^3,1 ^1,2/ 

/ r)d" \ / r>c' f>d \ I f>c" r>d' \ 

-C*3,2/ {^2,1 ^1,2/ \^3,1 ^2,2) 

I Rc tid" \ I t>c' Dd' \ I fic" fyd \ 
^3,2/ \^2,1 ^2,2) (-03,1 -^1,2/ 

self-energy terms. (A. 8) 



A. 2. Integrating out the light cone coordinates 
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The relative prefactors of the different graph types have now been adequately explained on the basis 
of our derivation of high-energy scattering in position space. As Table A.l shows, the prefactors are the 
same as in the derivation from the Feynman graph formulation but differ in their composition (confer 
Table 2.1). What remains to be done is to actually calculate the two-correlators (A. 7). This will be 
done in the following section. 



Type 


Expansion of Vf 


Colour 


Equivalence 


Ctype 


(aa) 

(ab) 

(ac) 

(bbl) 

(bb2) 

(be) 

(cc) 


1 

36 
1 

12 

1 
6 
1 
4 

1 
4 

1 
2 

1 


1 

4 

1 

8 

1 
4 

1 

16 
1 

16 
1 
8 
1 
4 


6 
6 
6 
2 
4 
2 
1 


1 (by definition) 

3 
2 

6 

3 
4 

3 
2 

-3 
6 



Table A.l: Prefactors of different graph types originating from the expansion of V°, the colour structure 
and the summing up of equivalent terms. 



A. 2 Integrating out the light cone coordinates 

In the previous section, wc obtained terms containing correlators of two ^ , the expansion coefficients 
of path integrals over the gauge potentials. We now want to relate them to an expression in perturbative 
QCD which we know, such as the gluon propagator. To that end, we first express the correlator in the 
expansion coefficients of the potentials themselves. 

(^^,1 4'2) = U W j^^ dy'' Bi{y)\ , (A.9) 



where 



Bl^iT^ = / dz'^ B^(z) and BI{z)t^ = B^(z) . (A.IO) 
The paths run along the light cone, in opposite directions for each scattering baryon. We obtain 

(^a,i4'2) = Q' dx+Bl{x'i,x+,0) J dy-Bt{xlO,y-)'j 
= J dx+ J dy-{Bl{x1,x+,0)Bi{xlO,y-)) 

= J dx+ J dx- J dy+ J dy-{Bl{xl,x+,x-)Bi{xly+,y-))6{x-)6{y+). (A.ll) 

Xi and X2 are the transverse coordinates of the paths F" and F2, respectively. Those paths are the 
trajectories of a quark in each baryon. See Figure 2.3 for a graphical representation of them. The 
formulation in the last line allows us to introduce into our calculations the gluon Green's function 

D;i{x -y,x+-y+,x--y-) = {-if (B^(a;, x+, x") B^iy, y+, y-)) . 
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Appendix A. From correlation functions to Feynman graphs 



It already contains the factors —i which in our derivation come from the expansion of V^. D is 
translationally invariant, ie it depends only on the relative coordinates of the gluon potentials. We 
will use the Fourier transform D'^{k, of D: 

f d^kdk+dk- 



= -Jdx-Jdx-Jdy-Jdy-J 2^2^^, 

• k+, k-) ei(-fe«-»=^)+ife+(^--!/-)+l'=-(^+-?/+)) 5{x-) 6{y+) 

= - jdx+ j dy- J f^^!^^Df_{k,k+,k-) ei(-fe(-?-4)-i'=+y-+ifc-x+) 



/ 



= -2 



2 (27r)4 '^"^ e-^'°(''?-'^^) 27r-2 5(fc+) 27r-2 5(fc-) 

y"^^D^l(fe,0,0)e-''=('=?-'=^). (A.12) 



The additional factor 2 in the denominator of the Fourier integral reflects the integral measure of the 
light-cone coordinates, dx~^ dx~ = 2 dx^ dx^. So we have arrived at the gluon propagator in transverse 
space: 

-2 / ^^;l(fc,0,0) e-^'=(<-==^)=-2 5'^'^ /■ ^ ^1^+^ e-ifeCa'?-!) 

= / 4^ e-*'=(-?-»'^) =: -,5'='' . (A.13) 

X is thereby defined as minus the spatial part of the gluon propagator in transverse space. The minus 
sign is conventional and represents a factor {—i)^ from the expansion (A. 2). A non-zero gluon mass was 
introduced to make the two-dimensional Fourier integral convergent. In gauge-invariant expressions 
which we will obtain later it can be set to without causing a divergence. All that remains to be done 
now is to perform the Fourier transformation: 



OO 271 





27r 

dip 6*'="°^*' 



OO 



Here n=\k\ and r = |a;" — ajj], the distance bridged by the gluon. To perform the second integration, 
the following integral relation between Bessel functions has been used: 

n^in-'+c^f J^.r{Kr)dK = 0) ^ll^^ K^+x{cr) (A.15) 



n 



OO 

= 2, A = -1 j k{k^ + c^)"Vo(Kr) dK = Ko{cr) . 



Appendix B 



The parameter values of the 
Donnachie-LandshofF fit 

This appendix will give all the parameters of the Donnachie-Landshoff fit. This will enable the inter- 
ested reader to reproduce all our fits with the help of the publication [DL84]. 

Most of the parameters were taken straight from [DL84] and never modified. All the same, they 
are listed for the sakc^ of compk'tc^nc'ss in Tabk^ B.l. 



Object exchanged 


Parameter 


Numerical value 


Pomeron 


ao — 1 = e 


0.08 


a' 


0.25 GeV-2 




^/2 mbV4 GeV-V2 


'^dipole 


0.71 GeV^ 


2 Pomerons 




1.608 


Reggeon 


ao 


0.44 


a' 


0.93 GeV-2 


A 


7.8 


B 


2.1 


Gluon 


to 


0.3 GeV2 


ti 


0.3 GeV^ 


A 


2/tl 


T 


3 i 
- 2 II 


3 Particles 


A 


109.961 fmi/4 GeV^/^ 



Table B.l: All parameters of the Donnachie-Landshoff fit. 

"^dipoie square of a dipole mass which appears in the dipole form factor in Section 2.4.1. 

For double pomeron exchange, the constant £> • is given instead of D since this combination is 
independent of /?. D is the constant that ensures that the pomeron and double pomeron contributions 
cancel in the dip region, and the double pomeron amplitude is proportional to while the pomeron 



111 



112 



Appendix B. The parameter values of the Donnachie-Landshoff fit 



is proportional to 

ti is the gluon propagator cutoff, ie the value for —t below which the hyperbola 1/t is replaced 
by a parabola (see Figure 2.6) or some other cutoff function (see Section 3.2.2). Its value depends on 
the ciitoff function used and the coupling constant. The values for ti are displayed in Table B.2. The 
cutoff parabola used by Donnachie and Landshoff is parametrised by A and r. Its form is Xt{t + r). 

The last parameter A represents the strength of the coupling of the triple-particle exchanges. It is 
determined by fitting the triple gluon exchange differential cross section to an asymptotic power law. 
It depends on ti and the coupling constant. It is also listed in Table B.2. 

The values given for to, ti and A in Table B.l and in the first line of Table B.2 apply to the 
unmodified Donnachie-Landshoff fit which was used in nearly all calculations. The others were used 
when different cutoff functions for the gluon propagator were chosen and for a running coupling. 
Finally, some of them apply to calculations in which we changed the coupling constant of the DL fit. 
They were done before we decided to use the DL fit with the unchanged coupling constant as a fixed 
framework for the other odderon contributions. All the same the relevant parameters are given in 
Table B.2. 



Cutoff 


a.s 


to [GeV2] 


ti [GeV2] 


A [fmi/4 GeV^/4] 


Parabola 


0.3 


0.3 


0.3 


109.961 


Parabola 


0.3 


0.0003 


0.31 


109.589 


Parabola 


0.32 


0.0003 


0.31 


104.411 


Parabola 


0.4 


0.0003 


0.31 


88.321 


Parabola 


0.5 


0.0003 


0.31 


74.7105 


Flat 


0.4 


0.0003 


0.37 


88.5942 


Linear 


0.4 


0.0003 


0.22 


88.2452 


Parabola 


running 


0.0003 


0.44 


91.048 



Table B.2: The values of ti and A for different cutoff functions and coupling constants. 



Appendix C 

Derivatives of the conformal 
eigenfunctions 



When proving the conformal transformation properties of the G function, we used the explicit form of 
derivatives of the conformal eigenfunction . These derivatives will be calculated in this appendix. 
As a starting point, here is the explicit form of the conformal eigenfunction E'^"'^^ in position space: 

= (-^) . (c.i) 

\P\aPla) \P\aP2aJ 

a is the external coordinate of the pomeron state which does not enter into our calculations and which 
we will not write as an argument of E^'^'"''^ . 

The differential operators which occur in the G function are: A1A2, V1A2 and a single Ai. 
Expressed as derivatives with respect to the complex coordinates, they become: 16 81818282, 8 8I8282 
and 49i9J'. Since none of the exponents of coordinates in (C.I) is —1, we can completely separate the 
conjugated from the unconjugated coordinates.^ 

We will first calculate the gradient with respect to one coordinate of the eigenfunction . 



l + n , ■ 

Plap2aJ 



PlaP2aJ \ 2 ; {plaPLJ \PlaP2a (plaPL) 



1+r 



\ —H — / -1 \ / * 



P12 \ 


-2 — I-''' 


PlaPla) 






( ' 




\P\aPl 




1 




* 2 
Pla 



P\ap2aJ \ ^ J \PlaP2a 

= {l-n + 2iu) E^'^'-Hpi,P2) (C.2) 

PlaPl2 

Performing the differentiation with respect to the non-conjugated coordinate yields an analogous result. 
Here it is together with the gradients with respect to the other coordinate:^ 

VJi?('''")(Pi,P2) = (l+n + 2z:.)-^ (pi,p2) 

PlaPl2 

V2^("'"Hpi,P2) = -{l-n + 2w) £;('^'")(pi,p2) (C.3) 

P2aPl2 



^Otherwise, we might have to take relation (5.29) into account. 

^The symbol V* is meaningful in the complex, though not in the vectorial notation. See Section 5.2.3 for its definition. 
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Appendix C. Derivatives of the conformal eigenfunctions 



VIE^-'^'^ (p, , P2) = - (1 + n + 2 iv) (pi, p2) 

P2aPl2 

The derivatives with respect to get a minus sign because the result is expressed in pi^ instead of p2i • 

Having successfully derived the gradient, we c;an immediately give the Laplacians of a pomoron wave 
function E^'^'"-\ Since the differentiations with respect to (un)conjugated variables add only factors 
of variables of the same type, the two differentiations contained in a Laplacian are quite independent. 
Despite the pia and pi2 in the denominator of the prcfactor, the relation concerning the derivative of 
1/p (5.29) must not be used. Since iS^"^'") contains additional powers of both these variables, the total 
power is not —1, so their derivative with respect to the complex conjugated variable vanishes. The 
results are: 

A2^('^'")(pi,P2) 

The combination of a gradient and a Laplacian requires some more work. Here the gradient affects 
also the additional coordinate moduli in (C.4). We use the product rule to obtain: 

Vi A2£('''") (pi, P2) = Vi ((1 + 2 i,yf - n^) ^'''f E^'^'"^ {p^^p^) 

\p2ar\p12r 

= ((1 + 2 iuf - n') (2 (d*, , + I „ (1 - n + 2 ii.) ^) £;(^'") {p, , P2) • 

V V |P2aP|pi2|V \p2a\^\pl2\^ PlaPl2 J 

The derivative of the first factor becomes: 

^ \p2a?\pi2? \p2a?Pl2 ^ P*i2 \p2a? Pl2 VPi2 Pl2^ ) \P2a?\pl2? Ph \P2a?\pl2? PlaP*2 ' 

This term can easily be added to the derivative of the second factor, giving the constant prefactor 
(— 1 — n + 2iv). After factorising the prefactor of (C.4) and combining it with the new factor, we 
obtain the result: 

Vi A2i;('''") (pi , P2) = - (4 y^ + {n + If) {l-n + 2iv) . ""^f „ E^^'^) (pi, P2). (C.5) 

PlaPl2 \P2aV\pi2\^ 

The additional differentiation necessary to apply a double Laplacian to is again independent 

of the one just performed and can be done analogously. The moduli of pia and p2a then cancel out. 
The result is: 

Ai A2ii;('^'") (pi, P2) = (4 + (n + 1)2) (4 ^^ + {n- if) -L. (p,, p^) . (c.6) 



((1 + 2 



■ \2 2\ \P2a 
lu) — n - 



Plap|pi2p 



[[l + 2ivf-n') 

\p2aY\p12V 



E^-'''\pi,P2) 
E^^'^\pr,P2) 



(C.4) 
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